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t».  ABSTRACT 


>The  articulated  autorotor  decelerator  is  a  device  intended  to  provide  retardation  and  stability  to  air  delivered 
stores.  This  unique  device,  combining  aerodynamics,  dynamics,  and  gyroscopics,  shows  promise  of  having  a 
significant  potential  in  the  decelerator  field.  Limited  wind  tunnel  tests  using  small  models  indicate  that  this 
device  gives  the  largest  drag  coefficient,  based  on  solidity  factor  and  tip  speed  ratio  considerations,  of  any 
known  decelerator,  including  parachutes.!' )The  system,  is  composed  of  a  number  of  blades  attached  to  a  hub  at 
the  rear  of  the  store  body.  Each  blade  is  freeHo  spin  about  its  long  axis  by  means  of  a  bearing  at  the  point 
where  the  blade  is  attached  to  the  hub.  The  hub,  in  turn,  is  free  to  spin  about  an  axis  through  the  store  center 
line.  As  the  store  falls  through  the  air,  the  air  flow  over  the  blades  will  cause  them  to  spin  (or  autorotate)  about 
their  long  axes.  The  autorotation  is  accomplished  by  proper  shaping  of  the  blade  cross  section.  The  aero¬ 
dynamic  lift  and  drag  forces  acting  on  the  autorotating  blades  result  in  a  torque  about  the  store  center  line, 
thus  imparting  a  spin  rate  to  the  hub  (and  attached  blades).  The  aerodynamic  forces  acting  on  the  blades, 
because  of  the  blade  spin  and  hub  spin,  result  in  a  large  force  opposite  to  the  direction  of  the  store  flight 
(parallel  to  the  store  center  line),  which  is  the  system  drag.  The  blades  are  also  free  to  travel  through  a  fore  and 
aft  flap  angle.  This  third  degree  of  freedom  eliminates  large  bending  loads  on  the  blade  and  also  results  in  a 
self  correcting  action  for  stability.  The  method  considered  for  initial  blade  deployment  is  to  utilize  small 
turbine  yanes  mounted  on  the  hub  which  will  initiate  hub  rotation.  The  blades,  which  are  initially  lying  back 
parallel  to  the  store  center  line,  are  thrown  out  and  forward  (deployed)  by  the  centrifugal  force  due  to  hub 
rotation.  As  the  blades  pivot  out,  they  will  spin  and,  in  turn,  will  sustain  the  hub  spin.  The  blades  will  final. > 
deploy  to  an  equilibrium  position  at  a  near  right  angle  to  the  store  center  line.  This  method  of  deployment  has 
the  added  feature  of  providing  a  positive  means  of  initiating  blade  autorotation  during  deployment.  The  study 
includes  the  derivation  of  the  equations  describing  the  dynamic  motion  of  the  decelerator  system.  These 
equations  involve  the  inertial,  gyroscopic,  and  aerodynamic  effects,  The  aerodynamic  terms  are  written  so  that 
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13.  ABSTRACT  (Continued) 

the  aerodynamic  sectional  characteristics  of  the  autorotating 
blades,  obtained  from  wind  tunnel  tests,  can  be  utilized  directly. 
These  equations  are  programmed  to  be  solved  on  a  digital  com¬ 
puter  so  that  deployment  transients  and  steady  state  performance 
characteristics  can  be  investigated.  Experimental  data  obtained 
from  wind  tunnel  tests  with  a  representative  full  scale  rotor 
system  are  presented.  These  results  are  used  to  determine  a 
blade  interference  factor,  which,  when  incorporated  in  the  basic 
theory,  gives  excellent  correlation  between  analytical  and 
experimental  results.  The  qualitative  and  quantitative  effects  of 
various  system  parameters  on  the  steady  state  and  deployment 
transient  performance  of  the  system  are  investigated.  The 
results  of  this  study  afford  an  understanding  in  the  mechanics  of 
the  overall  system  by  investigating  the  relative  effects  of  the 
various  physical  terms  on  the  system  performance.  This  study 
als  provides  an  insight  into  the  performance  potential  of  the 
basic  configuration  and  represents  a  basic  analytical  approach 
for  future  applications  of  the  device. 
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DIGEST 


The  articulated  autorotor  decelerator  is  a  device  intended  to  provide  retardation  and 
stability  to  air-delivered  stores.  This  unique  device,  combining  aerodynamics,  dynamics,  and 
gyroscopics,  shows  promise  of  having  a  significant  potential  in  the  decelerator  field.  Limited  wind 
tunnel  tests  using  small  models  indicate  that  this  device  gives  the  largest  drag  coefficient,  based  on 
solidity  factor  and  tip  speed  ratio  considerations,  of  any  known  decelerator,  including  parachutes. 

The  system  is  composed  of  a  number  of  blades  attached  to  a  hub  at  the  rear  of  the  store 
body.  Each  blade  is  free  to  spin  about  its  long  axis  by  means  of  a  bearing  at  the  point  where  the 
blade  is  attached  to  the  hub.  The  hub,  in  turn,  is  free  to  spin  about  an  axis  through  the  store  center 
line.  As  the  store  falls  through  the  air,  the  air  flow  over  the  blades  will  cause  them  to  spin  (or 
autorotate)  about  their  long  axes.  This  autorotation  is  accomplished  by  proper  shaping  of  the  blade 
cross  section.  The  aerodynamic  lift  and  drag  forces  acting  on  the  autorotatirig  blades  result  in  a 
torque  about  the  store  center  line,  thus  imparting  a  spin  rate  to  the  hub  (and  attached  blades).  The 
aerodynamic  forces  acting  on  the  blades,  because  of  the  blade  spin  and  hub  spin,  result  in  a  large 
force  opposite  to  the  direction  of  the  store  flight  (parallel  to  the  store  center  line),  which  is  the 
system  drag.  The  blades  are  also  free  to  travel  through  a  fore  and  aft  flap  angle.  This  third  degree  of 
freedom  eliminates  large  bending  loads  on  the  blade  and  also  results  in  a  self-correcting  action  for 
stability. 


The  method  considered  for  initial  blade  deployment  is  to  utilize  small  turbine  vanes 
mounted  on  the  hub,  which  will  initiate  hub  rotation.  The  blades,  which  are  initially  lying  back 
parallel  to  the  store  center  line,  are  thrown  out  and  forward  (deployed)  by  the  centrifugal  force  due 
to  hub  rotation.  As  the  blades  pivot  out,  they  will  spin  and,  in  turn,  will  sustain  the  hub  spin.  The 
blades  will  finally  deploy  to  an  equilibrium  position  at  a  near  right  angle  to  the  store  center  line. 
This  method  of  deployment  has  the  added  feature  of  providing  a  positive  means  of  initiating  blade 
autorotation  during  deployment. 

The  study  includes  the  derivation  of  the  equations  describing  the  dynamic  motion  of  the 
decelerator  system.  These  equations  involve  the  inertial,  gyroscopic,  and  aerodynamic  effects.  The 
aerodynamic  terms  are  written  so  that  the  aerodynamic  sectional  characteristics  of  the  autorotating 
blades,  obtained  from  wind  tunnel  tests,  can  be  utilized  directly.  These  equations  are  programmed 
to  be  solved  on  a  digital  computer  so  that  deployment  transients  and  steady  state  performance 
characteristics  can  be  investigated. 

Experimental  data  obtained  from  wind  tunnel  tests  with  a  representative  full  scale  rotor 
system  are  presented.  These  results  are  used  to  determine  a  blade  interference  factor,  which,  when 
incorporated  in  the  basic  theory,  gives  excellent  correlation  between  analytical  and  experimental 
results. 


The  qualitative  and  quantitative  effects  of  various  system  parameters  on  the  steady  state 
and  deployment  transient  performance  of  the  system  are  investigated.  The  results  of  this  study 
afford  an  understanding  of  the  mechanics  of  the  overall  system  by  investigating  the  relative  effects 
of  the  various  physical  terms  on  the  system  performance.  This  study  also  provides  an  insight  into 
the  performance  potential  of  the  basic  configuration  and  represents  a  basic  analytical  approach  for 
future  applications  of  the  device. 
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AN  ANALYSIS  OF  THE  DYNAMIC  AND  AERODYNAMIC  PERFORMANCE 
OF  A  SELF-DEPLOYING  ARTICULATED  AUTOROTOR  DECELERATOR 


I.  INTRODUCTION. 

The  articulated  autorotor  decelerator  is  the  invention  of  James  Brunk  and  Abraham 
Flatau.1  Limited  wind  tunnel  tests  at  Edgewood  Arsenal  using  small  models  indicated  an  extremely 
large  drag  coefficient  for  a  relatively  low  value  of  solidity  factor  and  tip  speed  ratio  when  compared 
with  other  decelerator  devices  which  utilize  some  form  of  autorotation  motion  in  their  mode  of 
operation.  These  devices  would  include:  vortex  ring  parachute,  rotochute,  flexirotor,  and  rotafoil 
parachute. 

f 

Limited  wind  tunnel  data  have  indicated  the  basic  potential  of  this  device  and  verify  the 
original  reasoning  that,  led  to  the  articulated  autorotor  decelerator  concept  in  that  a  large  system 
.  drag  coefficient  would  result  because  of  the  combination  of  both  a  large  lift  and  a  large  drag 
coefficient  characteristic  of  the  basic  autorotor  blade. 

This  report  presents  the  first  known  attempt  to  analyze  the  qualitative  and  quantitative 
performance  characteristics  of  the  articulated  autorotor  from  a  detailed  dynamic  and  aerodynamic 
viewpoint.  , 


The  articulated  autorotor  decelerator,  which  appears  in  figure  1 ,  is  particularly  applicable 
for  retardation  of  air-delivered  stores  because  of  the  stability  it  would  provide  in  addition  to  the 
basic  decelerator  function.  This  analysis  is  concerned  mainly  with  utilizing  the  device  for  the 
air-delivered  store  application.  Thus’  the  selection  of  the  general  analytical  model  configuration  and 
associated  dimensions  is  swayed  in  this  direction.  The  physical  system  analyzed  is  intended  to 
include  all  the  terms  basic  to  the  mechanics  of  the  system,  yet  be  as  simple  as  possible  to  facilitate 
an  understanding  of  the  fundamental  aspects  of  the  device. 

This  analysis  involves  the  dynamic  motion  and  aerodynamic  performance  of  the 
articulated  autorotor  decelerator  system  in  a  constant  and  axisymmetric  free  stream  velocity,  it 
does  not  include  trajectory  aspects  nor  consider  an  accelerating  body  system.  Nonetheless,  the 
performance  analysis  is  illustrative  of  the  behavior  of  the  system,  and  the  general  approach  and  the 
equations  thus  derived  are  valid  for  future  analyses  of  other  applications  of  the  basic  system. 

II.  GENERAL  DESCRIPTION  OF  THE  SYSTEM. 

The  articulated  autorotor  decelerator  system  as  applied  to  an  air-delivered  store 
application  is  illustrated  in  figure  1  and  functions  in  the  following  manner.  The  system  is  composed 
of  a  number  of  blades  (four  are  shown  in  the  illustration),  which  are  attached  to  a  hub  at  the  rear 
end  of  the  store  body.  Each  blade  is  free  to  spin  about  its  long  axis  by  means  of  a  bearing  at  the 
point  where  the  blade  is  attached  to  the  hub.  The  hub,  with  the  blades  attached,  is  free  to  spin 
about  an  axis  coincident  with  the  store  center  line  by  means  of  a  bearing  between  the  hub  arid  the 
store  body. 

As  the  store  falls  through  the  air,  the  air  flow  over  the  blades  will  cause  them  to  spin  or 
autorctate  about  their  long  axes.  This  spin  rate  is  shown  as  co  in  figure  1.  This  autorotation  is 
accomplished  by  proper  shaping  of  the  blade  cross  section.  The  aerodynamic  Eft  and  drag  acting  on 


•Flatau,  A.,  and  Brunk,  J.  E.  Free  Spinning  Articulated  Rotor.  US  Patent  3,291,418. 13  December  1966. 
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Figure  1 .  General  Configuration  of  Articulated  Autorotor  Decelerator  System 

the  blades  caused  by  their  autorotation  will  result  in  a  torque  about  the  store  center  line,  thus 
imparting  a  spin  rate  to  the  hub  (and  the  attached  blades).  This  hub  spin  rate  is  shown  as  SI  in  Figure 
1 .  The  aerodynamic  forces  acting  on  the  blades  caused  by  the  combination  of  blade  spin  and  hub 
spin  result  in  a  large  force  in  the  direction  opposite  to  the  store  flight  direction  (parallel  to  the  store 
center  line),  which  is  the  system  drag. 

The  blades  are  also  free  to  pivot  through  a  flap  angle,  shown  as  6  in  figure  1 .  This  third 
degree  of  freedom  eliminates  large  bending  loads  on  the  blade  and  also  results  in  a  self-correcting 
action  for  stability.  The  autorotor  blade  used  in  this  analysis  is  rectangular  in  plan  form  and  has  a 
cross-sectional  shape  similar  to  a  flattened  S.  This  shape  is  usually  referred  to  as  an  S-section 
autorotor. 

A  method  for  initial  deployment  of  the  blades  is  shown  in  figure  2.  Prior  to  deployment, 
the  blades  are  lying  back  parallel  to  the  store  center  line.  A  hub  spin  rate  is  achieved  by  some  means 
such  as  small  external  turbine  va..es  mounted  to  the  hub.  When  released,  the  blades  are  thrown  out 
and  forward  by  the  centrifugal  force  resulting  from  the  hub  spin  rate.  As  the  blades  pivot  out  to 
larger  values  of  6 ,  they  will  spin  about  their  long  axes;  the  resulting  aerodynamic  forces  on  the 
blades  will  sustain  the  hub  spin  rate.  The  blades  will  finally  deploy  to  an  equilibrium  position  at  a 
near  right  angle  to  the  store  center  line.  Thus,  the  flap  angle  0  will  hereafter  be  referred  to  as  the 
deployment  angle. 

This  aft-to-forward  deployment  method  not  only  minimizes  blade  bending  loads  that 
would  occur  in  a  forward-to-aft  deployment  when  the  blades  impacted  the  necessary  stop,  but  it 
also  has  the  added  feature  of  providing  a  positive  means  of  initiating  blade  autorotation  during 
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1.  Blade  position  prior  to  deployment: 

•  •  •  .  £1  —  0 
8=6  =  0=  0  to  =  a>  =  0  0  _  _ 

it  =  it  Q 

2.  Blade  position  during  deployment  transient:  all  terms  and  their  derivatives  have  values 

3.  Blade  position  at  steady  state  condition: 

0=0*=  0  d>  =  0  ft  =  0 

e=0ss  w=uss  ft  =  nss 


deployment.  It  is  this  automatic  initiation  of  blade  autorotation  along  with  the  unassisted  deploy¬ 
ment  motion  of  the  blades  when  released  (assuming  a  sufficient  hub  spin  rate)  that  are  the  reasons 
behind  the  term  “self-deploying”  in  the  system  description. 

III.  GENERAL  METHOD  OF  ANALYSIS. 

The  equations  describing  the  motion  of  the  system  are  developed  (section  IV)  for  the  case 
of  an  inertially  fixed  (i.e.,  nonaccelerating)  body.  The  aerodynamic  effects  are  treated  for  the  case 
of  constant  and  axisymmetric  flow  (i.e.,  parallel  to  the  store  center  line).  These  conditions  represent 
the  physical  situation  of  the  system  when  installed  in  a  wind  tunnel  (fixed  mount),  and  therefore 
the  results  are  for  a  condition  that  can  be  experimentally  verified. 

The  basic  blade  element  approach,  used  in  propeller  theory,  is  utilized  to  determine  the 
aerodynamic  load  distribution  over  the  blade,  covered  in  section  V.  This  approach  allows  the 
experimentally  determined  aerodynamic  sectional  characteristics  of  the  autorotor  blades  to  be  used 
directly  and  also  permits  a  detailed  understanding  of  the  aerodynamic  effects  over  the  blade. 

These  equations  were  programmed  to  be  solved  on  a  digital  computer,  as  described  in 
section  VI,  not  only  to  facilitate  the  solution  of  the  several  complex  equations,  but  also  to  allow  a 
detailed  insight  into  the  results  because  a  large  number  of  terms  could  be  computed  and 
investigated. 

Steady  state  system  performance  (for  a  constant  free  stream  velocity)  was  computed  and 
compared  with  experimental  results  obtained  in  a  wind  tunnel  test  using  a  full  scale  system 
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described  in  section  VII.  Comparison  between  the  initial  analytical  predictions  and  experimental 
results  showed  fair  correlation,  but  the  differences  were  unsatisfactorily  large. 

A  parametric  study  was  conducted  in  which  each  physical  property  and  condition  were 
varied  about  a  nominal  condition  to  determine  if  the  sensitivity  of  certain  terms  would  explain  the 
poor  correlation  between  the  theoretical  and  experimental  results  (section  VIII).  None  of  the 
parameters  studied  influenced  the  performance  either  qualitatively  or  quantitatively  toward  the 
experimental  results.  Attention  was  then  centered  on  the  aerodynamic  load  distribution  over  the 
blades.  From  an  initial  inspection  of  local  flow  directional  effects  and  an  empirical  analysis  using 
various  degrees  of  this  effect,  it  was  determined  that  mutual  blade  interference  effects  on  the  air 
flow  direction  over  the  blades  (i.e.,  local  angle  of  attack)  would  bring  all  the  analytically  predicted 
terms  more  in  agreement  with  their  experimental  values.  Differences  between  the  analytical  and 
experimental  results  were  effectively  eliminated  in  all  terms  by  adding  a  blade  interference  effect  to 
the  theoretical  analysis. 

With  the  theoretical  approach  verified,  a  parameter  study  was  conducted  to  indicate  the 
effects  of  the  various  terms  on  the  system  performance.  Also  presented  is  a  detailed  examination  of 
the  steady  state  values  of  specific  terms  of  interest.  These  results  are  included  in  section  IX. 

Finally,  the  transient  behavior  of  the  system  during  deployment  is  presented  in  sec¬ 
tion  X.  .... 

IV.  DERIVATION  OF  THE  EQUATIONS  OF  MOTION. 

A.  Formulation  of  the  Basic  Equations. 

The  basic  physical  model  representing  the  articulated  autorotor  decelerator  system  is 
shown  in  figure  3.  The  equations  describing  the  motion  of  this  system. will  be  developed  using  the 
modified  Euler’s  equations  approach.2 

Referring  to  figure  3,  X,  Y,  Z*  is  an  inertial  axes  system  fixed  to  ifle  store  body  with  the 
origin  ( O )  at  the  intersection  of  the  blade  spin  axis  and  hub  spin  axis,  /here  X  is  directed  along  the 
store  center  line,  and  Z  is  directed  downward.  The  body  axes  system  x,  y,  z  is  fixed  to  the  blade  and 
oriented  along  the  principal  axes  of  the  blade,  with  the  x-axis  oriented  along  the  blade. 

Let  «'  be  the  absolute  angular  velocity  of  the  blade  with  respect  to  the  X,  Y,  Z  system. 
Let  w  be  the  absolute  angular  velocity  of  the  x-,  y-,  z-axes  system  (fixed  to  the  blade)  with  respect 
to  the  X,  Y,  Z  system.  The  x-axis  is  oriented  along  the  blade  as  shown,  and  the  x-,  y-,  z- axes  system 
could  be  rotating  about  the  x-axis. 


Angular  momentum  =  H  =  Id' 


Expanding  (4-1): 


HX=WX 

Hy  =lyU>‘y 

Hz  =  Izu'z 


(4-1) 


w 

(4-2) 


2Greenwood,  D.  F.  Principals  of  Dynamics.  Prentice-Hall  Inc.,  Englewood  Cliffs,  New  Jersey. 
*The  mathematical  terms  are  defined  in  glossary  at  the  end  of  the  text. 
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4  =/i 

4=4=4 

Although  the  instantaneous  moment  of  inertia  of  the  blade  about  the  point  <?  would 
depend  on  whether  the  blade  were  oriented  with  the  flat  face  or  edge  on  relative  to  the  axis 
considered,  the  fact  that  the  blade  has  a  high  spin  rate  allows  a  mean  value  for  the  moment  of 
inertia  to  be  used.  This  value,  which  is  between  the  two  extreme  values  mentioned,  is  still  close 
enough  to  these  in  magnitude  to  allow  the  blade  to  be  considered  as  a  symmetrical  body. 

Substituting  these  into  (4-2)  through  (4-4): 

fix  =4«4 

Hy  =I0u'y  (4-3) 

*4=4  «*, 
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Let  the  blade  rotate  about  its  axis  of  symmetry  (x)  with  angular  velocity  s  (with  respect 
to  the  x-,  y-,  z-axes  system).  ' 


u’x=ux  +s 

CO  ‘ y  ~  Uy 

co',  =  CO. 


Substituting  these  results  into  the  angular  momentum  expression  (4-3), 

Hx  =  IX  (co*  +5) 

Hy  =lQO}y 

Hz  =Iquz 

Differentiating  equations  (4-5)  with  respect  to  time: 

Hx‘Ii(&x+*)‘ 

Hy  I0  ^  y 

Hz  =/0<i, 

where  the  dot  over  the  term  indicates  differentiation  with  respect  to  time: 


^dn 


and  (  ) = 


,♦  dH  ) 


The  time  rate  of  change  of  angular  momentum  referred  to  an  inertial  coordinate  system 
in  terms  of  a  rotating  coordinate  system  exp/e  sed  in  vector  notation  is 

•  • 

II  =  (H)r  +  co  X  H  (4-7) 

• 

where  ( H)r  is  the  time  rate  of  change  of  the  absolute  angular  momentum  referred  to  the  rotating 
coordinate  system,  H  is  the  absolute  angular  momentum  referred  to  the  rotating  coordinate  system, 
and  co  is  the  angular  velocity  of  the  rotating  coordinate  system  relative  to  the  inertial  coordinate 
system. 


Expanding  the  terms  of  (4-7), 

(//),  =Hxl+HyJ+Hzk 
co  =  coJ+  Uyj  +  uzk 

H  =HxT+HyT+Hzk 

where  Z  T kart  unit  vectors  for  the  x,  y,  z  system. 


(4-8) 

(4-9) 

(4-10) 


Evaluating  the  first  term  on  the  right-hand  side  of  (4-7),  and  substituting  (4-6)  into  (4-8), 


(H).  =  7j  (cox  +  s)i  +  I0uyj  +  I0uzk 


(4-11) 


Evaluating  the  second  term  on  the  right-hand  side  of  (4-7), 

7  7  k  I 


co  X  H  = 


H,  Hy 


co  X  H -  (c oyHz  -  u>zHy)i  +  (<0*77*  -  c oxHz)j  +  (c oxHy  -  uyHx)k 


(4-12) 
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Substituting  the  expressions  from  (4-1 1 )  and  (4-1 2)  into  (4-7)  gives: 

H  =  [/,  (oix  +  s )  +  uyHz  -  o ozHy  1/  +  (J0&y  +  iozHy  -  co * II, )/  +  (/0 oo2  +  o oxHy  -  c oyHz  )k  (4-13) 

Expanding  (4-13)  gives: 

H=  (^co^+Zji  +  uyI0uz  -  iozI0cjy)i 

+  [I0u>y  +  coz/1(co*  +  5)  -  ioxI0uz]j  (4-14) 

+  [/Ocoz  +ojxI0ioy  -  coyIl(ux  +  $)]£' 

The  time  rate  of  change  of  angular  momentum  is  equal  to  the  net  external  moment;  in 
vector  notation, 

M  =  H  (4-15) 

Expanding  (4-1 5)  and  substituting  in  the  expressions  from  (4-14): 

Mx  =/i^x  +  V 

My  =  I0U)y  -  (I0  -  Il)lOxOJ2  +  (4-16) 

Mz  -  /Ocoz  +(/0  -  ll)cjxuy  -  Ixuys 
B.  Conversion  of  Euler  Angles. 

_  _  _The  Euler  angles  (7,  5,  0),  which  describe  the  orientation  of  the  bladewith  lespect  to  the 
X-,  Y-,  Z- axes  system,  are  chosen  so  that  the  K-axis  always  remains  in  the  YZ  plane  (i.e..  7  is 
measured  in  the  YZ  plane),  as  shown  in  figure  3.  , 


Expressing  the  angular  rates  in  terms  of  these  Euler  angles, 

s  =  0 
ojx  ~  -7  sin  6 
uy=b 
03  z  =  7  cos  8 

Differentiating  (4-1 8)  with  respect  to  time: 

w*  =  78  cos  5  -  7  sin  5 

•  -  V 
<0^-6 

coz  =  76  sin  5+V  cos  6 

Substituting  (4-1 7),  (4-1 8),  and  (4-1 9)  into  (4-1 6): 

Mx  =  -  1 1 76  cos  5  -  /j  7  sin  8  +  /j 0 

My  =  IqS  +(Iq-  /j  ) 72  sin  8  cos  8  +  /j  70  cos  8 

Mz  =  -70  78  sin  6  +/0 7  cos  8  -  (70  -  /j  )yS  sin  8  -  /j06 

Expressing  the  moment  terms  about  the  Euler  axes  system  as  defined  in  figure  3: 

My  ~~MX  sin  8  +■  Mz  cos  8 

M&  =My 

M*=MX 


(4-17) 


(4-18) 


(4-19) 


(4-20) 


(4-21) 


Substituting  (4-20)  into  (4-2 1 )  results  in: 

My  =  V(/0  cos2  6  +  /[  sin2  6 )  -  2(/0  -  / j  )yS  sin  6  cos  8  -  /j  05  cos  6  -  /j  <5  sin  5 

4/s  =  /0*S  +  (/0  -  /j  )y2  sin  8  cos  6  +  /t  -y*0  cos  6  (4-22) 

4L,  =  / j  78  cos  8  -  /|  7  sin  8  + 1  j  0 

Note  that  the  equations  do  not  contain  the  terms  7  or  <t>.  nor  any  functions  of  them,  but  include 
only  the  first  and  second  derivatives  (rates  and  accelerations)  of  these  terms. 

C.  Changing  Descriptive  Terminology  Symbols. 

The  basic  symbols  used  thus  far  to  evolve  these  equations  are  now  changed  to  different 
symbols,  which  have  more  meaning  to  the  physical  system  being  analyzed.  These  changes  are  shown 
in  figure  4. 

The  symbol  denoting  rotational  rate  of  the  blade  about  the  store  center  line  7  is  replaced 
by  the  symbol  £2,  which  is  common  terminology  for  hub  rotation  rate.  All  other  terms  are  altered 
accordingly: 


7  =  £2  7  =  SI  My  -  Mn 

The  symbol  denoting  rotational  rate  of  the  blade  about  its  long  axis  <t>  is  replaced  by  the 
symbol  oj,  which  is  common  terminology  for  autorotor  spin  rate,  resulting  in: 

•  •  •  • 

0  =  w  0  =  0)  4/0  =MCJ 

The  angular  symbol  8 ,  which  describes  the  angle  of  the  blade  measured  from  a  plane 
perpendicular  to  the  store  center  line  (hub  spin  axis),  is  replaced  by  its  complementary  angle  0  so 
that: 


sin  8  =  cos  0 

•  • 

6  =-e 


cos  8  -  sin  0 

•  •  •• 

8  =  -0 


8  =  90°  -  O 
Ms=-Mg 


The  angle  0  is  then  termed  the  “deployment  angle,”  which  has  a  more  useful  physical  meaning  in 
this  application. 


Substituting  these  new  terms  in  equations  (4-22)  gives: 

4/n  =  n(/0  sin2  0  +  7j  cos2  0 )  +  2(/0  -•  )Sld  sin  0  cos 0  +  lv a>0  sin  0  -  Ixu  cos 0 
-Me  =  -IqO  +  (/q  -  1 1  )£22  sin  0  cos  0  +  I^Slco  sin  0 
4/w  =  /j SIO  sin  0  -  flSl  cos 0  +  /j  w 

Solving  for  the  angular  accelerations: 


(4-23) 


4/. 


cj  -  ■ 


-  SIO  sin  0  +£2  cosO 


4/, 


0  =—+  (/0  -  /j  )Sl2  sin  9  cos  0  +  /j  £2«  sin  0 


(4-24) 


SI  = 


(7osin20  +  /j  cos20)  [/W«  '  2(1°  '  !^l°  sin0  cos0  ‘  7iw*  sin  0  +  ^ <0  cos0] 
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Figure  4.  Basic  System  Terminology 


D.  Formulation  of  the  Aerodynamic  lvfoments. 

The  external  moments  M w ,  Mg ,  and  Mn  in  expressions  (4-24)  are  assumed  to  result 
solely  from  aerodynamic  forces.  Blade  weight  and  bearing  friction  for  both  the  blade  and  hub  are 
not  considered. 

The  aerodynamic  load  distribution  acting  over  the  blade  is  evaluated  by  the  standard 
blade  element  propeller  theory  method.3  In  this  method,  the  magnitude  and  directional  sense  of  the 
aerodynamic  force  is  evaluated  at  distinct  locations,  these  being  the  center  of  area  elements  located 
at  incremental  intervals  along  the  blade.  The  forces  are  a  function  of  the  basic  aerodynamic 
characteristics  of  the  blade  shape  and  the  local  flow  conditions  at  each  el  .-^ent. 

In  addition  to  the  determination  of  the  aerodynamic  load  distribution,  the  forces 
evaluated  at  each  element  are  summed  up,  resulting  in  the  net  force  and,  consequently,  the  net 
moment  acting  on  the  blade. 

In  this  analysis,  axisymmetric  flow  relative  to_the  store  body  is  assumed  where  the  free 
stream  velocity  V„  is  parallel  to  the  body  center  line  (A'-axis,  which  is  also  the  hub  spin  axis, ft). 
Consider  the  aerodynamic  forces  acting  on  an  element  of  the  blade  located  at  a  distance  r  from  the 


3Von  Mises,  R.  Theory  of  Flight.  Dover  Publications,  Inc.,  New  York,  New  York.  1959. 


point  where  the  blade  int<  rsects  the  body  center  line  axis  X ,  as  shown  in  5.  The  equations 
describing  the  motion  of  the  blade  derived  thus  far  do  not  depend  on  the  angular  position  of  the 
blade  as  it  rotates  about  the  store  center  line,  indicated  as  7  in  figure  4.  Thus!  to  simplify  the 
geometry  in  the  following  analysis,  the  blade  is  considered  to  be  oriented  in  the  XZ  plane. 

The  relative  velocity  at  this  element  is  V,  which  is  the  vector  sum  of  the  free  stream 
velocity  and  the  tangential  velocity  of  the  element  caused  by  the  hub  spin. 

V  =  sTvJ  +(rQ. i  6)2  (4-25) 

Thejelative  angle  of  attack  a ,  which  is  defined  here  as  the  angle  between  the  resultant 
local  velocity  V  and  the  free  stream  velocity  ,  is  as  follows: 


Figure  5.  Aerodynamic  Conditions  at  a  Typical  Blade  Element 
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The  velocity  component  at  the  element  resulting  from  the  blade  deployment  rate,  rO ,  is 
not  included  in  the  computation  of  the  relative  velocity  because  it  is  small,  relative  to  the  other 
velocity  terms;  its  inclusion  is  not  warranted  in  light  of  the  increased  complexity  that  would  be 
added  to  the  already  complicated  terms. 

Let  xw,  yw,  zw  be  a  wind  ajces  coordinate  system  centered  at  the  blade  element 
considered  (/),  where  xw  is  parallel  to  V,  zw  is  perpendicular  to  both  V  and  go,  and  yw  is 
perpendicular  to  the  plane  defined  by  xwzw.  The  aerodynamic  force  acting  at  this  element  can  be 
expressed  as  a  drag  force  component  (D),  acting  along  xw ,  a  lift  force  component  ( L )  acting  along 
zw ,  and  a  side  force  component  ( SF)  acting  along  yw .  The  directional  sense  of  the  aerodynamic 
forces  are  indicated  in  figure  5.  It  should  be  noted  that  the  wind  axes  system  used  is  a  “left-hand” 
system,  as  opposed  to  the  commonly  used  “right-hand  rule”  concerning  the  positive  direction  of 
orthogonal  vectors.  This  was  done  to  facilitate  the  derivation  of  the  angular  relations  for  the 
physical  system  considered. 

The  aerodynamic  forces  acting  on  the  blade  in  the  autorotor  decelerator  system  will  be 
expressed  in  such  a  way  that  the  autorotor  blade  sectional  aerodynamic  coefficients  can  be  utilized 
directly  in  the  same  form  in  which  they  have  been  determined  from  wind  tunnel  tests  of  the 
particular  autorotor  blade  sectional  configuration. 

The  aerodynamic  forces  are  expressed  in  the  standard  coefficient  form  where: 


D  =  Qj  qS 

L  =  CLqS 

(4-27) 

SF  =  CyqS 

7j  =  dynamic  pressure  = 

(4-28) 

S  =  planform  area  =  CA; 

(4-29) 

The  dynamic  pressure  used  is  the  local  value,  evaluated  at  the  center  of  the  element;  the 
area  of  the  element  CAr  is  used  as  the  reference  area,  as  shown  in  figure  6. 

The  aerodynamic  coefficients  CL,?Dt  and  CY  are  obtained  from  wind  tunnel  tests  on  a 
rotor  blade  of  identical  cross-sectional  shape.4 

It  should  be  noted  that  once  the  blade  spin  «  has  been  initiated,  the  aerodynamics  are 
not  a  function  of  the  angle  of  the  air  flow  relative  to  the  blade  in  a  plane  normal  to  the  blade  spin 

axis.  The  aerodynamic  coefficient  values,  however,  are  a  function  of  the  blade  tip  speed  ratio  -jp-  at 

the  element,  and  ’>  v  angle  < p  of  the  blade  to  the  local  velocity  vector  V,  as  shown  in  figure  6. 

A  mot  deiai  .  ascription  of  the  form  of  these  coefficients,  as  well  as  the  method  by 
which  they  are  obtain*'  included  in  section  V. 

The  aerodynamic  characteristics  of  the  autorotor  blade,  as  determined  from  wind  tunnel 
tests,  are  expressed  as  a  function  of  the  yaw  angle  ,  as  illustrated  in  figure  6. 


4Flatau,  A.  CRDL  Technical  Memorandum  15-49.  Preliminary  Investigation  of  Modified  Autorotating  Firmer.1 
Cylinders.  17  August  1961 .  UNCLASSIFIED  Memorandum. 
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Figure  6.  Details  of  Blade  Element  Terminology 


__  For  the  autorotor  decelerator  system,  the  true  yaw  angle  between  the  relative  velocity 
vector  V  and  the  blade  spin  axis  cj  is  termed  0,  as  shown  in  figure  6.  The  angle  0  is  the  complement 
of  the  angl j  'P,  ,  ' .  ■ 

i  ■  t 

(3  =  90°.  -  tfr  (4-30) 

From  figure  7  it  can  be  seen  that  the  angle  0  is  a  function  of  both.  6  and  a: 

0  =  cos-1  (cos  6  cos  a)  (4-31) 

Figure  7  also  illustrates  the  details  of  the  derivation  of  the  angular  relations  between  the 
wind  axes  and  the  body  axes  systems  that  are  required  to  resolye  the  aerodynamic  forces  in  the 
proper  form  for  inclusion  into  the  equations  of  motiop  of  the  system.  Let  the  wind  axes  system 
Xw,  Yw,  Zw  be  initially  oriented  to  coincide  with  the  X-,  Y-,  Z-axes,  with  XM,  directed  along  Y,  Yw 
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plane  Z^0zw  is  perpendicular  to  V; 
.*.  intersection  of  planes  Q0zw  and 

2w0zw  deffne  3  line  zw  which  is 
perpendicular  to  both  w  and  V 
where: 

tan£  =  tan (90°  -  0)cos(9O°  -  a) 
or:  tans  =  ctn0  sina 


Figure  7.  Derivation  of  True  Blade  Yaw  Angie  and  Wind  Axea  Angular  Relariona 
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directed  along  Z,  and  Zw  directed  along  the  negative  A'-axis.  Rotate  the  X„,  Yw,  Zw  system  about 
the  Tw-axis  through  the  angle  (90°  -"a).  Now  Xw  is  parallel  to  V,  and  Z„  is  perpendicular  to  V. 
Next,  rotate  the  axes  X'w,  Y'w,  Z^  about  the  axis  X'w  through  the  angle  £,  where: 

tan  £  =  tan(90°  -  0)  cos(90°  -  a) 

or:  £  =  tan-1(cfn0  sine?)  (4-32) 

Now  zw  is  perpendicular  to  both  V  and  u  (i.e.,  the  long  axis  of  the  blade).  The  yw  -axis  is  likewise 
perpendicular_to  V,  but  lies  in  the  plane  of  V  and  to .  The  drag  force  (£>)  is  directed  along  xw ,  which 
is  parallel  to  V ;  the  lift  force  (I)  is  directed  along  zw,  which  is  perpendicular  to  the  plane  of  0.  The 
side  force  ( SF)  is  directed  along  yw ,  which  is  perpendicular  to  V  and  lies  in  the  plane  of  0.  The 
aerodynamic  forces  relative  to  the  blade  and  to  the  relative  velocity  are  now  identical  in  form  to  the 
wind  tunnel  data;  the  relations  between  these  forces  and  the  autorotor  decelerator  system  body 
axes  are  defined. 

_  _  The  aerodynamic  forces  acting  along  the  xw-,yw-,  zw-axes  will  not  be  resolved  along  the 
X-,  Y-,  Z- axes,  figures  5  and  7. 

Considering  figure  8,  which  contains  a  view  at  the  blade  element  looking  along  theZ-axis 
toward  the  inertial  axes  origin,  and  figure  9,  which  contains  a  view  looking  along  the  T-axis  toward 
the  origin,  it  can  be  seen  that: 

Fy  =  L  cos  £  sin  a  +  D  cos  a  -  SF  sin  £  sin  a 

Fy  =  L  cos  £  cos  a-  D  sin  a  -  SF  sin  £  cos  "a  (4-33) 

Fy  -  L  sin  £  +  SF  cos  £ 

Now  consider  figure  9,  which  contains  a  view  of  the  blade  element  looking  along  the 
T-axis  toward  the  element.  The  forces_of  (4-33)  can  be  resolved  into  a  force  component  that  is 
normal  to  the  blade  spin  axis  in  il'a-  X-,  Z- plane. 

Fn  -  Fy  sin  0  +  Fy  cos  0  *  (4-34) 

Substituting  the  expressions  from  (4-33)  into  (4-34)  results  in: 

Fn  =  ( L  cos  £  sin  cT+ Z)  cos  a  -SF  sin  £  Sin'S)  sin  0  +  (Z,  sin  £  +  SF  cos  £)cos  0  (4-35) 

_Jhus,  the  component  of  the  aerodynamic  force  acting  normal  to  the  blade  in  the  0-plane 

(i.e.,  the  X-,  Z- plane)  for  a  particular  blade  element  is: 

| F9  =  FN  =  L{ cos  £  sin  a  sin  0  +  sin  £  cos  0)  +  Z)  cos  a  sin  0  +  SF( cos  £  cos  0  -  sin  £  sina  sin  0)(4-36)  ; 

The  net  normal  force  acting  on  the  blade  in  the  0 -plane  is  the  summation  of  the  normal 
forces  evaluated  at  all  the  elements  along  the  entire  length  of  the  blade.  Using  the  aerodynamic 
expressions  of  (4-27)  through  (4-29)  in  (4-36),  and  integrating  over  the  length  of  the  blade  results  in 
the  following  expression  for  the  net  normal  force : 


Fe  =f^\ Cx^t*C(cos  £  sin  a  sin  0  +  sin  £  cos  0)  +  CD  C  cos  a  sin  0 

r0 

0v2  _  I 

+  CYJ^2~ C(cos£  cos0  -  sin  £  sina  sin0)  j  dr 


(4-37) 
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Figure  8.  Aerodynamic  Forces  Contributing  to  Hub  Spin  Moment  and  System  Drag 


to 


z 


Note:  View  along  Y  axis  toward  X,Y,Z  origin 


Figure  9.  Aerodynamic  Forces  Contributing  to  Biade  Deployment  Moment 
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Where  r0  indicates  the  location  of  the  initial  blade  element,  nearest  the  root  of  the  blade, 
and  Tf  indicates  the  final  outermost  element,  near  the  blade  tip. 


Thus,  the  net  aerodynamic  moment  acting  on  the  blade,  about  the  blade/hub  pivot  point, 
in  the  0 -plane  is: 


-f'\ 


CL^~. ~(cos  £  sin  a  sin  0  +  sin  £  cos  0 )  +  CD 


pV2  _ 
tL-y-  C  cos  a  sin  0 


+  C 


pV 2 


C  (cos  £  cos  6  -  sin  £  sin  a  sin  0)  rdr 


(4-38) 


The  component  of  the  aerodynamic  force  at  each  blade  element  acting  normal  to  the 
blade  spin  axis  and  in  the  ft  -plane  (i.e.,  the  plane  normal  to  the  Ma  vector  in  figure  4)  is  shown 
from  figure  8  to  be: 


Fn  =  Fy  =  I.  cos  £  cos  a  -  D  sin  a  -  SF  sin  £  cos  a 


(4-39) 


The  net  force  acting  to  give  a  “torque”  to  the  blade,  about  the  hub  spin  axis,  is  obtained 
by  the  summation  of  these  load  forces  evaluated  at  all  the  elements  along  the  entire  length  of  the 
blade.  Substituting  the  aerodynamic  terms  of  (4-27)  through  (4-29)  into  the  expression  (4-39),  and 
integrating  over  the  length  of  the  blade,  results  in  the  following  expression  for  the  net  hub  “torque” 
force: 


'.-/'I 


C  cos  £  cos  a  -  C  sin  a  -  Cy^y-  C  sin  £  cos  a  I 


<4-40) 


Thus,  the  aerodynamic  moment  acting  on  the  blade  about  the  blade/hub  pivot  point,  in 
the  ft -plane  is: 


Ma  =  f  j  Qyp C cos  £  cos  a  -  CD^~Cs\n  a  -  Cy^-C sin  £  cos  a  J  rdr 


(4-41) 


element  is: 


The  aerodynamic  moment  acting  to  spin  the  blade  about  its  long  axis  at  a  particular  blade 


~  cm  ^5“  SC  (4-42) 

where,  as  before,  the  dynamic  pressure  and  reference  area  are  for  the  local  conditions  at  the 
element.  The  spin  moment  coefficient  Cm ,  like  the  force  coefficients,  is  a  function  of  the  blade  tip 
speed  ratio  and  yaw  angle. 

The  total  moment  acting  to  spin  the  blade  about  its  long  axis  is  the  summation  of  the 
moments  at  all  elements  (i.e.,  the  integral  of  the  moment  distribution  over  the  blade). 


M 


CJ 


pV 2 
2 


C2dr 


(4-43) 
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Substitution  of  the  moment  expressions  (4-38),  (4-41),  and  (443)  into  the  equations  of 
motion  (4-24)  results  in  the  general  equations  of  motion  for  the  articulated  autorotor  decelerator. 
These  equations  are  shown  summarized  in  figure  10.  along  with  a  simplified  sketch  of  the 
decelerator  system  and  the  main  variables. 

The  equations  of  (4-24)  represent  the  equations  of  motion  for  one  blade.  The  analysis  to 
this  point  assumes  that  all  the  blades  act  as  independent  entities,  with  no  mutual  effect  on  each 
other.  Further,  it  is  assumed  that  all  blades  act  in  unison  in  that  the  motion  of  one  blade  is  identical 
to  all  others. 

The  first  equation  of  (4-33)  represents  the  component  of  the  aerodynamic  force  at  the 
blade  element,  which  contributes  to  the  system  drag  (i.e.,  the  force  opposite  in  direction  to  the  free 
stream  velocity  V^).  The  total  drag  force  Ds  resulting  from  one  blade  is  the  integral  of  this  force 
component  over  the  length  of  the  blade: 


DS  ~  I  f{CL~^r  C  cos  £  sina  +  Cd^y  C  cos  a  -  Cy-^C sin  £  sin  a) dr 


' y~T 1 


(444) 


The  total  system  drag  is  obtained  by  multiplying  the  drag  determined  for  a  single  blade 
by  tne  number  of  blades  included  in  the  system.  The  autorotor  decelerator  system  considered  here 
has  four  blades;  therefore: 


Total  system  drag  =  D  =  4 Dt  (445) 

E.  Additional  Considerations. 

It  should  be  noted  that  in  an  actual  physical  system,  the  blade  spin  axis  might  not 
intersect  directly  with  the  hub  spin  axis  because  the  blade  pivot  point  would  probably  be  located  on 
the  hub  at  some  radial  distance  from  the  store  center  line.  This  factor  has  not  been  included  in  this 
analysis  because  it  would  incorporate  an  additional  term  into  an  already  complicated  system.  The 
purpose  of  this  analysis  is  to  obtain  a  general  understanding  of  the  system  from  as  simplified  a  view 
as  possible.  Future  effort  could  indude  this  feature  without  altering  the  general  results  of  this 
study. 


Inspection  of  the  equations  of  (4-24)  reveals  that  the  equation  describing  the  deployment 
angle  is  a  second-order  differential  equation  and,  as  such,  could  have  an  oscillatory  form  of  solution. 
For  this  type  of  equation,  a  term  containing  0,  usually  referred  to  as  a  damping  term,  is  required. 
Even  if  the  system  described  is  basically  stable,  the  creation  of  spurious,  positive  damping  effects  in 
the  sclution  process  could  result  in  an  invalid  buildup  of  the  basic  terms  and  eventual  “blowing-up” 
of  the  solution.  This  was  seen  to  be  the  case  in  this  instance,  and  an  aerodynamic  damping  term  was 
therefore  added  to  the  moment  expression.  The  damping  term  was  assumed  to  result  from  the 
action  of  the  blade  deployment  rate  B  relative  to  the  free  stream  velocity  expressed  in  the  following 
form: 

Me  =  cmg  (fr)  Sb  sin2  6  (446) 

6  b 

where-jT  is  the  nondimensional  deployment  rate  of  the  blade  referred  to  the  free  stream  velocity. 
Cm  •  is  the  damping  moment  derivative  with  aspect  to  the  nondimensional  deployment  rate.  The 
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Figure  10.  Summary  of  the  Equations  Describing  the  Motion  of  the 
Articulated  Autorotor  Decelerator  System 

NOT  reproducible 
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Where: 


5  «  tan-1 (ctne  slna) 

tan ■'(  r  l°l  5,1,6  ) 

V. 
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c0  =  cD(  a,  o  ) 

Cy  *  Cy (  B,  0  ) 
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Figure  10.  Continued 
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term  b  is  the  distance  from  the  store  center  line  to  the  blade  midpoint,  measured  along  the  blade 
span.  The  term  sin20  introduces  the  effect  of  the  blade  angle  relative  to  the  free  stream  velocity 
(normal  component  of  b  is  b  sin  0 ). 

In  addition  to  the  moments  of  inertia  of  the  blade,  the  moment  of  inertia  of  the  hub  lH 
must  also  be  considered  in  the  £2  equation  (4-24).  Therefore,  a  factor  of  /w/4  (for  u  four-bladed 
system)  is  added  to  the  basic  moment  of  inertia  term  for  the  blade. 

V.  AUTOROTOR  SECTIONAL  AERODYNAMIC  CHARACTERISTICS. 

A.  S-Blade  Autorotor  Configuration. 

Figure  1 1  illustrates  the  type  of  autorotor  utilized  in  this  analysis.  The  shaping  of  the 
blade  cross  section  in  an  S-shape  causes  the  biade  to  autorotate.  This  cross-sectional  shape  offers 
excellent  lift  and  drag  coefficients,  and  its  geometry  allows  good  stowage  potential  (prior  to 
deployment)  along  with  good  structural  qualities  resulting  from  the  corregated  effect. 

B.  Testing  Methodology. 

In  using  the  blade  element  propeller  theory  method  outlined  in  section  IV,  the  sectional 
aerodynamic  characteristics  of  the  autorotor  blade  are  used  directly.  No  satisfactory  method  has 
been  developed  to  analytically  predict  the  characteristics  for  autorotating  shapes;  thus,  these 
characteristics  must  be  determined  experimentally  from  wind  tunnel  tests  of  the  autorotor  blade. 
Although  somewhat  different  in  testing  procedure,  because  of  the  unique  spinning  properties  of  the 
autorotor,  the  general  method  and  results  are  analogous  to  the  approach,  used  for  airfoil  sectional 
characteristics  and  their  use  in  standard  propeller  analysis. 

The  aerodynamic  characteristics  of  the  autorotor  blade  cross  section  are  determined 
experimentally  in  the  wind  tunnel  using  the  testing  arrangement  shown  in  figure  i  2.  The  rotor  is 
mounted  on  a  yoke  type  mount,  which  is  attached  to  the  tunnel  balance  system.  Each  end  of  the 
rotor  spin  axis  shaft  is  attached  to  the  upright  yoke  arms  with  high-speed  precision  bearings.  This 
allows  the  rotor  to  spin  about  its  long  axis  with  a  minimum  amount  of  friction.  The  rotor  can  be 
tested  at  various  yaw  angles  to  the  air  flow  by  rotating  the  yoke  about  its  single  center  support. 
Figure  13  shows  an  S-section  rotor  blade  model  installed  in  the  wind  tunnel.  The  yoke  is  mounted 
to  a  pyramidal  type  strain  gage  balance  located  beneath  the  test  section  floor.  This  balance  measures 
the  three  forces  and  the  three  moments  directed  along  and  about  the  three  wind  axes  directions. 
The  strain  gage  deflections  are  read  out  as  microvolts  and  are  calibrated  to  determine  the  actual 


Figure  1 1 .  S-Section  Autorotor  Configuration 
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loads  acting  on  the  autorotor.  A  strip  chart  recorder  is  used  so  that  all  six  load  components  can  be 
recorded  simultaneously  during  the  run  and  so  that  the  transient  effects  can  be  determined. 

Autorotation  spin  rates  were  obtained  by  means  of  an  optical  tachometer.  This  basically 
consists  of  a  light  source  that  projects  a  light  beam  on  the  autorotor  normal  to  the  autoroior  spin 
axis.  The  light  beam  is  positioned  so  that  the  beam  is  interrupted  by  the  surface  of  the  autorotor.  A 
light  pick  up  senses  the  reflection  of  the  beam  over  one-half  cycle.  Then  the  periodic  reflection  of 
the  light  beam  is  counted  by  the  counter  part  of  the  device,  and  the  spin  rate  of  the  autorotor  is 
recorded  on  the  strip  chart  along  with  ihe  balance  loads. 
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Figure  l'*-.  Autorotor  Blade  Installed  in  Aerodynamics  Research  Group 
14-  by  20-Inch  Subsonic  Wind  Tunnel 

The  procedure  of  only  testing  the  autoroto*-  over  a  range  of  yaw  angles  in  the  plane 
defined  by  the  rotor  spin  axis  and  the  free  stream  velocity  will  provide  the  complete  relationships 
between  the  aerodynamic  characteristics  and  the  angle  of  the  free  stream  velocity  to  the  autorotor 
for  all  combinations  of  pitch  and  yaw  angles.  This  happens  because  once  spin  is  initiated,  the 
aerodynamic  characteristics  are  not  a  function  of  the  angle  between  the  autorotor  spin  axis  and  the 
free  stream  velocity,  measured  in  the  plane  perpendicular  to  the  spin  axis.  Thus,  the  autorotor  can 
be  considered  to  be  a  purely  symmetrical  body,  aerodynamically,  thereby  permitting  this  simplified 
testing  method. 

C.  Aerodynamic  Characteristics  Nomenclature. 

The  axes  system  and  the  associated  terms  and  symbols  used  in  the  wind  tunnel  data  are 
illustrated  in  figure  1 4. 

D.  Aerodynamic  Characteristics  at  Steady  State  Spin  Rate. 

If  the  autorotor  is  set  at  a  fixed  yaw  angle  \p  in  a  constant  velocity  air  stream  F„ ,  it  will 
assume  a  constant  spin  rate  co.  The  “steady  state”  spin  characteristics  of  the  rotor  can  be  expressed 
in  nondimensional  form  as  the  steady  state  tip  speed  ratio  u>ss  where: 


A  _ 

“ss~  2  Foe 


(5-1) 


The  subscript  ss  refers  to  the  terms  evaluated  at  the  steady  state  condition.  When  discussing  the 
aerodynamics  of  the  basic  autorotor,  the  tip  referred  to  is  the  leading  edge  of  the  autorotor  shape  as 
viewed  along  the  spin  axis  of  the  autorotor.  Thus,  the  tip  speed  ratio  is  the  ratio  of  the  tangential 
velocity  of  .he  rotor  tip,  caused  by  spin,  to  the  free  stream  velocity.  The  steady  state  tip  speed  ratio 
is  a  constant  for  the  particular  rotor  cross-sectional  shape  and  the  particular  yaw  angle. 
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Figure  14.  Definitions  of  Wind  Tunnel  Axes  Systems  and  Aerodynamic  Nomenclature 


Another  parameter  that  innuences  the  steady  state  tip  speed  ratio  and  the  other 
aerodynamic  characteristics  ;s  the  rotor  aspect  ratio  (8/0-  1°  this  sense,  the  assumption 
that  the  aerodynamic  characteristics  obtained  in  the  wind  tunnel  on  a  rotor  with  an 
arbitrary  aspect  ratio  can  be  utilized  in  the  blade  element  method  for  the  analysis  of  the 
autorotor  decelerator  has  the  same  rationale  as  using  general  airfoil  section  data  in  evolving 
propeller  performance  using  the  same  method.  In  addition,  for  the  case  of  the  autorotor, 
the  combination  of  yaw  angle  with  the  consequent  span  wise  flew  would  seem  to  be 
aflected  b>  amect  ratio.  However,  other  than  using  autorotors  of  relatively  high  aspect 
ratios  in  the  wind  tunnel  tests  to  determine  the  sectional  characteristics,  no  aspect  ratio 
corrections  have  been  utilized  in  this  analysis. 

It  might  also  be  mentioned  that  the  aerodynamic  characteristics  of  autorotors  are 
also  a  function  of  Mach  number.  However,  this  analysis  does  not  involve  situations  outside 
of  subsonic  velocities,  and'  compressibility  effects  are  not  considered. 

At  a  steady  state  spin  condition,  the  autorotor  has  an  aerodynamic  lift,  drag,  and 
side  force  that  can  be  expressed  in  standard  coefficient  form  (QJS,  Cds Cys,  )  where  the 
subscrip;  ss  denotes  terms  evaluated  at  a  steady  state  spin  condition.  The  autorotor  also 
experiences  a  yawing  moment  and  rolling  moment.  But  in  the  case  of  the  S-blade  with  no 
end  plates,  considered  here,  the  yawing  and  rolling  moments  are  extremely  small  and  have 
no  meaning  because  they  are  a  result  of  side  flow  and  end  effects  for  the  blade  aspect 
ratio  and  constant  test  conditions  of  the  tunnel  that  cannot  be  applied  to  the  variable 
conditions  along  the  blade  for  the  autorotor  deceleration  system.  Thus,  they  are  not 
utilized  in  this  analysis. 

At  steady  state  spin  rate,  the  moment  about  the  spin  axis,  or  spin  moment,  is 
zero  (otherwise,  the  rotor  would  experience  spin  acceleration,  and  the  steady  state 
condition  would  not  exist).  During  the  condition  when  the  autorotor  is  spinning  up  to  the 
steady  state  spin  rate,  the  spin  moment  does  have  a  value,  and  this  will  be  treated  in 
detail  in  the  section  on  spinup  transient. 

The  steady  state  values  for  the  force  coefficients  and  tip  speed  ratio  are  shown  as 
functions  of  yaw  angle  in  figure  15.  These  parts  of  this  figure  contain  actual  wind  tunnel 
data  for  an  S-section  autorotor.  As  can  be  seen,  the  steady  state  characteristics  can  be 
expressed  as  a  function  of  yaw  angle  as: 


A  A  9  . 

“ss  -“W()C0S2^ 

(5-2) 

Ci  =  Ci  cos2  ^ 

,  ss  c  ss o 

(5-3) 

* 

c*»  =c0ssocos2* 

(5-4) 

Cy  =  Cy  sin2  2^ 

SS  SSq 

(5-5) 

where  the  subscript  (0)  refers  to  the  terms  maximum  value  (i.e.,  CV  „  and  CD  are 
evaluated  at  \p=0  ,  and  Cy  is  evaluated  at  \Ij  -  45°).  Because  the  aerodynamic 

characteristics  of  the  autorotor  can  be  expressed  in  standard  coefficient  form,  wind  tunnel 
tests  at  one  velocity  are  sufficient  to  determine  these  coefficients.  Actually,  all  wind  tunnel 
tests  were  conducted  at  several  free  stream  velocities  to  evaluate  such  effects  as  Reynolds’ 
number,  bearing  friction,  etc.  In  the  case  of  the  general  autorotor  shape,  inertial  properties, 
Reynolds’  numbers,  etc.,  used  in  this  analysis,  the  constant  coefficient  properties  assumed  in 
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A  =  6.67  in. 
C  =  2  in. 


Note:  symbol  O indicates  wind  tunnel  data 
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Figu-e  1 5.  Aerodynamic  Characteristics  as  a  Function  of  Angle  of  Attack 
for  an  S-Section  Autorotor  at  Steady  State  Spin 


the  analysis  were  verified  by  wind  tunnel  tests.  Thus,  the  wind  tunnel  data  are  presented  for  a 
single  set  of  wind  tunnel  conditions  merely  to  illustrate  a  typical  case  or  to  evaluate  a  general 
expression. 

E.  Aerodynamic  Characteristics  During  Spinup  (Transient). 

If  the  autorotor  is  set  as  a  fixed  angle  of  yaw  in  a  constant  free  stream  velocity 
and  released,  it  will  spinup  to  the  particular  steady  state  spin  rate  corresponding  to  the 
steady  state  tip  speed  ratio  as  described  in  the  previous  section.  Figure  16  indicates  how 
the  various  terms  vary  during  a  typical  spinup  in  the  wind  tunnel.  As  can  be  seen,  the 
drag  coefficient  and  side  force  coefficient  remain  fairly  constant  during  spinup;  however, 
the  lift  coefficient  and  spin  moment  coefficient  (indicated  by  the  relative  slope  of  the  tip 
speed  ratio  versus  time  curve)  vary  with  the  spin  rate  in  the  transient  condition. 

If  the  ratio  of  the  lift  coefficient  to  its  steady  state  value  is  plotted  versus  the 
ratio  of  the  tip  speed  ratio  to  its  steady  state  value,  as  shown  in  figure  17,  it  can  be  seen 
that  the  generalized  lift  coefficient  can  be  expressed  as  a  second-order  function  of  the 
generalized  tip  speed  ratio  as  follows: 


where  Q  j  and  Q,  are  constants  for  the  general  autorotor  configuration,  and  C7  __  and 
u>ss  are  steady  values  for  the  particular  autorotor  shape  and  for  a  specific  yaw  angle.  The  data  curve 
(solid  line)  in  figure  17  was  formed  using  a  second-order  curve  fit  of  the  wind  tunnel  data. 


Substituting  expressions  (5-2)  and  (5-3)  into  (5-6)  results  in: 


(5-7) 


Therefore,  the  lift  coefficient  during  the  transient  spinup  condition  can  be  expressed  as  a 
function  of  the  tip  speed  ratio,  yaw  angle,  and  the  maximum  steady  state  values  of  the  lift 
coefficient  and  tip  speed  ratio.  The  specific  lift  characteristics  of  the  S-section  autorotor,  utilized 
for  the  blade  in  the  autorotor  decelerator  system  and  presented  in  th<®  form  of  expression  (5-7),  are 
shown  in  figure  1 8.  These  data,  which  were  obtained  in  wind  tunnel  tests,  show  the  maximum 
steady  state  tip  speed  ratio  uts.  to  be  1.0,  and  the  maximum  steady  state  lift  coefficient  Cr  to 
be  1.55.  0  “° 

Equation  (5-7),  which  describes  the  transient  behavior  of  the  lift  coefficient  for  any  yaw 
angle  and  tip  speed  ratio,  is  based  on  the  assumption  that  the  transient  coefficients  Qj  and 
Ci  2  -which  were  determined  from  spinup  tests  at  a  0°  yaw  angle-are  valid  for  all  yaw  angles.  This 
assumption  was  verified  by  wind  tunnel  tests  in  which  the  autorotor  was  allowed  to  spinup  while 
being  held  at  various  fixed  yaw  angles.  Figure  18,  in  addition  to  containing  the  lift  coefficient  values 
versus  tip  speed  ratio  for  several  yaw  angles  using  expression  (5-7),  also  shows  data  points 
determined  directly  from  a  spinup  test  for  the  autorotor  at  a  yaw  angle  of  30°  as  a  means  of 
illustrating  the  justification  for  expression  (5-7). 

All  the  aerodynamic  forces  and  moments  are  measured  directly  during  the  wind  tunnel 
test  using  the  standard  tunnel  balance  system,  except  for  the  spin  moment.  The  spin  moment  is 
determined  indirectly  from  spin  rate  versus  time  traces  obtained  during  the  spinup  process.  The 
slope  of  this  curve  at  any  instant  of  time  is  the  spin  acceleration  at  that  time.  The  product  of  the 
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spin  acceleration  and  the  autorotor  moment  of  inertia  about  the  spin  axis  at  a  particular  instant  of 
time  is  equal  to  the  spin  moment  acting  on  the  autorotor  at  that  instant.  The  spin  moment  can  then 
be  reduced  to  a  spin  moment  coefficient  in  the  standard  way: 

r  -  Ik 
m  pV&SC 

where 

/  =  moment  of  inertia  about  spin  axis 
to  =  ~  =  spin  acceleration 

Thus,  the  plot  of  spin  moment  coefficient  as  a  function  of  tip  speed  ratio  can  be 
determined  by  following  this  procedure  for  a  large  number  of  points  from  the  spin  rate  versus  time 
data. 

The  most  common  method,  however,  is  to  use  a  three-point  curve  fit  from  the  spin  rate 
versus  time  data  to  form  a  second-order  equation  of  spin  moment  coefficient  as  a  function  of  tip 
speed  ratio.  The  three  points  chosen  to  evaluate  the  coefficients  are  the  initial  slope,  the  final  steady 
state  condition,  and  some  intermediate  point,  usually  the  point  where  the  spin  rate  is  one-half  its 
steady  state  value.  The  resulting  expression  is  extremely  accurate.5 

The  spin  moment  coefficient  can  be  plotted  versus  the  tip  speed  ratio  (similar  to  the 
method  used  with  the  lift  coefficient),  as  shown  in  figure  19,  and  can  be  expressed  in  the  following 
form: 


ZZL.  =  i  +  r 
Cm~  l 
m  o 


+  C, 


If) 

?“) 


m,l  / 

V 


(5-8) 


where  Cmt  and  Cmj  are  constants  for  the  general  autorotor  configuration  and  coss  is  the  steady 
state  value  of  the  tip  speed  ratio  for  the  particular  autorotor  shape  and  for  a  particular  yaw  angle. 
The  data  curve  in  figure  19  was  determined  from  wind  tunnel  data  using  the  second-order  curve  fit 
technique. 

Cm ^  is  the  value  of  the  spin  moment  coefficient  at  zero  spin  rate  (the  value  at  steady 
state  is  zero)  and  is  a  constant  for  the  particular  autorotor  shape  and  angle  of  yaw. 


=  Cmo  cos2* 

where  CmQ  is  the  value  of  the  spin  moment  coefficient  (evaluated  at  *  =  0°). 
Substituting  (5-2)  and  (5-8)  into  (5-7)  results  in: 


(5-9) 


Cm  =Cm0  COS2*  +Cmi-7 r 


& 


+  C„ 


U3* 


«o  m2<*>«2ocos2  * 


(5-10) 


5Brunk,  J.  E.  Final  Summary  Report.  Contract  DA-18-108  AMC  236(A).  Analytical  and  Aerodynamic  Studies  cf 
Rotating  Bomblet  Motion.  31  January  1965.  UNCLASSIFIED  Report. 
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Therefore,  the  spin  moment  coefficient  during  transient  spinup  can  be  expressed  as  a 
function  of  the  tip  speed  ratio,  yaw  angle,  steady  state  tip  speed  ratio,  and  the  maximum  spin 
moment  coefficient  evaluated  at  zero  spin. 

The  spin  moment  characteristics  of  the  S-section  autorotor,  illustrated  in  figure  1 1 ,  are 
presented  in  the  form  of  expression  (5-10)  in  figure  20.  These  data,  which  were  obtained  in  wind 
tunnel  tests,  show  the  spin  moment  coefficient  at  zero  spin  and  zero  yaw  angle  Cm  to  be  equal  to 
0.20.  0 


As  was  the  case  with  the  lift  coefficient,  the  transient  coefficients  Cm  0,  Cmi,  and  Cm  2 
for  the  spin  moment  coefficient,  which  were  determined  from  tests  with  the  autorotor  set  at  0°  of 
yaw,  are  valid  in  equation  (5-10)  for  all  yaw  angles  as  demonstrated  by  the  wind  tunnel  data  points 
included  in  figure  1 1  for  the  example  case  of  a  30°  yaw  angle. 

For  the  autorotor  decelerator  system  mathematical  model  evolved  in  section  IV,  the 
relative  angle  between  the  velocity  and  the  autorotor  blade  spin  axis  is  0,  where: 

0  =  90°  -  0  (4-30) 

Hence: 


sin  - 1>  -  cos  /} 
cos  =  sin  0 

Thus,  in  terms  of  the  system  angle  0,  the  expressions  for  the  steady  state  aerodynamic 
coefficients  (5-2)  through  (5-5)  become: 
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=  Ci  sin2/} 

(5-11) 
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=  CD  sin2/? 
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=  CY  sin2  2/? 
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(5-13) 

become: 


The  expressions  for  the  transient  values  of  the  aerodynamic  coefficients  (5-7)  and  (5-10) 
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F.  Aerodynamic  Characteristics  for  an  Overspin  Condition. 


(5-14) 

(5-15) 


As  described  in  the  previous  section,  the  autorotor  aerodynamics  can  to  basically 
expressed  as  a  function  of  the  tip  speed  ratio.AUp  Jo  tni^  point,  only  the  spin  regime  between  zero 
and  steady  state  has  been  considered,  i.e.,  0  £  co  £  coss. 

In  using  the  blade  element  method  in  the  analysis  of  the  rotor  decelerator  system,  if  the 
tip  speed  ratio,  yaw  angle,  local  dynamic  pressure,  etc.,  are  known  at  a  particular  element  of  the 
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blade,  the  aerodynamic  coefficient  for  that  element  can  be  determined  by  the  functional 
relationships  for  the  blade  shape.  For  the  case  of  the  rotor  retardation  system,  the  blade,  being  a 
rigid  body,  has  the  same  spin  rate  at  any  point  along  its  length.  However,  the  velocity  varies  along 
the  length  of  the  blade.  The  tip  speed,  therefore,  varies  along  the  blade,  decreasing  in  value  from 
root  to  tip.  At  some  element’s  on  the  blade,  a  condition  exists  in  which  the  tip  speed  ratio  is  greater 
than  the  steady  state  value  corresponding  to  the  conditions  at  the  element;  this  condition  is  termed 
“overspin.”  Therefore,  an  understanding  of  the  autorotor  characteristics  for  an  overspin  condition 
(i.e.,  cd  >  &ss)  is  required. 

A  wind  tunnel  test  was  conducted  in  which  the  aerodynamic  loads  were  measured  on  the 
S-section  autorotor,  which  was  forced  to  overspin  by  an  eledtric  motor.  Figure  21  illustrates  the 
rotor  and  mount  used  in  the  test.  The  results  of  this  test  are  found  in  figure  22.  As  can  be  seen,  the 
lift  and  drag  remain  at  their  respective  steady  state  values,  even  when  the  rotor  spin  rate  is  40% 
above  the  steady  state  value  for  the  test  conditions. 

Therefore,  if  to  >  to,, : 


C,  =  C, 


*«’  CD  “cz>ss>andcy  -C,Yss 


The  lift  and  drag  values  were  also  measured  with  the  motor  holding  the  rotor  spin  rates  at 
values  less  than  steady  state.  The  agreement  of  these  data  with  data  obtained  from  rotor  transient 
tests  using  the  free-spinning  rotor  method,  as  indicated  in  figure  17,  demonstrates  that  the  motor 
mounting  installation  did  not  introduce  invalid  or  spurious  conditions. 

,  »  »  i 

The  overspin  condition  also  affects  the  blade  spin  moment  coefficient  in  that  in  an 
overspin  condition,  the  blade  will  experience  a  negative' spin  moment  coefficient  tending  to  decrease 
the  'spin  rate.  A  wind  tunpel  test  was  conducted  to  determine  the  nature  of  the  spin  moment  acting 
on  the  autorotor  in  an  overspin  conditu-..  Tne  S-section  autorotor  and  yoke  mount  arrangement 
shown  in  figure  21  was  used.  An  ex  -.mal  motorized  friction  wheel  was  utilized  to  give  the 


Figure  21 .  Autorotor  Oveispin  Wind  Tunnel  Test  Mounting  Arrangment 
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s  From  Wind  Tunnel  Test  to  Determine  Effect  of  Autorotor 
Overspin  on  the  Lift  and  Drag  Coefficients 


autorotor  a  higher  spin  rate  than  the  steady  state  value  for  the  particular  free  stream  velocity  at 
which  the  test  was  conducted.  When  a  sufficiently  high  spin  rate  had  been  achieved,  the  friction 
driving  wheel  was  disengaged,  thus  allowing  the  autorotor  blade  to  freely  spin  down  to  the  steady 
state  spin  rate.  The  spin  rate  '  ersus  time  data  trace  was  analyzed  in  the  same  manner  as  for  the 
underspin  case.  The  results  of  this  test  revealed  that  the  expression  for  the  spin  moment  coefficient 
(5-10)  previously  determined  for  the  underspin  condition  remained  val'd  for  the  overspin  condition, 
at  least  up  to  a  blade  tip  speed  ratio  of  1 .45,  which  was  the  highest  value  attained  during  the  test. 
Thus,  a  very  strong  negative  spin  moment  coefficient  is  present  for  the  overspin  condition,  as 
indicated  in  figure  23. 

G.  Autorotor  Spin  Initiation. 


Almost  all  autorotating  shapes  have  the  common  characteristic  that  they  will  not 
automatically  initiate  spin  for  all  attitudes  of  the  rotor  chord  to  the  direction  of  the  free  stream 
velocity  (which  is  actually  the  definition  of  the  angle  of  attack  in  the  commonly  used  sense).  The 
angle  of  attack  in  this  instance,  as  defined  in  figure  24,  is  important  only  to  the  action  of  the  spin 
initiation.  Once  the  autorotor  is  spinning,  the  angle  of  attack  loses  meaning  in  the  description  of  the 
magnitude  of  the  blade  section  aerodynamic  characteristics.  Actually,  the  aerodynamic  forces  act  in 
a  cyclic  manner  over  a  complete  rotation.  This  periodic  nature  of  the  lift,  drag,  etc.,  can  be  observed 
a+  low  spin  rates  or  in  water.6  At  the  higher  spin  rates  common  to  the  rotors  in  air,  the  cyclic 
variances  are  too  rapid  to  be  measurable,  and  only  an  average  value  is  indicated. 


Wind  tunnel  tests  of  a  nonspinning  S-section  rotor  blade  at  various  angles  of  attack 
indicate  that  there  exists  angle  of  attack  ranges  in  which  the  rotor  is  statically  stable  in  pitching 
moment  (dead  zone)  and  ranges  where  it  is  statically  unstable  (live  zones).  These  conditions  are 
apparent  in  figure  24,  which  shows  the  pitching  moment  coefficient  as  a  function  of  angle  of  attack, 
where  a  live  zone  extends  from  - 105°  to  -60°,  and  a  dead  zone  extends  from  -60°  to  +75°,  with  a 
stable  trim  point  at  -25°.  Thus,  the  S-blade  autorotor  considered  has  four  alternating  zones  (two 
live  and  two  dead)  in  a  one  complete  360°  rotation. 


This  means  that  if  the  rotor  were  oriented  at  an  angle  of  attack  within  the  dead  zone  and 
released,  it  would  assume  a  very  stable  trim  altitude  at  an  angle  of  attack  of  25°.  If  oriented  ar  an 
angle  in  the  live  zone  and  released,  it  would  initiate  spin  and  proceed  to  spin  up  to  the  corresponding 
steady  state  spin  rate.  In  actuality,  releasing  the  rotor  in  the  outer  fringes  of  the  dead  zones  still 
might  result  in  spin  initiation  because  of  overshoot  of  the  blade  past  the  trim  point  and  into  the  live 
zone.  This  results  in  more  realistic  live  and  dead  regions,  each  covering  a  90°  range  as  indicated  in 
figure  25. 


6Dupleich,  P.  NASA  TM  1201.  Rotation  in  Free  Fall  of  Rectangular  Wings  of  Elongated  Shape.  April  1949. 
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Figure  23.  Results  From  Wind  Tunnel  Test  To  Determine  Effect  of  Autoroior 
Overspin  on  the  Spin  Moment  Coefficient 
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igure  24.  Pitching  Moment  Coefficient  as  a  Function  of  Angle  of  Attack  for 
a  Nonspinning  S-Section  Autorotor 


The  inability  of  the  autorotor  shape  to  initiate  spin  at  any  altitude  is  of  particular 
importance  in  the  autorotor  decelerator  system  application.  Experience  has  shown  that  if  the  blades 
of  the  autorotor  decelerator  system  are  manually  deployed  (i.e..  thrown  out  a  near  right  angle  to 
the  body  center  line)  without  blade  spin,  the  blades  would  assume  a  “stall”  attitude  or  angle  in 
which  they  would  not  autorotate  about  their  long  axes.  Because  this  stall  angle  results  in  the  blades 
being  at  a  slight  cant  angle,  they  will  initiate  a  steady  state  hub  spin  rate  opposite  in  direction  to 
that  which  would  be  obtained  if  the  blades  were  autorotating.  Although  the  resulting  hub  spin  rate 
is  greater  for  the  stalled  condition  than  the  autorotating  condition,  the  system  drag  of  the  former  is 
less. 


This  stalled  condition  occurs  whether  the  blades  are  deployed  aft  to  forward  or  forward 
to  aft.  It  is  extremely  difficult  to  get  all  the  blades  to  autorotate;  usually  one  or  more  are  in  the 
stalled  condition  with  a  resultant  loss  in  system  performance. 


Figure  25.  “Live”  and  “Dead”  Angle  of  Attack  Zones  for  Spin  Initiation  for 

S-Section  Autorotor 
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The  aft  to  fonvard  blade  deployment  method  (from  an  initial  hub  spin  rate)  considered 
here  always  results  in  all  the  blades  achieving  the  proper  autorotation.  This  is  a  result  of  several 
factors.  Frist,  the  blades  start  out  with  an  initial  spin  rate  caused  by  the  hub  rotation;  secondly, 
their  angular  movement  during  deployment  from  a  deployment  angle  of  0°  to  an  angle  of  near  60° 
causes  portions  of  the  blades  to  see  angles  of  attack  in  the  live  zones.  Most  important,  as  will  be 
shown  in  section  X,  a  gyroscopic  moment  due  to  the  hub  spin  rate  and  the  blade  deployment  rate 
would  seem  to  be  the  main  cause  of  blade  spin  initiation  for  the  autorotor  decelerator  system. 

Figures  26  and  27  contain  the  lift  coefficient  and  drag  coefficient,  respectively,  that  were 
obtained  during  the  wind  tunnel  test  on  the  nonspinning  autorotor.  These  data  are  included  for 
completeness  and  also  to  give  an  appreciation  for  the  effects  of  the  spin  on  the  aerodynamic  coefficients 
by  comparing  these  nonspin  values  with  the  values  at  steady  state  spin  (figure  15).  It  should  be 
noted  that  the  lift  coefficient  at  steady  state  spin  is  about  twice  the  value  of  the  maximum  lift 
coefficient  for  the  nonspinning  case.  It  is  also  interesting  to  note  that  the  maximum  spin  moment 
coefficient,  as  determined  from  the  spinup  type  tests,  figure  20,  is  much  larger  than  the  maximum 
spin  moment  coefficient  measured  on  the  nonspinning  autorotor.  This  also  emphasizes  the 
significance  of  the  effects  of  the  induced  flow  field  created  by  the  autorotation. 

VI.  METHOD  OF  SOLUTION. 

A.  General  Approach. 

The  equations  describing  the  motion  of  the  autorotor  decelerator  system  (figure  10)  were 
programmed  to  be  solved  on  a  digital  computer  The  physical  characteristics  of  the  system  (blade 
weight,  blade  length,  etc.),  the  free  stream  conditions  (velocity,  air  density,  etc.),  and  the  autorotor 
aerodynamic  sectional  characteristics  (lift  coefficient,  drag  coefficient,  etc.)  are  read  into  the 
program  for  each  run.  Starting  from  a  set  of  initial  conditions  (blade  deployment  angle,  hub  spin 
rate,  etc.),  the  program  solves  the  equations  and  computes  the  system  motion  as  a  function  of  time. 

The  aerodynamic  load  distribution  acting  over  the  blade  span  and  the  resultant  net 
moments  are  determined  by  computing  the  aerodynamic  loads  acting  at  several  elements  over  the 
blade  span  and  summing  these  to  get  the  total  load.  Thus,  all  the  local  conditions  at  each  blade 
element,  as  well  as  their  distribution  over  the  blade  span,  can  be  determined,  in  addition  to  the  net 
effect,  at  each  instant  of  time.  The  differential  equations  are  solved  using  the  improved 
Euler-Cauchy  method  of  numerical  solution.7 

The  steady  state  solution  is  obtained  by  initiating  the  program  with  arbitrary  values  for 
the  initial  conditions  for  each  term.  The  program  will  then  solve  for  the  system  motion  as  a 
function  of  time  as  the  system  experiences  a  transient  motion  convening  toward  a  steady  state 
condition.  The  steady  state  condition  is  achieved  when  the  higher  derivatives  of  9,  co,  and  fi  are 
negligibly  small,  indicating  the  absence  of  nonsteady  effects.  The  closeness  of  the  initial  condition 
values  to  the  final  steady  state  values  deterrnines  the  rapidity  with  which  the  steady  state  solution  is 
achieved. 

B.  Description  of  the  Computer  Program. 

The  form,  as  well  as  the  nomenclature,  of  some  of  the  equations  evolved  in  sections  IV 
and  V  to  describe  the  motion  of  the  autorotor  decelerator  system  have  been  altered  to  facilitate 
programming  on  the  digital  computer.  Thus,  where  applicable  in  this  section,  the  computer  program 
expressions  have  been  designated  with  the  corresponding  equation  number  from  sections  IV  and  V. 


7Kreyszic,  E.  Advanced  Engineering  Mathematics.  John  Wiley  and  Sons  Inc.,  New  York,  New  York.  1962. 
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Coefficient  as  a  Function  of  Angle  of  Attack  for  a  Nonspinning  S-Section  Autorotor 


^ASi-^-  ^ ^ ’%1' i'*4-'  i'  » r>^ s  ^  r?y 


Figure  27,  Drag  Coefficient  as  a  Function  of  Angle  of  Attack  for  a  Nonspinning  S-Section  Autorotor 


The  digital  computer  program  essentially  solves  the  three  basic  differential  equations  that 
describe  the  motion  of  the  system  in  terms  of  the  angular  accelerations  of  the  three  main  system 
variables  0,  £2,  and  co. 

.  •  *  .  Mu> 

w  =  £2  cos  0  -  0.0  sm  0  +~r~ 

7l 

O  - - I - [/.  u  cos  0  -20.0  (In  -  /, )  sin  0  cos  0  -  0<o/.  sin  0  -  M n  1 

(fa  14 +  i i  c°s20  +/0  sin20) 

0  =-p-(cof2/i  sin  0  +  O2  (70  -  /t  )sin  0  cos  0  +M'0  +Mg  ]  (4-24) 

'0 

These  equations  are  coupled  and  contain  terms  that  are,  themselves,  functionsof  various 
time  dependent  system  conditions. 

The  method  by  which  the  equations  are  solved  is  to  start  with  a  known  or  assumed  value 
for  each  motion  term  corresponding  to  conditions  at  a  given  instant  of  tinie(/0).  Although  the 
terms  actually  vary  as  a  function  of  time,  it  is  assumed  that  they  remain  constant  over  a  small  time 
increment,  (At),  providing  the  time  increment  is  extremely  small.  Using  this  assumption,  the  values 
of  the  main  system  variables  and  their  time  derivatives  can  be  determined  at  the  end  of  the  selected 
time  increment  (t  =  t0+  At)  by  using  the  improved  Euler-Cauchy  numerical  method  as  follows: 

6  =  b0+00At 

AO  -  (00  +b)nf 
0  — 0q  + AO 

O  =  O0+O  At 
co  =  OJq  +  <jj  At 

where  the  subscript  0  indicates  the  term  value  at  time  t0 . 

Using  the  values  for  the  main  system  variables  and  their  time  derivatives  computed  for  the 
new  time  (/),  the  other  system  terms  can  be  computed,  including  the  values  of  u>,  &,  and  0.  thus 
totally  describing  the  system  performance  at  time  t.  These  term  values  then  are  used  as  the  initial 
conditions  for  computing  the  values  at  the  end  of  the  next  time  increment.  This  procedure  is 
repeated,  thus  providing  a  continual  solution  as  a  function  of  time.  The  small  time  element  or 
integradon  time  used  in  this  analysis  was  At  =  0.001  second.  Check  transient  runs  made  with 
integration  time  elements  of  0.0001  seconds  did  not  indicate  any  appreciable  change  in  the  results. 

C.  Aerodynamic  Load  Distribution  Over  Blade  Span. 

The  aerodynamic  load  distribution  over  the  blade  span  is  determined  in  the  following 
manner.  The  blade  span  or  length  is  divided  into  10  equal  spar. wise  area  elements  (Ar).  The  center 
of  each  element  is  termed  a  blade  spanwise  station  and  is  identified  by  the  distance  (r)  it  is  located 
from  the  hub  spin  axis  (i.e.,  store  center  line)  as  measured  along  the  blade,  illustrated  in  figure  28. 

The  aerodynamic  conditions  and  the  associated  aerodynamic  coefficients  arc  evaluated  at 
each  blade  spanwise  station  as  follows: 


V  +  (fir  sin  0  )7- 


(4-25) 


55 


57 


The  remaining  system  performance  terms  are  computed  using  the  following  expressions: 


D  -  irJs 

n  _  £l(b  +  2/2)  sin  9 
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Figure  29  contains  a  summary  of  these  terms  and  illustrates  their  directional  sense  relative 
to  the  physical  system. 

Appendix  A  contains  a  listing  of  the  actual  computei  program  along  with  a  listing  of  the 
definitions  of  the  computer  program  terms.  Also  in  append >x  A  is  an  example  of  a  typical  solution 
printout  format.  The  first  block  of  data  for  each  run  a  ntains  the  physical  and  aerodynamic 
characteristics  of  the  various  elements  of  the  autorotor  decelerator  system,  the  free  stream 
conditions,  and  the  initial  conditions  for  the  particular  run  or  case  to  be  solved. 

The  solution  giving  the  various  term  values  describing  the  system  motion,  and  other 
conditions,  are  printed  out  at  preselected  time  intervals.  At  each  time,  two  blocks  of  data  are 
printed  out.  The  first  data  block  describes  the  aerodynamic  effects  acting  over  the  blade  span.  The 
values  of  the  various  terms  at  each  blade  element  station,  as  well  as  the  net  effect  of  certain  of  these 
terms,  are  presented.  The  second  block  of  data  printed  out  at  each  time  interval  contains  the  various 
system  motion  and  force  terms. 

It  should  be  pointed  out  that  the  computer  program  includes  terms  for  Mach  number 
effects  and  trajectory  determination;  however,  these  areas  are  not  included  in  this  analysis. 

VII.  EXPERIMENTAL  RESULTS. 

A.  General. 

A  large-scale  autorotor  decelerator  system  was  tested  in  the  University  of  Maryland  7-  by 
1 0-foot  subsonic  wind  tunnel.  The  unit  tested  actually  represented  realistic  full-scale  dimensions 
and  structuring  for  application  to  a  typical  air  delivered  store  that  was  used  as  the  forebody  in  these 
tests. 


The  hardware  was  fabricated,  and  the  test  was  conducted  under  an  Edgewood  Arsenal 
project.8  Although  this  test  was  not  conducted  specifically  for  the  thesis  purposes,  it  afforded  an 
excellent  opportunity  to  obtain  performance  data  on  a  representative  full-scale  system,  which  could 
then  be  compared  to  the  theoretical  predictions  based  on  the  thesis  material. 

B.  Model  Configurational  Details. 

Figure  30  contains  the  structural  details  of  the  unit  tested.  The  autorotor  blades  were  of 
S-cross-sectional  type  and  were  1 1  inches  long  by  3  inch  chord,  as  shown  in  figure  1 1.  The  blades 
were  made  of  aluminum  of  0.0625-inch  thickness. 

The  maximum  forebody  diameter  (i.e.,  hub  diameter)  was  5  inches.  The  blades  were  sized 
to  fit  within  this  5-inch  diameter  envelope  when  in  the  stowed  (undeployed)  condition  (i.e.,  blades 
lying  back  parallel  to  the  store  center  line). 

C.  Test  Setup  and  Procedure. 

Figure  31  shows  the  unit  mounted  in  the  wind  tunnel  test  section.  A  single  strut  mount 
was  used  to  attach  the  model  to  the  beam  balance  system  beneath  the  test  section  floor.  The  unit 
was  tested  in  0°  angle  of  attack  and  sideslip.  Steady  state  drag  load  measurements  were  made  by 
means  of  the  tunnel  balance  system.  Transient  and  steady  state  spin  rates  of  the  blades  and  hub 
were  obtained  from  a  high  speed  camera.  Blade  deployment  angles  were  measured  from  a  low  speed 
camera  viewing  the  unit  hub  from  the  side. 


8 Autorotor  Decelerator  Retardation  System.  All  American  Engineering  Company,  Wilmington,  Delaware.  Contract 
DA  18-035  AMC-356(A).  1969. 
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Figure  31 .  Full-Scale  Autorotor  Decelerator  Mounted  in  University  of 
Maryland  7- by  10-foot  Subsonic  Wind  Tunnel 


Tests  were  conducted  at  three  velocities:  100,  141,  and  173  fl/sec.  The  procedure  for 
each  test  was  as  follows.  The  autorotor  blades  were  initially  pivo.ed  back  to  lie  parallel  to  the  body 
center  line.  The  blades  were  held  in  place  by  a  single  strand  of  break  cord  looped  around  the  folded 
blades  at  their  tips.  The  break  cord  was  kept  from  slipping  off  the  blades  by  a  small  notch  in  each 
blade  edge  located  about  an  inch  from  the  blade  tip.  The  break  cord  was  of  1 1 -pound  break 
strength  and  would  break  because  of  centrifugal  force  tending  to  deploy  the  blades  at  a  hub  spin 
rate  of  '<bout  900  rpm.  This  initial  hub  spin  rate  was  obtained  by  four  small  turbine  vanes  located 
on  the  hub. 

The  turbine  vanes  initially  were  folded  down  inside  the  hub  structure  and  held  there  by  a 
cord  loop  with  a  simple  pull  release  knot,  which  cc  ild  be  released  by  pulling  it  outside  the  tunnel 
test  section. 

After  the  tunnel  had  been  brought  up  to  the  desired  velocity,  the  turbine  vanes  were 
deployed,  thus  initiating  hub  rotation.  The  hub  spin  rate  would  inc.^asc  until  the  break  cord 
holding  down  the  autorotor  blades  would  break.  The  blades  would  then  deploy,  and  the  system 
would  undergo  a  deployment  transient  motion  and  finally  achieve  a  steady  state  condition.  Figure 
32  shows  details  of  the  blades  and  turbine  vanes. 

D.  Test  Results. 

Figures  33  through  35  contain  the  steady  state  performance  results  obtained  from  these 
tests.  Two  separate  units  of  identical  configurations  were  tested  as  indicated.  The  spin  rate  of  the 
blade  (w)  and  the  hub  (ft)  are  linear  with  velocity,  as  shown  in  figure  34.  The  deployment  angle 
(0),  the  ratio  of  blade  spin  rate  to  hub  spin  rate  (oj/ft).  drag  coefficient  based  on  swept  area  (CD  ), 
and  consequently  the  effective  drag  area,  seem  to  be  constant  with  velocity.  J 
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Figure  32.  Blade  and  Hub  Details  of  Full-Scale  Autorotor  Decelerator 

A  summary  of  these  characteristic  term  values  for  the  configuration  tested  as  determined 
from  the  wind  tunnel  test  is  as  follows: 


Term 

Steady  state  v. 

e 

76.5° 

oj/q 

11.3 

q 

1.01 

Cdt 

0.58 

D/q 

2.6  sq  ft 

in  addition  to  the  quantitative  results,  which  could  then  be  compared  with  the  theoretical 
results,  the  test  indicated  other  important  factors  concerning  the  design  of  this  type  of  decelerator 
system.  Comparison  between  t!."*  results  obtained  in  the  two  separate  but  identical  units  indicates 
that  consistent  performance  can  be  achieved  with  a  large  number  of  these  items  in  spite  of  their 
seemingly  complex  nature. 

The  linear  nature  of  the  blade  and  hub  spin  rates  with  velocity  indicates  that  no  frictional 
or  unbalance  effects  are  present.  Thus,  initial  fears  that  precision  or  advanced  state-of-the-art 
bearings  would  be  necessary  are  not  substantiated.  With  regard  to  the  bearing  problem,  it  is  realized 
that  for  actual  application,  the  relatively  short  times  of  flight  mean  that  the  operational  time  is 
measured  in  seconds.  Each  run  in  the  tunnel  lasted  an  average  60  seconds.  Yet  at  the  end  of  the  test, 
the  bearings  showed  no  indication  of  wear  or  degradation.  Analysis  of  the  test  film  record  showed 
that  all  four  blades  had  the  same  spin  rates  at  any  given  time. 

It  should  also  be  mentioned  that  the  original  concept  of  the  deployment  angle  degree  of 
freedom  was  to  relieve  the  large  bending  moments  on  the  blade/hub  attachment  point.  Early 


64 


Rate  (a>)  ^  rpm 


Figure  35.  Steady  State  Drag  Coefficient  and  Effective  Drag  Area  as  Determined  From 
Wind  Tunnel  Tests  on  a  Full-Scale  Autorotor  Decelerator 


autorotor  decelerator  designs  all  had  this  “flapping”  freedom,  but  rotation  of  the  blade  in  a  plane 
normal  to  the  hub  spin  axis  referred  to  as  “lead-lag”  was  restricted.  Tests  with  these  early  units 
resulted  in  the  blade  shafts  fracturing  completely.  Analysis  of  the  failure  indicated  that  the  shafts 
were  experiencing  large  bending  loads  because  of  the  restrictions  of  the  movement  in  the  lead-lag 
direction.  Any  unbalance  or  extraneous  loads  of  the  blade  results  in  small  precession  about  the 
blade  spin  axis.  If  this  ihotion  is  restricted,  the  resultant  moment  is  taken  out  by  the  blade  shaft.  In 
addition,  these  .loads  occurred  at  a  high  frequency  (i.e.,  that  of  the  blade  spin  rate),  resulting  in  a 
fatigue  effect  on  the  blade  shaft  and  consequent  structural  failure.  The  units  used  in  the  test 
described  here  were  designed  to  give  the  blade  a  few  degrees  of  angular  motion  freedom'  in  the 
lead-lag  direction.  This  relieved  the  high  loads  resulting  from  the. small  precession,  but  allowed  the, 
basic  hub  torque  moment  to  be  realized. 

VIII.  CORRELATION  BETWEEN  EXPERIMENTAL  AND  ANALYTICAL  RESULTS. 

■  The  physical  properties  of  the  rotor  system  and  the  test  conditions  for  the  wind  tunnel 
test  were  used  in  the  computer  program  described  in  section  VI,  and  the  analytical  steady  state 
solution  was  obtained  based  on  the  equations  derived  in  section  V.  , 

i  •  < 

‘  '  /  ' 

The  conditions  simulated  and  the  resultant  output  are  included  in  appendix  B.  A 
comparison  between  the  cpmputed  theoretical  results  and  the  experimental  results  are  shown 
below  for  a  free  stream  velocity  of  100  ft/sec.  As  can  be  seen,  the  agreement  is  not  satisfactory. 


Term 

Experimental 

results 

Analytical 

results 

n 

850  rpm 

1 

1 330  rpm 

03 

9650  rpm 

1 1200  rpm 

c o./Sl 

11.3 

8.4 

e 

76.5° 

80.8° 

cdt 

0.58 

0.84 

DM  ' 

•  2.6  sq  ft 

3.8  sq  ft 

SI 

1.01 

1.64 

A  general  parameter  study  was  instituted  to  determine  if  any  term  was  so  sensitive  that  a 
very  accurate  measurement  was  necessary  and  also  to  indicate  any  term  or  terms  that  could  be 
altered  to  give  better  results.  The  procedure  was  to  solve  for  the  steady  state  results  for  conditions 
in  which  a  single  parameter  was  varied  on  . incremental  value  above  and  below  the.  nominal  value. 
These  parameters  included:  blade  weight,  blade  length,  blade  chord,  free  stream  velocity,  blade 
aerodynamic  coefficients,  blade  aerodynamic  tip  losses,  and  bearing  friction. , 

r 

The  procedure  in  solving  for  the  steady  state  conditions  was  to  start  with  assumed  initial 
conditions  for  each  term.  The  motion  of  the  system  would  then  be  solved  for  a  function  of  time;  by 
1  or  2  seconds,  the  values  were  essentially  constant,  this  being  the  steady  state  condition. 

’  ■  I 

A  summary  of  the  results  of  this  .study  are  shown  in  table  I.  This  table  presents  the 
various  terms  of  interest  along  the  left  side  with  the  trend  (increasing  value  or  decreasing  value) 
required  to  bring  the  analytical  value  more  in  line  with  the  experimental  value.  The  most  important 
term  is  the  hub  spin  rate  in  that  it  has  the  greatest  effect  on  the  other  term  values.  As  can  be  seen, 
altering  of  none  of  the  parameters  results  in  the  desired  trends  for  all  terms  considered.  Even  tip 
losses  were  considered.  This  can  be  done  by  assuming  a  linear  lift  distribution  from  the  value  at  a 
selected  blade  element  to  a  value  of  zero  at  the  tip.  Although  the  trends  of  the  terms  are  all  in  the 
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Table  I.  Summary  of  Results  of  Initial  Parameter  Study  to  Determine  Reasons 
for  Discrepancies  Between  Theoretical  and  Experimental  Data 


1 

Term 

Parameter 

Required* 

trend 

ft** 

CO 

CO 

SI 

0 

cdt 
_ _ 

D 

<7 

A 

ft 

Blade  weight  (vv^;) 

t 

\ 

- 

t 

t 

1 

1 

t 

Blade  length  (C) 

t 

i 

- 

t 

1 

l 

t 

i 

Blade  chord  ( C) 

t 

i 

i 

i 

t 

t 

t 

i 

t 

Blade  tip  speed 
'  ^tio  <$)SS() 

i 

i 

i 

i 

t 

1 

1 

I 

Tip  losses 

t 

■  i 

i 

t 

1 

l 

i 

i 

Blade  bearing 
,  friction 

t 

i 

i 

i 

I 

l 

i 

Desired  trend  to 
match  analytical^ 

■  and  experimental 
results 

i 

t 

l 

ir 

i 

‘This  column  indicates  whether  the  parameter  value  had  to  be  increased  (t)  or  decreased  (I)  relative  to  the  nominal 
parameter  value  for  the  analytically  determined  value  of  the  huh  spin  rate  to  equal  the  experimental  value. 

“Columns  3  through  7  indicate  whether  the  analytically  determined  value  of  the  indicated  term  is  greater  ft)  or  smaller 
(f)  than  the  experimental  value  when  the  analytical  hub  spin  rate  equals  the  experimental  value. 


desired  direction,  matching  the  analytical  with  the  experimental  values  for  all  terms  could  not  be 
done.  Frictional  effects  were  simulated  by  limiting  the  blade  spin  rate  to  various  fractions  of  the 
steady  state  value  of  the  no  friction  case. 

With  the  inability  of  the  initial' parameter  study  to  provide  an  answer  to  the  discrepancies 
between  the  theoretical  and  experimental  results,  attention  was  centered  on  the  aerodynamic  load 
distribution  over  the  blade.  It  should  be  noted  that  the  basic  method  utilized  in  analyzing  the 
performance  of  the  autorotor  system  was  to  determine  the  motion  and  forces  resulting  from  a  single 
blade  and  then  to  multiply  this  by  the  number  of  blades  where  required  to  obtain  system  effects.  It 
was  assumed  that  the  blade  is  not  influenced  by  the  presence  of  the  other  blades. 

■Figure  36  shows  in  graphic  form  the  conditions  at  several  stations  along  the  blade  at  the 
condition  of  steady  state  spin.  It  illustrates  how  the  magnitude  and  direction  of  the  relative  velocity 
changes  from  the  elements  near  the  blade  root  to  elements  near  the  blade  tip.  It  also  shows  ■  *  me 
lift  and  drag  forces  combine  to  give  a  resultant  force  that  acts  to  spin  the  blade  about  the  ’mb  -.»d 
also  to  provide  the  system  drag.  It  should  be  noted  that  elements  near  the  root  result  in  c  jm,'j  spin 
force  tending  to  spin  the  hub  faster,  whereas  elements  near  the  tip  tend  to  resist  hub  spin  rate.  This 
illustration  shows  how  the  lift  force  contributes  to  increasing  the  hub  spin  rate,  whereas  the  drag 
force  acts  to  decrease  the  hub. spin  rate.  Also  clearly  shown  is  how  both  the  lift  and  drag  forces 
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View  normal  to  Conditions  at  stations  indicated 

blade  spin  axis  viewed  along  blade  spin  axis,  from 

tip  toward  root 

Figure  36.  Aerodynamic  Conditions  at  Various  Blade  Spanwise  Stations  at  Steady  State 
Condition  for  Analytical  Approach  With  Interference  Factor  (AT  =  0) 
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contribute  to  the  system  drag,  and  the  advantageous  effect,  of  the  high  lift  and  high  drag  inherent  in 
the  autorotor  blade,  toward  producing  the  high  drag  of  the  autorotor  dccelerator  system. 

Inspection  of  the  aerodynamic  load  distribution  over  the  blade  reveals  that  a  rotation  of 
the  velocity  vector  at  each  element  (i.e.,  an  incremental  change  in  the  relative  angle  of  attack  a) 
would  cause  a  rotation  of  the  lift  and  drag  directional  sense  at  each  element.  The  general  result 
would  be  to  increase  the  magnitude  of  the  force  acting  to  spin  down  the  hub,  as  shown  in  figure  37. 
It  should  also  be  pointed  out  that  a  decrease  in  hub  spin  rate  would  also  decrease  the  magnitude  of 
the  relative  velocities  at  each  element,  thus  decreasing  the  aerodynamic  forces,  and  influencing  all 
the  other  terms  in  the  system,  and  in  particular,  reduce  the  system  drag.  All  of  the  resulting  effects 
of  changing  the  angle  of  attack  at  each  element  by  an  incremental  amount  Aa  would  seem  to  bring 
the  analytical  and  the  experimental  results  more  in  agreement. 


Runs  were  then  made  in  which  the  angle  of  attack  at  each  blade  element  was  altered  a 
fixed  amount  Aa 


a  =  tan" 1 


I r  1  Q.  I  sin  0 ' 


+  Aa 


(8-1) 


The  Aa  value  was  a  constant  at  all  elements.  Runs  were  made  with  various  values  for  Aa. 


It  should  be  pointed  out  that  the  viewing  sense  along  the  blade  spin  axis,  indicated  in 
figure  36,  does  not  present  the  lift  and  drag  forces  truly  relative  to  the  hub  torque  and  system  drag 
directions  due  to  the  angle  £,  at  which  they  are  inclined.  At  the  steady  state  condition  being 
considered,  however,  when  the  blade  deployment  angle  is  large  and  the  relative  velocity  at  all  points 
on  the  blade  span  is  essentially  normal  to  the  blade  spin  axis,  the  effect  of  the  angle  £  can  be 
neglected  for  the  purpose  of  this  discussion.  At  the  steady  state  condition,  the  effect  of  the 
aerodynamic  side  force  is  negligibly  small;  thus,  only  the  lift  and  drag  forces  are  considered. 

Excellent  agreement  was  obtained  fora  Aa  of  about  15°.  A  possible  physical  explanation 
for  this  interference  effect  is  that  it  results  from  the  flow  field  about  a  blade  being  influenced  by 
the  flow  fields  of  the  other  blades  in  the  system. 

Figure  38  shows  a  typical  flow  field  about  a  single  autorotating  shape9  and  illustrates  the 
strong  positive  and  negative  vortices  each  shed  periodically  during  one  revolution.  A  relatively  large 
wake  is  produced  with  a  pronounced  downwash  angle.  These  flow  field  features  are  typical  of  all 
autorotating  shapes.  It  is  easy  to  see  how  the  flow  Fields  of  the  preceding  and  following  blades  could 
have  a  turning  effect  on  the  flow  over  a  blade  in  the  autorotor  decelerator  system.  The  physical 
reason  for  this  angular  change  over  the  original  approach  was  probably  an  interference  effect  or 
induced  effect  caused  by  the  blade  flow  fields.  It  then  was  decided  to  express  it  as  a  function  of  the 
lift  coefficient.  This  follows  the  usual  downwash  or  interference  practice  and  also  makes  it,  in 
effect,  proportional  to  the  spin  rate  of  the  blade,  which  would  seem  reasonable.  In  addition,  it 
allows  the  effect  to  be  physically  valid  during  the  deployment  sequence  in  that  this  induced  angle 
would  be  zero  when  the  blade  did  not  have  a  spin  rate.  The  new  expression  for  angle  of  attack  is 
therefore: 


Aa  =  KC,  (8-2) 

where  K  is  an  arbitrary  constant  of  proportionality. 


9Flow  Analysis  and  Dynamics  of  Autorotating  Lifting  Bodies.  University  of  Notre  Dame,  Notre  Dame.  Indiana. 
Contract  DA  18-035  AMC-356(A).  1964. 
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Figure  38.  Typical  Flow  Field  Around  a  Single  Autorotor 


Figures  39  through  41  show  the  steady  state  results  computed  for  various  values  of  K  for 
the  three  velocities  used  during  the  wind  tunnel  test.  Also  indicated  are  the  data  values  obtained 
during  the  test  to  indicate  the  value  of  K  that  best  matches  the  experimental  data.  A  value  of 
K  =  1 0  results  in  the  values  shown  below  as  compared  to  the  experimentally  determined  values: 


Term 

Experimental 

Analytical 
results 
(K=  10) 

9 

850  rpm 

830  rpm 

GO 

9650  rpm 

9200  rpm 

co/9 

11.3 

1 1.1 

e 

76.5° 

76.6° 

cdt 

0.58 

0.576 

Die, 

2.6  sq  ft 

2.55  sq  ft 

9 

1. 01 

1.00 

Table  II  compares  the  experimental  results  with  the  analytical  results  obtained  without  an 
interference  effect  (K  =  0)  and  the  analytical  results  obtained  with  an  interference  effect  ( K  =  1 0). 
The  i’t'odynamic  conditions  over  the  blade  span  are  shown  in  figure  42  for  the  K  =  10  case,  and  can 
be  compared  with  K  =  0  case  shown  in  figure  32.  Note  in  figure  42  that  the  angle  of  attack  at  each 
station  is  slightly  larger  than  those  in  the  A'  =  0  case  (figure  36).  In  other  words,  the  new  angles  of 
attack  are  not  merely  the  old  plus  the  correction  angle  (Aa).  The  effect  of  adding  the  Aa  term  alters 
all  other  terms  in  a  feedback  type  effect,  resulting  in  all  new  conditions.  The  complete  steady  state 
conditions  for  the  K  =  10  case  are  included  in  appendix  B. 

NOT  REPRODUCIBLE 
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Note:  symbol  Vindicates  the  value  of  the  particular  term 
obtained  experimentally  from  the  full  scale  autorotor 
decelerator  wind  tunnel  test 
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Figure  39.  Steady  State  Performance  of  Full-Scale  Autorotor  Decelerator 
Computed  as  a  Function  of  Interference  Factor  (K)  for  a  Free  Stream 
Velocity  of  100  Ft/Sec 
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Figure  39.  Continued 


The  value  of  K  could  be  a  function  of  other  terms,  such  as  number  of  blades,  autorotor 
blade  shape,  solidity  factor,  lift  coefficient,  etc.  It  is  most  certainly  a  function  of  the  lift  coefficient. 
The  constant  nature  indicated  in  this  case  could  be  a  result  of  the  nea-  constant  lift  coefficient  over 
the  blade  span. 

Figures  43  through  45  indicate  the  transient  motion  during  deployment  determined  from 
the  computer  program  using  K=  10,  compared  with  the  data  obtained  during  the  wind  tunnel  test. 

It  should  be  noted  that  the  aerodynamic  damping  coefficient  was  initially  selected 
arbitrarily  to  be  the  lowest  value  required  to  keep  the  system  from  “blowing  up.”  Without  a 
sufficiently  large  value  of  Cm  .,  0  was  found  to  diverge  after  the  system  had  achieved  a  near  steady 
state  condition.  This  diverging  oscillatory  behavior  of  0  resulted  in  the  other  system  terms  diverging. 

The  existence  of  an  aerodynamic  moment,  which  is  a  function  of  0,  is  physically  realistic, 
especially  in  view  of  the  large  values  of  9  realized  during  deployment. 

Several  values  of  Cm  •  were  used  ir.  figures  43  through  45.  As  can  be  seen,  the  lower  the 
damping  value,  the  more  oscillatory  is  the  transient  motion.  A  low  value  of  damping  would  seem  to 
be  present  during  the  initial  phase  of  the  deployment,  evidenced  by  the  large  angular  motion  of  0  in 
figure  43.  The  absence  of  additional  oscillatory  motion  after  this  initial  phase  may  indicate  that  the 
damping  effect  may  be  a  function  of  factors  other  than  those  assumed  in  the  analysis  and  may  not 
be  of  a  constant  nature. 
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An  appreciation  of  the  rapidity  of  motion  involved  is  gained  by  noting  that  the  system 
achieves  steady  state  conditions  in  about  1 Vi  seconds  after  deployment  initiation. 

It  should  be  noted  from  figure  43  that  initial  deployment  motion  of  the  blade  takes 
about  0.1  second.  From  figures  44  and  45,  it  can  be  determined  from  the  average  hub  and  blade 
rotation  rates  during  this  time  interval  that  the  blade  makes  only  one  complete  revolution,  and  the 
hub  makes  about  3/4  of  a  revolution. during  this  initial  deployment  phase. 

IX.  STEADY  STATE  PERFORMANCE  PARAMETER  STUDY. 

A.  Method  of  Analysis,. 

A  study  was  conducted  to  evaluate  the  effects  of  the  various  physical  and  aerodynamic 
properties  on  the  steady  state  performance  of  the  autorotor  decelerator  system.  A  nominal  set  of 
values  for  the  rotor  physical  characteristics,  aerodynamic  characteristics,  and  air  flow  conditions 
were  selected  that  were  representative  of  a  realistic  physical  case.  The  steady  state  performance  of 
the  rotor  system  for  this  nominal  case  was  computed  using  the  equations  of  motion  programmed  on 
the  computer. 


81 


Note:  symbol  V  Indicates  the  value  of  the  particular  term 
obtained  experimentally  from  the  full  scale  autorotor 
decelerator  wind  tunnel  test 


Figure  41 .  Steady  State  Performance  of  Full-Scale  Autorotor  Decelerator 
Computed  as  a  Function  of  Interference  Factor  ( K )  for  a  Free  Stream 
Velocity  of  1 73  Ft/Sec 


K  ^  Deg. 
Figure  41.  Continued 


The  steady  state  performance  then  was  calculated  for  the  cases  in  which  a  selected 
parameter  was  given  a  value  larger  than  the  nominal  value  and  smaller  than  the  nominal  value.  Table 
III  contains  a  summary  of  the  terms  and  values  used  in  the  parameter  study.  The  high  and  low 
values  selected  for  each  term  represent  realistic  ranges.  Appendix  C  contains  the  computer  printout 
showing  the  steady  state  conditions  for  each  parameter  run. 

To  facilitate  discussion,  the  performance  of  the  rotor  system  is  plotted  versus  a  specific 
parameter  for  a  constant  value  of  K.  A  value  of  K  =  10  was  selected  because  this  value  of  K  was 
determined  for  the  full-scale  autorotor  decelerator,  which  closely  resembles  the  nominal  case  for  the 
parameter  study. 

As  stated  previously,  the  interference  effect,  and  thus  the  value  of  K,  may  be  a  function 
of  several  terms.  In  this  parameter  study,  however,  it  is  assumed  that  the  value  K  =  10  holds  for  the 
values  of  the  parameters  used.  These  data  will  be  presented  in  three  parts:  rotor  system  physical 
properties,  flow  conditions,  and  rotor  blade  aerodynamic  characteristics. 

The  terms  with  values  plotted  as  a  function  of  each  parameter  include  the  following: 
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Table  II.  Comparison  Between  Experimental  Results  and  Analytical  Results 
Computed  With  and  Without  an  Interference  Factor  ( K )  for  a  Free 
Stream  Velocity  of  100  Ft/Sec 


Term 

Experimental 

results 

Analytical 
results 
(*  =  0) 

Analytical 
results 
(A  =  10) 

SI 

850  rpm 

1330  rpm 

830  rpm 

CO 

9650  rpm 

11200  rpm 

9200  rpm 

oj/SI 

11.3 

8.4 

11.1 

e 

76.5° 

80.8° 

76.6° 

CDr 

0.58 

0.84 

0.576 

D/q 

2.6  sq  ft 

3.8  sqft 

2.55  sq  ft 

SI 

1.01 

1.64 

1.00 

a  J2(i+^)sinfl 

1 .  Tip  speed  ratio,  SI  = - r; - ,  is  the  ratio  of  the  tangential  velocity  of  the  blade 

V  oo 

tip  in  a  plane  perpendicular  to  the  hub  spin  axis  and  the  free  stream  velocity.  This  is  the  inverse  of 
the  advance  ratio  and  is  a  common  nondimensional  term  used  in  describing  the  rotational 
characteristics  of  autorotating  or  spinning  type  devices. 

2.  Deployment  angle,  6 ,  is  the  angle  that  the  blade  spin  axis  (i.e.,  long  axis)  makes  with 
the  body  center  line. 

3.  Blade  spin  rate,  w ,  is  the  spin  rate  of  the  blade  about  its  long  axis. 

4.  Hub  spin  rate,  SI ,  is  the  spin  rate  of  the  blade  long  axis  about  the  body  center  line 
(i.e.,  the  hub  spin  rate). 

5.  Blade/hub  spin  ratio,  u/Sl ,  is  the  ratio  of  the  blade  spin  rate  (w)  to  the  .rub  spin  rate 

(«). 

6.  Effective  drag  area,  D/q,  is  the  total  drag  of  the  system  expressed  in  terms  of  the 
free  stream  dynamic  pressure. 

7.  Drag  coefficient,  CD  =-5— ,  the  drag  coefficient  of  the  system  based  on  the  total 

qsT 

projected  area  swept  over  by  the  blades  where: 

ST  =  ir(b  +-|^2  sin2  0 

This  coefficient  represents  the  efficiency  of  the  system  to  produce  drag  for  a  given  projected  frontal 
area. 

Although  only  three  data  points  are  used  to  construct  the  curves  for  each  parameter,  the 
curves  would  seem  to  present  valid  trends  because  of  the  absence  of  any  severe  nonlinearities  under 
the  assumptions  used  in  the  analysis  and  the  relative  closeness  of  the  parameter  value  (i.e.,  same 
order  of  magnitude)  used  to  construct  the  curve  in  each  case. 
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Figure  42.  Aerodynamic  Conditions  at  Various  Blade  Spanwise  Stations  at  Steady 
Steady  State  Condition  for  Analytical  Approach  With  Interference 
Factor  (A-  =10) 
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Table  III.  Summary  of  Term  Values  Used  in  Articulated  Autorotor  Deceierator 
Steady  State  Performance  Parameter  Study 


VO 

tO 

z 

z 

z 

z 

z 

z 

CO 

CO 

ro 

00 

z 

z 

00 

to 

00 

o 

2 

z 

o 

o 

<N 

00 

00 

ro 

o 

to 

CO 

CO 

z 

z 

CM 

CO 

Tj- 

co 

z 

z 

o 

o 

o 

o 

O 

R 

to 

OO 

O 

r- 

to 

CM 

f-H 

oo 

00 

VO 

vO 

CM 

00 

vn 

VO 

CM 

1— « 

CM 

Ov 

CO 

r- 

O 

O 

o 

o 

o 

o 

o 

8 

O 

O 

8 

o 

o 

o 

o 

o 

d 

d 

oo  m  o 

VO  CO  CO 
00  On  f 
tr>  CN 

O  O  -h 

O  O  O 
O  ©  O 
O  O  O 

odd 


o  o 

VO 
Ov  co 

o~ZZZZZZZZZZZ 

©  O 

©  o 
©  o 


oo  to  n 

8zz£S3zzzzzzzzzzzzz 

O  O  O 


CO  ©  ©  Ov  00 

SZZ^S^^ZZZZZZZZZZZ 
o  dodo 


^ZZZZ'^'tZZZZZZZZZZZ 

Szz^^zzzzzzzzzzzzz 


O  VI  O  O  V)  Tf  r~ 

t  O'  O  t"  Os  00 

to^IScs^rs^ZZZZZZZZZZZ 

o  o  o  d  o  o  o 

§zr-vzzz£gzzzzzzzzz 


S';  o  o 

r**  cs  \o 

co  w  ~  r- 

oz  zzzzzzzSSzzzzzzz 

©  O  O 

o  o  o 


ni,fi\ONcoo\o*H^ffit^vot,'Oo 


91 


‘  v*fSp****r'  -*•<?  *  rp^wgQfi+qy  ”rWSa9C^*»m*TW«»*^mv*^W^^  ■'■* 


% 


% 

i 

c? 


The  main  purpose  of  these  curves  is  to  indicate  the  general  effect  that  varying  a  single 
parameter  term  has  on  the  decelerator  system  performance.  Of  special  interest  is  the  sensitivity  of 
the  decelerator  performance  on  a  particular  parametric  term.  Thus,  the  curves  developed  should 
provide  vaiid  and  meaningful  results. 

B.  i  Effect  of  Rotor  System  Physical  Properties. 

1 .  Effects  of  Blade  Weight. 

t 

The  effects  of  the  weight  of  the  autorotor  blade  on  the  steady  state  performance  are 
shown  in  figures  46  through  49.  It  should  be  noted  that  the  nominal  value  for  blade  weight  would 
represent  a'  blade  fabricated  from  aluminum.  The  heavier  weight  would  represent  a  steel  blade,  and 
the  lighter  weight  would  be  a  boron-fibre  type  material.  The  most  significant  effect  of  blade  weight 
is  that  the  blade  deployment  angle  increases  with  increasing  blade  weight  because  of  the  mass  effect 
on  the  centrifugal  force.  The  hub  spin  rate  increases  slightly  with  decreasing  blade  weight.  This 
makes  up  for  the  reduced  projected  frontal  area  from  the  lower  deployment  angle,  and  the  result  is 
J  that  the  total  drag  produced  is  independent  of  blade  weight.  The  drag  coefficient  increases  with 
reduced  weight  because  the  drag  results  from  a  reduced  projected  area. 

It  should  be  remembered  that  changes  in  the  system  weight  and  dimensions  are  also 
included  in  the  mpments  of  inertia,  which,  in  turn,  affect  the  inertial  and  gyroscopic  terms. 

2.  Effects  of  Blade  Length. 

The  effects  of  blade  length  on  the  performance  are  included  in  figures  50  through  53.  It 
should  be  noted  that  the  changes  in  blade  length  also  affect  corresponding  changes  in  the  mass  and 
inertial  properties  of  the.  blade.  The  deployment  angle  decreases  slightly  with  blade  length.  The 
blade  spin  rate  remains  the  same  regardless  of  blade  length.  The  hub  spin  rate  decreases  with 
increasing  blade  length;  however,  the  tip  speed  ratio  remains  practically  constant.  The  system  drag 
increases  with  increased  blade  length,  but  the  drag  coefficient  decreases,  indicating  a  loss  of 
efficiency  with  increased  blade  length.  The  blade-to-hub  spin  ratio  increases  with  increased  blade 
length,  reflecting  the  reduced  hub  spin  rate. 

3.  Effects  of  Blade  Chord. 

Figures  54  through  57  show  the  effects  of  the,  blade  chord  length  on  the  performance 
characteristics.  The  mass  and  inertial  values  of  the  blade  correspond  to  the  changes  in  the  blade 
chord.  Here,  the  hub  spin  rate  remains  constant  \Vith  chord  length,  but  the  blade  spin  rate  decreases 
markedly  with  increasing  chord  length.  The  effect  of  chord  length  on  the  blade  spin  ratio  is 
consistent  with  basic  autorotor  performance.  Interestingly,  the  deployment  angle  exhibits  a 
minimum  value  near  the  nominal  conditions,  which  is  a  result  of  the  interactions  of  the  inertial 
gyroscopic  and  aerodynamic  effects.  These  deployment  angle  differences,  however,  do  not  result  in 
■  any  effect  on  the  tip  speed  ratio,  which  remains  constant  in  accordance  with  the  hub  spin  rate.  The 
effective  drag  area  increases  linearly  with  increasing  chord,  as  does  the  drag  coefficient.  Note  that  a 
chord  of  4  inches,  which  is  not  ari  unrealistic  size  for  the  configuration  considered,  results  in  a 
relatively  high  drag  coefficient  of  0.78. 

4.  Effects  of  Solidity  Factor. 

The  effects  of  blade  chord  and  blade  length  also  can  be  expressed  as  a  function  of  the 
solidity  factor  (a),  which  is  defined  as  the  ratio  of  the  total  areas  of  the  blades  to  the  maximum 
projected  swept  qrea. 
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Figure  46.  Effect 
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Figure  53.  Effect  of  Blade  Length  on  Steady  State  Effective  Drag  Area 
and  Drag  Coefficient 
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Figure  55.  Effect  of  Blade  Chord  on  Steady  State  Deployment  Angle  and 
Blade/Hub  Spin  Ratio 
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where 


_  4%C 

a~T- 

t 

S-  X  C2£1~-S83)2  =  maximum  swept  area 
4  (blades  at  0  =  90°) 

These  data  are  contained  in  figures  58  through  61. 

In  effect,  an  increase  in  solidity  factor  gives  the  same  results  as  an  increase  in  the  blade  or 
chord  length,  as  shown  previously.  However,  it  is  interesting  to  note  the  linear  nature  of  the  drag 
coefficient  with  solidity  factor.  The  nonlinear  effects  indicated  can  be  attributed  to  the  influence  of 
the  inertial  terms  (gyioscopic  terms)  that  are  present  even  at  steady  state  conditions.  It  should  be 
emphasized,  however,  that  the  interference  factor  K  could  be  a  function  of  the  solidity  factor, 
among  other  things,  and  could  influence  the  shape  of  this  curve. 

C.  Effect  of  Air  Flow  Conditions. 

1 .  Effects  of  Free  Stream  Velocity. 

The  data  showing  the  effects  of  free  stream  velocity  are  included  in  figures  62  through 
65.  Although  the  blade  and  hub  spin  rates  increase  linearly  with  increasing  free  stream  velocity,  the 
remaining  terms  are  constant,  including  the  blade-to-hub  spin  ratio. 

2.  Effects  of  Air  Density. 

The  effects  of  air  density  are  shown  in  figures  66  through  69.  The  densities  investigated 
were  half  the  nominal  density  that  would  correspond  to  an  altitude  of  2C,)00  feet  and  double  the 
density  that  admittedly  is  not  a  realistic  value  for  air  but  may  have  meaning  for  another  fluid  if 
Reynolds’  number  effects  are  not  considered.  The  results  are  surprising  in  that  the  increased  density 
results  in  a  reduced  steady  state  deployment  angle,  probably  due  to  the  large  aerodynamic  loads 
relative  to  the  system  inertial  loads.  This  causes  a  slight  increase  in  hub  spin  rate  with  increasing 
density,  resulting  in  about  the  same  effective  drag  area  as  for  the  lower  density  values.  Since  the 
frontal  area  is  reduced  slightly,  the  drag  coefficient  is  constant  with  air  density. 

D.  Effect  of  Autoro'or  Blade  Aerodynamic  Characteristics. 

1 .  Linearized  Aerodynamic  Effect. 

As  described  in  section  V,  the  aerodynamic  characteristics  of  the  autorotor  blade  are 
described  as  a  function  of  the  steady  state  t;p  speed  ratio  evaluated  at  zero  yaw  angle  (uSSq).  This 

term  not  only  describes  the  blades  spin  moment  coefficient,  but  also  its  sectional  lift  and  drag 
characteristics.  As  can  be  seen  in  figures  18  and  20,  the  aerodynamic  coefficients  are  actually 
expressed  by  a  quadratic  function  of  the  tip  speed  ratio.  To  simplify  the  expressions  used  in 
investigating  the  effects  of  the  aerodynamic  characteristics,  the  spin  moment  coefficients  are 
expressed  as  a  linear  function  of  tip  speed  ratio,  as  shown  in  figure  70.  This  linearization  makes  it 
possible  to  understand  the  effects  of  vaiying  the  aerodynamic  coefficients  in  a  general  manner 
without  complicating  the  results  with  second  order  effects.  The  difference  in  the  steady  state 
performance  between  the  linear  and  nonlinear  nominal  case  is  small,  as  shown  in  figures  70  through 
73.  This  is  due  to  the  relative  symmetry  and  the  similar  slope  values  of  the  coefficient  curves  in  the 
vicinity  of  the  steady  state  tip  speed  ratio.  Inspection  of  the  computer  output  for  runs  1  and  12  in 
appendix  C  shows  that  the  effects  across  the  blade  span  are  similar.  The  linear  case,  with  its 
decreased  lift  coefficient  in  the  underspin  region,  results  in  a  slightly  lower  hub  spin  rate  and  lift 
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Fipire  58.  Effect  of  Solidity  Factor  on  Steady  State  Tip  Speed  Ratio 
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Figure  60.  Effect  of  Solidity  Factor  on  Steady  State  Blade  Spin  Rate 
and  Hub  Spin  Rate 


Figure  62.  Effect  of  Free  Stream  Velocity  on  Steady  State  Tip  Speed  Ratio 
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Figure  66.  Effect  of  Air  Density  on  Steady  State  Tip  Speed  Ratio 
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Figure  69.  Effect  of  Air  Density  on  Steady  State  Effective  Drag  Area  and  Drag  Coefficient 
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Figure  70.  Effect  of  Blade  Sectional  Steady  State  Tip  Speed  Ratio  on  Steady  State  Tip  Speed  Ratio 
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Figure  73.  Effect  of  Blade  Sectional  Steady  State  Tip  Speed  Ratio  on  Steady  State 
Effective  Drag  Area  and  Drag  Coefficient 
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coefficient.  The  reduced  spin  moment  coefficient  of  the  linear  case  results  in  a  slightly  reduced 
blade  spin  rate. 

2.  Effects  of  Blade  Steady  State  Tip  Speel  Ratio  (w„Q  ). 

The  effects  of  the  blade  sectional  steady  state  tip  speed  ratio  are  included  in  figures  70 
through  73.  The  lift  coefficient  is  assumed  to  be  directly  proportional  to  the  tip  speed  ratio,  as 
shown  in  figure  70. 

The  final  system  blade  spin  rate  is  proportional  to  the  blades  section  tip  speed  ratio.  All 
other  terms  remain  constant  except  the  deployment  angle,  whose  differences  in  values  reflect  the 
effect  of  the  gyroscopic  term  that  is  proportional  to  the  blade  spin  rate. 

3.  Effects  of  Blade  Sectional  Maximum  Steady  State  Lift  Coefficient. 

The  effects  of  the  blade  sectional  lift  coefficient  are  shown  by  varying  the  maximum  lift 
coefficient  value.  As  before,  the  lift  coefficient  is  expressed  as  a  linear  function  of  the  blade  tip 
speed  ratio.  The  performance  data  in  figures  74  through  77  show  that  increasing  the  maximum  lift 
coefficient  value  results  in  increasing  hub  rotation.  The  increased  hub  rotation  results  in  a  higher 
blade  spin  rate  and  a  slightly  larger  deployment  angle;  the  blade-to-hub  spin  rate  ratio,  however, 
decreases.  The  effect  of  increased  hub  spin  rate  and  lift  coefficient  values  results  in  increased  system 
drag  and  drag  coefficient.  The  hub  tip  speed  ratio  increases  with  increasing  lift  coefficient. 

4.  Effects  of  Blade  Sectional  Maximum  Steady  State  Drag  Coefficient. 

The  effects  of  the  blade  sectional  drag  coefficient  are  shown  in  figures  78  through  81. 
Increasing  the  sectional  ^  ig  coefficient  has  the  most  marked  effect  of  decreasing  the  hub  spin  rate 
and,  consequently,  the  h.  i  speed  ratio. 

The  consequence  of  this  effect  is  to  decrease  the  blade  spin  rate,  increase  the  blade  to  hub 
spin  ratio,  and  decrease  the  deployment  angle.  The  interesting  result  of  all  of  these  effects  is  that  an 
increase  in  the  drag  coefficient  of  the  blade  results  in  a  slight  decrease  in  the  total  system  drag  as 
shown  in  the  effective  drag  area  and  drag  coefficient.  The  change  in  system  drag  is  small  compared 
to  the  relatively  large  change  in  the  drag  coefficient  value. 

5.  Effects  of  Blade  Steady  State  Sectional  Lift  to  Drag  Ratio. 

The  effects  of  the  lift  and  drag  coefficient  can  also  be  expressed  as  the  ratio  of  the  lift 
coefficient  to  the  drag  coefficient  fi  e.,  the  lift-to-drag  ratio).  These  cases  are  shown  in  figures  82 
through  85.  Essentially  increasing  the  lift-to-drag  ratio  results  in  an  increased  drag  for  the  system. 
Whether  the  increased  lift-to-drag  ratio  results  from  increasing  the  lift  or  decreasing  the  drag,  the 
basic  result  is  to  increase  the  hub  spin  rate  (i.e.,  hub  tip  speed  ratio)  with  the  resultant  increased 
dynamic  pressure  acting  on  the  blade.  The  value  of  the  lift  coefficient  results  in  a  much  more 
pronounced  effect  than  the  value  of  the  drag  coefficient  for  the  same  lift-to-drag  ratio  because  of 
the  effect  of  the  value  of  the  lift  coefficient  in  addition  to  the  increased  hub  spin  rate. 

E.  Aerodynamic  Effects  Distribution  Over  the  Blade. 

Figure  86  is  a  graphic  representation  of  the  manner  in  which  the  various  aerodynamic 
terms  vary  at  different  spanwise  locations  (stations)  on  the  blade.  The  example  shown  is  the 
nominal  case  at  steady  state  conditions.  The  velocity  at  each  station  is  the  vector  sum  of  the 
free  stream  velocity  and  the  velocity  caused  the  blade  spin  about  the  hub.  Thus,  the  velocity  and 
the  angle  of  attack  increases  from  root  to  tip.  The  increased  local  velocity  results  in  an  increased 


Figure  75.  Effect  of  Blade  Sectional  Maximum  Lift  Coefficient  on  Steady  State 
Deployment  Angle  and  Blade/Hub  Spin  Ratio 
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Figure  82.  Effect  of  Blade  Sectional  Maximum  Lift-to-Drag  Ratio  on  Steady  State  Tip  Speed  Ratio 
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Figure  83.  Effect  of  Blade  Sectional  Maximum  Lift-to-Drag  Ratio  on  Steady  State 
Deployment  Angle  and  Blade/Hub  Spin  Ratio 
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Figure  85.  Effect  of  Blade  Sectional  Lift-to-Drag  Ratio  on  Steady  State 
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dynamic  pressure;  consequently,  the  aerodynamic  loads  are  larger  at  the  blade  tip,  as  indicated  in 
figure  87,  which  shows  the  aerodynamic  force  distribution  over  the  blade  span.  The  angle  of  attack 
of  the  local  velocity  at  each  station  determines  the  directional  sense  of  the  lift  and  drag  forces,  as 
shown  in  figure  86.  The  resultant  force  due  to  the  lift  and  drag,  shown  as  a  dashed  vector  in  figure 
86,  is  seen  to  act  to  spin  the  hub  up  at  stations  near  the  blade  root  and  to  act  to  de-spin  the  hub  at 
stations  near  the  tip.  This  is  also  shown  in  figure  87.  The  component  of  this  resultant  force  in  the 
direction  of  the  free  stream  velocity  is  the  system  drag  for  the  blade,  shown  in  figure  87. 

The  essentially  constant  nature  of  the  lift  coefficient  distribution  over  the  blade  is  shown 
in  figure  88.  Figure  88  also  shows  how  the  local  tlade  tip  speed  ratio  decreases  from  root  to  tip. 
This  results  in  the  blade  spin  moment  acting  to  spin  the  blade  up  near  the  tip  and  to  de-spin  it  near 
the  hub.  ' 


The  angle  of  attack  distribution  (including  the  effect  of  the  induced  angle  of  attack)  is 
included  in  figure  88,  as  is  the  true  yaw  angle  of  the  blade,  0,  which  is  used  to  determine  the 
aerodynamic  coefficients  values.  The  relative  tip  speed  ratio,  uR ,  which  is  a  function  of  0,  as 
indicated  in  figure  88,  is  also  shown.  The  relative  tip  speed  ratio  ojr  indicates  whether  an  overspin 
condition  exists,  and  this  also  determines  the  aerodynamic  coefficient.  It  will  be  noted  that  60%  of 
the  blade  length  is  in  a  true  overspin  condition. 

X.  DEPLOYMENT  TRANSIENT  CHARACTERISTICS. 

A.  Detailed  Discussion  of  Equations  of  Motion  Terms. 

Before  presenting  the  results  of  the  analysis  of  the  full  deployment  transient 
characteristics  of  the  rotor  system,  it  would  be  helpful  to  consider  the  meaning  of  each  of  the  terms 
in  the  equations  of  motion.  A  familiarity  with  the  basic  effects  of  each  term  is  helpful  in 
understanding  the  consequent  transient  (and  steady  state)  behavior  of  the  system. 

Consider  the  equations  for  the  three  degrees  of  freedom  of  the  system  shown  in  figure  1 0. 

1.  0 -Equation. 

The  equation  describing  the  angular  acceleration  in  the  0 -sense  includes  four  basic  terms, 
which  are  as  follows. 

a.  cjn/j  sin  0-An  inertial  moment  term,  which  could  be  termed  a  “gyroscopic” 
moment.  Because  an  autorotor  blade  spinning  with  a  positive  u,  as  shown  in  figure  10,  will  create 
aerodynamic  forces  that  cause  a  hub  spin  rate  in  the  negative  S2  sense,  this  moment  will  always  act 
in  a  negative  0  sense  (i.e.,  to  force  the  blade  back  to  a  smaller  deployment  angle).  Note  that  the 
term  becomes  larger  with  increasing  values  of  0.  This  moment  is  present  even  at  steady  state 
condition. 


b.  £S2  (70  -  /j)  sin  0  cos  0  -The  moment  resulting  from  centrifugal  force.  This 
moment  always  acts  to  pivot  the  blades  out  normal  to  the  hub  spin  axis.  It  has  a  maximum  value 
whcn0  =45°. 


c.  Me  A  -The  external  aerodynamic  moment  acting  on  the  blade  in  the  0- 

direction.  It  always  acts  in  the  negative  0 -sense  and  is  composed  of  the  blade  lift,  drag,  and  side 
force  coefficients. 

Note  that  all  these  moment  terms  are  divHed  by  the  moment  of  inertia  of  the 
blade  (/0)  to  give  the  blade  angular  acceleration. 
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Figure  87.  Steady  State  Aerodynamic  Force- Distribution  Over  Blade  Span  for  the  Nominal  Case 
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Figure  88.  Steady  State  Conditions  Spanwise  Distribution  for  the  Nominal  Case 


d.  M$- An  aerodynamic  “damping”  moment.  This  moment  is  proportional  to, 
and  acts  in,  a  directional  sense  opposite  to  the  rotational  rate  of  the  blade  in  the  deployment  sense 
(0). 

2.  fi-Equation, 

The  equation  describing  the  angular  acceleration  of  the  blade  about  the  store  center  line 
(i.e.,  hub  spin  acceleration)  is  composed  of  four  terms. 

• 

a.  -8coIx  sin  8- A  .“gyroscopic”  moment  that  is  effective  only  during  the 
deployment  transient  because  of  the  8  term.  From  a  consideration  of  the  terms  and  their  signs,  it 
can  be  seen  that  this  term  results  in  a  “spin  acceleration”  effect  on  the  hub  as  the  blade  deployment 
angle  is  increasing  (i.e.,  positive  0),  because  the  hub  spin  is  in  the  negative  sense. 

• 

b.  -2 S18(I0  -  A, )  sin  8  cos  0-This  is  a  Coriolis  moment..which  acts  to  decelerate 
the  hub  spin  rate  as  the  blade  deployment  angle  is  increasing  (i.e.,  positive  8). 

c.  -The  external  aerodynamic  moment  acting  to  spin  the  blade  about  the 
hub.  As  can  be  ^ n ,  a  positive  value  of  the  moment  results  in  hub  spin  moment  in  the  negative 
sense.  This  moment  is  a  function  of  the  blade  sectional  lift,  drag,  and  side  force  coefficients. 

d.  Ix  co  cos  0-This  is  a  kinematic  term  that  evolves  from  the  derivation  of  the 
equation  in  the  form  where  the  blade  spin  rate  co  can  be  considered  as  a  relative  term. 

3.  co-Equation. 

The  equation  describing  the  blade  spin  acceleration  about  its  long  axis  (i.e.,  autorotation 
spin)  contains  three  terms. 

a.  M  -An  external  aerodynamic  spin  moment.  This  net  moment  is  determined 

A 

by  summing  the  spin  moments  acting  at  the  various  blade  sections.  It  is  only  a  function  of  the  blade 
sectional  spin  moment  coefficient. 

b.  S18IX  sin  8- A  “gyroscopic”  moment  that  acts  to  accelerate  the  blade  spin  rate 
as  the  blade  is  deploying  (i.e.,  8  is  positive). 

c.  Sllx  cos  8 -A  “kinematic”  term  that  is  a  result  of  the  derivation  that  makes  co 

a  relative  term. 

B.  Relative  Contributions  of  Terms  to  the  Transient  and  Steady  State  System  Motion. 

The  relative  contributions  of  the  various  terms  in  the  equations  of  motion  to  the  motion 
of  the  system  during  the  deployment  transient  is  shown  in  figures  89  through  91  for  the  8,  SI,  and 
co  equations  respectively.  In  these  figures,  the  terms  representing  the  moments  resulting  from 
various  means  are  plotted  as  a  function  of  time  during  the  deployment  transient.  The  term  values 
used  for  the  system  are  the  same  as  the  nominal  case  of  the  parameter  study  of  section  IX, 
specifically,  run  1 .  The  nominal  hub  spin  rate  at  blade  release  is  750  rpm.  Note  that  the  nominal 
interference  factor  of  K  =  10  is  used.  The  motion  of  the  system  during  the  deployment  transient  is 
shown  in  figures  92  through  95.  The  complete  computer  printout  for  the  nominal  deployment  is  in 
appendix  D.  The  purpose  of  these  plots  is  to  show  the  relative  contributions  of  the  inertial, 
gyroscopic,  and  aerodynamic  moments  to  the  transient  and  steady  state  motion  of  the  system  so  as 
to  facilitate  a  physical  understanding  of  the  system. 

Figure  89  shows  the  three  moments  acting  in  the  blade  deployment  (0)  direction  as  a 
function  of  time  during  a  typical  deployment  transient.  The  moments  are  caused  by  centrifugal 
(inertial),  gyroscopic,  and  aerodynamic  forces.  The  aerodynamic  blade  damping  moment  Mg  is  also 
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I,)  sine  cose  =  Centrifugal  Moment 


Figure  89.  0-Equation  Term  Values  as  a  Function  of  Time  During  a  Typical  Deployment 


Aerodynamic  Moment 


Figure  90.  ^-Equation  Term  Vaiues  as  a  Function  of  Time  During  a  Typical  Deployment 
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Figure  91 .  co-Equation  Term  Values  as  a  Function  of  Time  During  a  Typical  Deployment 
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Figure  92.  Effect  of  Initial  Hub  Spin  Rate  on  the  Transient  Values  of  the  Deployment  Angle  During  Deployment 


Spin  Rate  on  the  Transient  Values  of  the 


ililtli 


Transient  Values  of  the  Blade  Spin  Rate  During  Deployment 


Figure  95.  Effect  of  Initial  Hub  Spin  Rate  on  the  Transient  Values  of  the 
Total  System  Drag  During  Deployment 


included.  From  figure  89  it  can  be  seen  that  the  centrifugal  moment  always  acts  in  a  positive  sense 
in  that  it  acts  to  pivot  the  blade  out  to  larger  values  of  8.  The  net  aerodynamic  and  gyroscopic 
moments  always  act  to  force  the  blade  back  toward  the  undeployed  position.  In  the  initial 
deployment  phase,  the  centrifugal  and  aerodynamic  moments  are  prevalent;  but  when  the  system 
has  assumed  a  steady  state  condition,  the  gyroscopic  and  aeroi-vnamic  moments  are  of  about  the 
same  magnitude  and  together  balance  out  the  centrifugal  moment,  thereby  holding  the  blade  at  the 
steady  state  deployment  angle.  The  sharp  increase,  decrease,  and  steady  increase  of  the  values  of  the 
centrifugal  and  aerodynamic  moments  during  the  initial  deployment  phase  reflect  the  corresponding 
fluctuation  in  the  hub  spin  rate  during  this  time. 

The  time  histories  of  the  moments  in  the  hub  spin  (£2)  direction  for  the  nominal 
deoloyment  transient  are  shown  in  figure  90.  The  moments  considered  are  the  Coriolis  (inertial), 
gyroscopic,  and  aerodynamic.  It  should  be  noted  that  because  of  the  defined  positive  sense  of  the 
blade  spin  rate,  u,  the  resulting  aerodynamic  moments  acting  on  the  blade  will  cause  the  hub  to 
spin  in  a  negative  sense  (as  defined  in  figure  10).  Thus,  a  negative  moment  value  tends  to  spin  the 
hub  faster,  whereas  a  positive  value  tends  to  resist  the  hub  spin  rate. 

From  figure  90  it  can  be  seen  that  the  gyroscopic  moment  is  small  relative  to  the 
aerodynamic  and  Coriolis  terms.  Thus,  the  hub  spin  acceleration  is  essentially  due  entirely  to  the 
aerodynamic  and  Coriolis  moments.  The  net  aerodynamic  moment  always  tends  to  increase  the  hub 
spin  rate.  The  large  Coriolis  moment  resulting  from  the  initial  rapid  deployment  motion  of  the 
blade  is  the  main  reason  for  the  abrupt  decrease  in  hub  spin  rate.  The  consequent  oscillation  of  the 
blade  deployment  angle  in  reaction  to  the  hub  spin  rate  change  results  in  the  Coriolis  moment 
fluctuating  between  a  negative  and  positive  value.  With  the  transient  effect  of  8  diminished,  the  net 
aerodynamic  moment  then  continues  to  accelerate  the  hub  spin  up  to  its  steady  state  value  where 
the  net  value  of  the  aerodynamic  moment  is  zero. 

In  addition  to  the  external  moments  acting  on  the  system,  the  moment  of  inertia  for  the 
£2 -equation  also  varies  during  the  deployment  phase,  becoming  larger  as  the  blades  deploy  to  larger 
angles. 

Figure  91  shows  the  moments  affecting  the  blade  spin  acceleration.  As  can  be  seen  in 
figure  91,  tne  initial  blade  deployment  motion,  combined  with  the  hub  spin  rate,  results  in  a  large 
positive  value  for  the  gyroscopic  moment,  which  acts  to  give  a  large  positive  spin  moment  to  the 
blade.  This  is  the  most  significant  effect  in  initiating  blade  spin  (autorotation).  The  aerodynamic 
moment,  acting  to  spin  the  blade,  increases  as  the  blade  is  deployed  to  larger  values  of  8 . 

As  will  be  noted  in  figures  89,  90,  and  91 ,  the  deployment  transient  effects  are  essentially 
overby  1  sec  jnd  after  deployment  initiation. 

C.  Illustrative  Examples  of  Deployment  Transients. 

In  order  to  investigate  the  general  transient  behavior  of  the  system  during  deployment, 
deployments  were  simulated  using  the  computer  program.  Selected  variables  were  varied  to  show 
their  effect  on  the  resulting  motion,  using  the  nominal  case  terms  (table  III,  run  1)  as  the  basic 
input.  Figures  92  through  95  show  the  effects  of  varying  the  initial  hub  spin  rate  on  the  system 
motion  as  a  function  of  time.  Figure  92  shows  the  characteristic  deployment  motion  of  the  blade 
where  the  blade  is  initially  thrown  out  to  a  large  deployment  angle,  hesitates  momentarily,  and  then 
deploys  to  the  final  steady  state  value.  Figure  93  illustrates  the  severe  hub  spin  deceleration  that 
occurs  during  the  initial  deployment  phase.  Note  that  the  hub  spin  raie  decreases  to  about  the  same 
minimum  value  regardless  of  the  initial  spin  rate  value.  Figure  94  shows  how  the  blade  spin  rate 
steadily  increases  throughout  the  deployment  transient.  The  total  system  drag  is  shown  as  a 
function  of  time  in  figure  95.  Although  these  data  are  for  nominal  conditions  selected  to  illustrate 


the  system  performances,  they  show  he  same  characteristic  behavior  as  was  obtained  experi¬ 
mentally  in  figures  43  through  45. 


The  effect  of  free  stream  velocity  (or  free  stream  dynamic  pressure)  is  illustrated  in 
figures  96  through  99.  It  should  be  noted  that  a  steady  state  condition  is  achieved  more  quickly  at 
the  high  velocities,  although  the  associated  spin  rates  are  higher. 

XI.  SUMMARY. 

Equations  have  been  evolved  that  describe  the  aerodynamic  and  dynamic  performance  of 
a  self-deploying  articulated  autorotor  decelerator  system.  These  equations  describe  the  transient 
motion  during  the  deployment  sequence  whereby  the  decelerator  system  goes  from  a  stowed  (low 
drag)  configuration  to  a  fully  deployed,  high  drag  configuration;  they  also  fully  describe  the 
motions  of  the  various  system  components  and  the  overall  system  performance  during  both  the 
transient  phase  and  the  final  steady  state  conditions.  These  equations  were  formulated  using  the 
modified  Euler’s  equations  approach  and  include  the  inertial,  gyroscopic,  and  aerodynamic  forces. 

The  aerodynamic  force  distribution  and  the  resulting  moment  acting  on  the  decelerator 
blades  are  determined  by  means  of  the  blade  element  propeller  theory.  The  sectional  aerodynamic 
characteristics  of  the  decelerator  blade  are  determined  from  wind  tunnel  tests  of  the  autorotor 
section  utilized  in  the  blade.  The  wind  tunnel  data  thus  obtained  can  be  utilized  directly  in  the 
equations. 


General  expressions  were  evolved  to  describe  the  transient  and  steady  state  aerodynamic 
characteristics  of  the  autorotor  configuration  and  were  verified  by  wind  tunnel  tests.  The 
expressions  derived  permit  the  complete  transient  behavior  of  an  autorotor  to  be  determined  with  a 
minimum  of  wind  tunnel  testing. 

A  condition  of  overspin  exists  at  certain  spanwise  locations  of  the  decelerator  blade, 
where  the  term  overspin  denotes  a  condition  in  which  the  tip  speed  ratio  of  the  blade  at  that 
location  is  greater  than  the  steady  state  tip  speed  ratio  corresponding  to  the  conditions  at  that 
location.  Wind  tunnel  tests  on  an  autorotor  in  an  overspin  condition  were  conducted  to  determine 
the  aerodynamic  characteristics  of  the  autorotor  for  inclusion  in  the  decelerator  system  analysis. 
These  tests  revealed  that,  for  an  autorotor  in  an  overspin  condition  (1)  the  lift  and  drag  coefficients 
remain  at  their  steady  state  values,  and  (2)  the  spin  moment  coefficient  follows  the  general 
expression  determined  for  the  underspin  condition. 

Wind  tunnel  data  obtained  on  a  nonspinning  autorotor  are  presented  to  explain  the  spin 
initiating  difficulties  characteristic  of  autorotors  in  general  and  to  illustrate  the  effect  of  spin  in 
increasing  the  lift  and  spin  moment  coefficients  of  the  basic  shape. 

Hie  complete  equations  were  programmed  to  be  solved  on  a  digital  computer.  In  addition 
to  the  transient  and  steady  state  performance,  the  distribution  of  the  various  terms  over  the  blade 
span  is  also  determined. 

Data  from  a  wind  tunnel  test  of  a  representative  full  scale  autorotor  decelerator  are 
presented  and  used  as  the  basis  of  comparison  with  the  analytical  results.  Excellent  correlation 
between  analytical  and  experimental  results  are  achieved  when  a  mutual  blade  interference  effect  is 
included  in  the  analytical  equations.  The  blade  interference  effect  utilized  is  in  the  form  of  a  local 
flow  directional  change,  similar  to  a  downwash,  which  would  seem  to  be  physically  valid. 

The  analytical  results,  which  were  verified  by  the  experimental  data,  indicate. the 
following  for  a  particular  autorotor  decelerator  system  under  steady  conditions. 
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Figure  96.  Effect  of  Free  Stream  Velocity  Magnitude  on  the  Transient  Values  of 
the  Deployment  Angle  During  Deployment 


Figure  97.  Effect  of  Free  Stream  Velocity  Magnitude  on  the  Transient  Values 
of  the  Hub  Spin  Rate  During  Deployment 


Figure  99.  Effect  of  Free  Stream  Velocity  Magnitude  on  the  Transient 
Total  System  Drag  During  Deployment 


1 .  The  blade  deployment  angle  is  independent  of  free  stream  velocity. 

2.  The  blade  spin  rate  and  hub  spin  rate  are  linear  with  velocity. 

3.  The  ratio  of  the  blade  spin  rate  to  hub  spin  rate  is  a  constant  with  velocity. 

4.  The  drag  coefficient  is  a  constant  with  velocity. 

A  parametric  study  is  included  in  which  the  steady  state  performance  of  the  decelerator 
system  is  investigated  analytically  by  investigating  the  effect  on  the  decelerator  performance  by 
varying  a  single  term  value  about  a  nominal  case  value.  Plots  of  these  effects  indicate  specific 
qualitative  trends  and  illustrate  the  contribution  of  the  inertial  and  especially  the  gyroscopic  effects 
on  the  system  performance.  The  presence  of  the  gyroscopic  effects,  which  are  a  function  of  the 
inertial  properties  of  the  system,  does  not  permit  the  non-dimensionalization  of  the  system 
performance  that  would  be  possible  it  only  the  aerodynamic  terms  were  used. 

A  nominal  deployment  transient  is  presented  in  detail  to  illustrate  the  relative 
contribution  of  the  various  terms  of  the  basic  equations  of  motion. 

Only  one  gyroscopic  moment  is  present  under  steady  state  spin  conditions;  namely,  a 
moment  that  acts  to  push  the  blades  back.  This  moment,  which  is  a  function  of  the  blade  spin  rate 
and  hub  spin  rate,  is  of  the  same  magnitude  and  directional  sense  as  the  aerodynamic  moment. 
During  the  deployment  transient,  only  one  gyroscopic  moment  is  of  consequence.  This  moment, 
which  is  a  function  of  the  hub  spin  rate  and  the  blade  deployment  rate,  acts  to  initiate  blade 
autoro.ation  and  is  the  major  r:<j3cn  for  the  autorotor  blade  self-starting  feature  during  the 
deployment  sequence. 
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Blade  aspect  ratio 

Angle  of  attack  for  nonspinning  autorotor 
Angle  of  attack  resulting  from  hub  spin  rate 
Flow  angular  change  (interference  effect) 

Distance  between  blade  center  of  gravity  and  hub  spin  axis-measured  along 
blade 

Yaw  angle  used  in  decelerator  system 

Blade  chord 

Drag  coefficient 

Steady  state  drag  coefficient 

Maximum  steady  state  drag  coefficient 

Lift  coefficient 

Transient  lift  coefficient  constants 
Steady  state  lift  coefficient 
Maximum  steady  state  lift  coefficient 
Spin  moment  coefficient 

Spin  moment  coefficient  evaluated  at  condition  cf  zero  spin  and  zero  yaw 

Transient  spin  moment  coefficient  constants 

Combined  transient  spin  moment  coefficient  constants 

Deployment  damping  moment  coefficient 

Drag  coefficient  based  on  swept  area 

Side  force  coefficient 

Steady  state  side  force  coefficient 

Maximum  steady  state  side  force  coefficient 

Drag  force 

Blade  drag  force  directed  along  wind  axes 
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Aerodynamic  force  acting  on  blade  in  the  direction  of  the  store  center  line 
Effective  drag  area 

Euler  angle  defining  blade  orientation  in  plane  containing  the  store  center  line 
Total  force  along  blade  spin  axis 
Total  force  acting  at  blade  pivot  point 

Aerodynamic  force  acting  normal  to  blade  and  in  a  plane  perpendicular  to  the 
store  center  line 

Aerodynamic  force  acting  normal  to  blade  and  in  a  plane  containing  the  store 
center  line 

Total  force  normal  to  blade  spin  axis 

Aerodynamic  force  acting  on  blade  referred  to  inertial  axes  system 
Gravitational  acceleration  constant 

Euler  angle  defining  blade  orientation  in  plane  perpendicular  to  store  center  line 
Angular  momentum 

Absolute  angular  momentum  referred  to  body  axes  system 
Moment  of  inertia 

Moment  of  inertia  of  blade  about  transverse  axis  through  blade/hub  pivot  point 
Moment  of  inertia  of  blade  about  the  blade  long  axis 
Total  moment  of  inertia  referred  to  body  axes  system 
Moment  of  inertia  of  hub 

Subscript  indicating  conditions  at  a  particular  blade  element  or  spanwise  station 
Unit  vectors  for  body  axes  system 
Blade  interference  factor 
Blade  length 

Distance  between  store  center  line  and  root  edge  of  blade-measured  along  blade 
Lift  force 

Blade  lift  force  directed  along  wind  axes 
Blade  mass 
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m 


MB  Spin  moment 

My,Ms,  M#  Total  moment  acting  on  blade-referred  to  the  Euler  axes  system 

Mu  ,  Me  ,  Mn  Aerodynamic  moment  acting  on  blade-referred  to  the  modified  Euler  axes 
system 

Afw,  Me,  Mn  Total  moment  acting  on  blade-referred  to  the  modified  Euler  axes  system 
M  Total  moment 

Mx ,  My ,  Mz  Total  moment- referred  to  body  axes  system 

Mj  Deployment  damping  moment 

S2  Hub  spin  rate 

A 

Q,  Ratio  of  tip  velocity  of  blade  about  the  store  center  line  to  the  free  stream 

velocity 

<P  Euler  angle  defining  blade  orientation  in  plane  perpendicular  to  the  blade  long 


Yaw  angle  used  in  wind  tunnel 
Ratio  of  circle  circumference  to  radius 
Dynamic  pressure 
Free  stream  dynamic  pressure 
Total  dynamic  pressure 

Distance  from  store  center  line  to  particular  blade  element-measured  along  blade 

Radius  of  curvature  of  S-section  autorotor 

Blade  station  located  nearest  root  of  blade 

Blade  station  located  nearest  tip  of  blade 

Blade  center  of  pressure  location  for  force  FN 

Blade  center  of  pressure  location  for  force  Ts 

Spin  rate  of  the  blade  relative  to  the  x-component  of  the  body  axes 

Blade  aerodynamic  reference  area 

Total  blade  area 

Spanwise  interval  used  in  blade  element  method 
Air  density 
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a 

Solidity  factor 

SF 

Side  force 

SFb 

Blade  side  force  directed  along  wind  axes 

t 

Time 

T 

Blade  thickness 

At 

Integration  time  interval 

0 

Deployment  angle 

V ,  V 

Total  velocity 

vm 

Free  stream  velocity 

wgt 

Blade  weight 

CJ 

Blade  spin  rate  about  the  blade  long  axis 

-  # 

CO 

Total  angular  velocity 

dj* ,  >  w^ 

Total  angular  velocity 

u'X,u'y,  (o'Z 

Total  angular  velocity-referred  to  body  axes  system 

0)X,Uy,  C02 

Angular  velocity  of  the  body  axes  system,  relative  to  the  inertial  axes  system- 
referred  to  the  body  axes  system 

CO 

Total  angular  velocity  of  the  body  axes  system 

co/fi 

Ratio  of  the  blade  spin  rate  to  the  hub  spin  rate 

A 

co 

Ratio  of  tip  velocity  of  blade  about  the  blade  long  axis  to  the  total  velocity  (tip 
speed  ratio) 

A 

w/c 

Ratio  of  actual  blade  sectional  tip  speed  ratio  to  the  steady  state  tip  speed  ratio 
corresponding  to  the  local  sectional  F.ow  conditions 

A 

wss 

Steady  state  tip  speed  ratio 

A 

CO  — 
jS  0 

Maximum  steady  state  tip  speed  ratio 

Angle  used  to  relate  the  wind  axes  to  the  inertial  axes  in  the  decelerator  system 

*,  ^,2 

Inertial  axes  system 

x,y,z 

Body  axes  system 

xu>yu,*z» 

Wind  axes  system 
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Subscripts: 

X,Y,Z 

x, y,z 

y, S,<t> 

&,  $,  w 
ss 
B 
i 

Superscripts: 


(  ), 


Term  referred  to  inertial  axes  system 
Term  referred  to  body  axes  system 
Term  referred  to  Euler  angle  system 
Term  referred  to  modified  Euler  angle  system 
Steady  state  condition 

Refers  to  aerodynamic  force  or  moment  acting  on  blade  in  the  decelerator  system 
Denotes  condition  at  a  particular  blade  element 

Denotes  absolute  value  of  angular  momentum 
Denotes  a  vector  (except  in  the  case  of  q,  S,  and  a) 

Denotes  the  first  derivative  with  respect  to  time 
Denotes  the  second  derivative  with  respect  to  time 

Denotes  absolute  value  of  term  expressed  in  terms  of  the  body  axis  components 
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APPENDIX  A 

COMPUTER  PROGRAM  DETAILED  LISTING  AND  SYMBOL  DEFINITIONS 


Onjnn 

1* 

C 

50 AO  PARAMETER  5TU0* 

ooioo 

2* 

f  « 

etAINF  J  wr5TF*>VflT 

nomo 

**. 

c 

PARCH  1969 

onmi 

1  3* 

0TREN6T0M  ARACH(20)  »V56R(20)  .OFN(  20 1  .MOT  ( 201  *"9 !  1 101  •  ALMMO)  .05 1 ( 

onmi 

5* 

I101*T0*T|10J*A98T||"| ,5FBT(10I  *F»I (10) *FNI (10) 

noios 

6* 

DATA  pAD.* .REV. SS/.O! 73532*2. 32. 173. 9. 5392963. 11*0.2 

oonn 

7t 

N*»Tl 

oono 

8* 

c 

iXOUT 

ooni 

9* 

READ  lOm.NRtlN 

00113 

10* 

10 

Rf»o  tana. wgt. mgt*.  »lto.vtm6,6m 

0"1?3 

11* 

nEAD  1000.*.IL.6.C.r'HO 

00132. 

12* 

READ  1000.AI0.AJl.AIu.6V.0V 

00131 

13* 

RE *0  t0n0.(*L65.C055.Cr66.C0S.068 

00160 

13* 

READ  1000.C80.C91  .C97.CL1.CL7- 

00157 

15* 

# 

8r*0  1000. C8V1 .CHV7.r8OT.CP00.KI 

00166 

16* 

RFAO  10O0.THfTA.THET* .OMG.OPGT .M 

001-»6 

17* 

RE *0  IOOO.VLIm.TO.XO.tc.ThFT to 

00203 

18* 

READ  1000.ALFT.RTHT.flP.OT.OT2 

0  *17 1  3 

19* 

. RFAO  1000. TIME  1  *  TIMC7 .TI8f3.T6TP.A5T9 

00222 

20* 

RE  «0  1010»H(..HP7.HM3.MO*.HRAH.tftfH.TVrL.TEH0.7L 

na?36 

1  7 1  * 

■  .  .. 

RrAO.JOOPt.AM.AK 

00731 

27c 

RRTHT  7070  ■ 

0070  * 

73* 

RRTNT  7080 

0  0235 

23* 

CALL  C6U*7( APACH.HSrR. APCH.0.V568* 11 .Mil 

00236 

75* 

IF(f  I0EN.FO.01GO  TO  l® 

0"760 

2** 

POTHT  20*0 

00752 

27* 

- 

CAU  X5UB7(HGT,0EH.ALT.0.®H0.I2.82) 

"0753 

28* 

15 

"RtNT  2DI0.VGT »HGTS.*l TO.VTNF.<*MO 

00262 

2** 

•RTNT  2020.5.C.AL.6.GH 

00771 

303. 

RRTNT  2030, ATM. Alt. ATO.SV.OV 

0030" 

31* 

PRINT  7n30.CL56.C056.CYS5.08G.oPnT 

00(07 

37* 

PRINT  2n5n.nf8.C8Vl,O8V2.C«0T.rNn0 

00316 

33* 

.OMGZflHGVRr*, 

00317 

33* 

RRTNT  7035.CPn.CP!  »C**?.CL1  .CL7 

00376 

35* 

PRINT  7n60. THETA. TMrrT.v.vT. CO  - .THF1 IO.RLI W 

00337 

36* 

PRINT  2*00. RCNt. OR.»N,Air.IL 

00337 

37* 

c 

INITIAL  rpNOITTONS 

00336 

38* 

I6T0P=n 

ori**  39*  Tin 

nrmn  *"•  I7=l._  .  ... 

0(1351  "  THFTAsTHf TA*»40 

00362  37t  GH=GH*RAD 

0  "35  3  .  J.V _ UL;YU«-t£OS  (GH I 

00363  M*  7Trvi»^*5TNtr,H» 

0"366  55*  •  A»US5=HGT5/6 

o itw  35*  _ A»*saJMS*/G  .  . 

00357  *▼•  WrU/RFV 

0056"  50*  WlT«l=WLlH/«»rv 

0Q36J - 39 t- - TSU . -  ■ 

note?  5n*  xrto 

00363  51*  ,  2-70  ' 

003*3  52*  ..  .  ALUAIUO-  .. 

003*5  i  53*  ,THFII=TH2TTn 

00765  53*  C  i»mtJ*ro3 

003&f_ _ S5*.. .  sa _ .nxjoo.ji/.ifn.  . 

0O771  5**  Tl55I*tGM> 

9"3T?  57*  T7ZC05IGHI 

.^50* 

003^3  59*  C7rC05«THfT*l 

00375  G"*  CHV=CHVI*C*‘V?*0«G*07/(?,*VI*(f  1 

on37f  .  61*__  _ jr.i4DC»L.Cl»^lG0.I0  50 

00300  62*  C*LL  C5U07(MC7#0CN*«I  T,  j#OM0*iT'*'#N?l 

0030J  67*  *0  OiOH0*VlNF**2/?, 

0O3"7  $•*.  C7:U0-*U 

003"3  65*  C«T0*6**.**?*C1**2/VT8F 

00*"3  65*  03rO« 

00*05.  67*- _ AL*^a»  .... 

0C3"6  '6**  06  =  0. 

003"7  69*  T0*ft=ni 

.0031.0  TO*  .  AH9:0. 

0031!  71*  ANOrO. 

0031?  77*  ARTtO. 

0Q313  .  733-  - AHTSFJO. 

00313  73*  <r*=0. 

003t5  75*  fj?0. 

00316  .  .76*  .  FN:0. 

00317  77*  FN6=0. 

0"3?0  73*  '  flrOIWT 

00321  .  79*.  —  lF(J*A»f,NNIfi0  TO  65 

00373  80*  RRTNT  2300*7 

00376  81*  *5  0"  75  TSI.HO* 

00331  92*  .  AL®HArATANUS6tO"Gl*Cl*P/¥WF» 

00337  83*  C5rSIN(AL®M*» 

0"33T  83*  C6-C0*( AL°HA> 

Q0333...  95*  .  „  eETA=ATAHCSeRTfl.-C2**Z«C6**?)/CC?*C61l 

0033*  86*  C3t6TN(3FTA> 

00336  67*  VB350»TIVT»(F**7*0«r.**7*Cl**2****?l 

0"337  88*  Q*:RH0»V8**?/7« 

00330  $6*  C3:08«C«0* 

00331  «  9«*  «66rW*C/l7.*V81 

0033?  9t*  A*CW=V8*C3/SS 

00333  92*  CALL  C«U37(  l*ACH.  US**'.  A8CN.1  .**69.  tl  .*1  » 

nn*3«  *3*  COrC066*r3*«7 

00*36  93*  «0=W$66*C3**2 

00*36  96*  CL7CL56*CCLl*WS5/W66'»*fL?*«65**?/«358**?*C3**?1 

00337-  96*  tF(HS6.GC.W*KL7CL66*C9**2 


Ml  RtPHOOUOBU 
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*  «**•*»& 


nO«M 

IT* 

mCL.i  ft.n.irt  ro. 

on»^l 

9** 

Tr  (AMC**.nc.l.  irtro. 

nnm 

90* 

Al  PMAT«l.r»HA*»4*»*l.*p*o 

00444 

inn* 

rsrSiw(  *L'»wi  1 

00*57 

101* 

CftcCOM  AI.PMAJ 

Omen 

10?* 

isootf  ].-c’’**?#Cft  **?>/(r?*ccn 

oo*R  1 

101* 

f*:lTM|4rrii 

0n4f  9 

1 04  • 

CM  -  IC-o.C4**-»*rHl  ♦W**/«S54*o»?*W*,ft**2/US5«**2*C4**?) 

OPlftl 

104* 

C7:CTft<*5iw»2, *«rT*i*** 

004*4 

ICC* 

0«T<  T  >:C0*o7 

00*45 

107*. 

ALBlt I) :rt *f 1 

nn«sr 

in** 

ftfft! <  n  roy.ri 

004*7 

IQo* 

14  II.ft7.TL)  At.  OH  1  >5*tOi7(WDP-l*U*«L8T(IU/ri«4T«N0O*1-lL> 

no,7j 

:in* 

CTH*rr?/Cl 

r«47? 

111* 

■J*,5”S0°7  {  1  .  *C7NT  **7*rft**?l 

004*3 

!!*• 

Cl irC7vT*04/50ft 

00*74 

113* 

C»'*=l  ./ft°5 

004*4 

114* 

0«im?M.oMTt*Cl?*rr*o&I»T»*rc-<F8!«Ti*Cl  |*CS 

0047C 

lift* 

T""!(  »»  r  At  MIT  »*Cl?*oft-nBT  m*o5-«T8I»  n*Cll*CS 

004?* 

ns* 

A-oifti :C“*C3*C 

0*i«O" 

117* 

f*Kn:DOM  n*fft*C‘»-ft4Rj(Ti.Cl‘»*ri-5F4lrn*cil*C5*C2 

0O40J 

ii«* 

fVJCI  m*r  !?*o** C1*AL4T<  M *C 1 1 *C7*08 II l » *Cft*Cl *5*81 1 ! » *C12*C 

OOcnj 

Ho* 

17-*.rBTI  T»  *o|l*r4*ci 

00407 

l>o* 

04:0**oft! r  f 1 

005*3 

121* 

•  t  4r  Al  4  ♦  Al  4  T  |  T  ) 

0  0<O4 

177* 

* 

0r:0S«o4i *  r 1 

nf)40r 

1>3* 

TOPQr TOP0*  TOOT  1 T ) 

00506 

1*4* 

A"4:»44*AMOf ft » 

00507 

1?4* 

4r4:S4B*5rqriTi 

oqc  |  a 

I’ft* 

FlTf  4.C  Aim 

00411 

177* 

F*»:FN*rwim 

ontp 

1?*. 

FN5:FN<*CMT  IT  1  *o 

one  I  T 

1?«* 

AMT-«M**04T It) *0 

0051* 

1  in* 

A**prA»*0*7oo»tt  )*p 

0051* 

111* 

A**T4Fr%HT*F*44  4T  f  T  I  •" 

OOejf 

|17* 

TMj.pr.NMK.n  to  74 

00*70 

HI* 

A|  PM*:At«»Mi/o*n 

op471 

1  J4* 

POINT  7l00,T,e,ve.Wftft ♦ORItt) *A14II I) ,051 f 1 1 * TOOIl 1 1 **«9I <T 1 tSFBI II 

005?1 

P«* 

1  1  trt,FNTm.AlPH« 

OO440 

1  1ft* 

75 

p?p*op 

00*47 

15-7* 

IF  |J*  *,r  .N»*tfiO  TO  40 

00444 

114* 

point  7?00tne#U4.0'-«Top0»AMB*4«’*  *fa,fm 

00*44 

11<»* 

40 

0**0: av*o5 

0044  T 

140# 

0*  :Clift*0*5v 

oo«eo 

1M* 

QM-pMr 

0044  1 

147* 

0**f  =  P"r  f 

oocc  7 

141* 

TH^THC*  t 

00547 

144* 

T 1 r 7M4*  A 

004  .4 

144* 

Vr  :w 

0  044  4 

144* 

M143VT 

00444 

147* 

71*51  -0PAC-D4)  *T?7A»*45C 

0O447 

144* 

r*t:XT*7TT*nT 

00470 

144* 

orzlX»!*t»!*n7*.5 

OP47| 

140* 

tr»*07 

00477 

1M* 

TTsttl 

00571 

IS?* 

7JI:l  -orAF*7l  ♦WGT6-o«‘*T11/A4*ft4 

00474 

157# 

2MS77*7rT.OT 

0O474 

154. 

07TIZ»!  *7T)*n**#5 

f)0*T5 

1SS* 

7t  7*  0* 

00477 

15ft* 

’*r?U 

OOfOO 

157* 

M7TAMB/A7  1-0M.TH.C1.OPT*''? 

nocot 

15** 

wrw*wi*n* 

nofOj 

15** 

IF  I T«  Gf  ,T  T*F3  MirWl  T** 

00ftO4 

ISO. 

iriT.or.TTPFllWltO. 

one  Oft 

161  • 

0*MS  f  -  4**0«AT|  *VIE*C2v?«*OM*TH*C7*r  I*C?—7H*VC  *4 11  *F1  «C**V*0* 

OOftOR 

16?* 

lSV*0V*C40o*C8 I/cATM/4**AI1*C?**?*At0*F1»*2) 

00507 

163* 

OPOrOMft.oMor.OT 

oncm 

164  • 

TM*-!  |:|  -AMT*CI->4T4r*c7*wC*/'N*AIl  *n  *0**  *2*C  7*C  1  *r2*C“0T  *TM*C8*  *81 

OQfi»o 

164* 

I54*0EB*9*OMG**2*C?l/ftI0 

OOf  1 1 

16ft* 

TMFM-7M4TT*7MrTT*0T 

OOftl  7 

167* 

0*MJ  1  TM4  T  1  *  TMf  TM«nT*,1 

onci  * 

168* 

TM*T  A  r  T  Wf*  4  «OTH 

004  1* 

164. 

TncTIr*tiEIl 

00414 

170* 

CMrATAP|*T/Tr ) 

on*  14 

171. 

At  t:ai*o-* 

00417 

17** 

lr  1  Tv4t  to  •< 

00471 

171* 

VT4Fr 40 QT  ftr*TT*  7I*o*  ) 

0067* 

174* 

45 

14  <  JFVn  .fo.oioo  to  00 

004  74 

17ft* 

IrIAL*.  14.  *570)4,0  to  174 

00474 

17ft* 

*JO  TO  4ft 

00477 

177. 

40 

Tr IT.64.T5TP1CO  TO  1*5 

00411 

174. 

45 

Tr**p* 

00*17 

174* 

ion 

roNTlwtjr 

00$14 

140* 

I4«T.C4.Tlo£l  )N**:V*? 

no4 14 

141* 

IF|T.l7.7T«r?|ro  *0  i?S 

0*14*0 

147* 

N*»;n“7 

0"5«I 

183* 

07:0*7 

00447 

144* 

1  ?s 

04 “T INC  *MT  *J 1 

004*1 

1 4  ft  * 

f JO;C0* ( *mo  T  A ) 

004** 

144  * 

OHCP:OHfi*orv 

004*4 

187* 

0w«4:ow6* |4*IL/i. 1.04/VT4F 

one  *c 

144* 

WP;tf *orv 

0  04*1 

184* 

TmFTO:ThE  ti/oao 

0  0C<o 

140* 

G“0:GM/o*o 
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oorr* 

131* 

oocs** 

10’* 

nn«<T 

H3* 

oorr* 

1M  • 

OBfiS* 

MS* 

oorrr 

1AF* 

0"F«7 

!•*• 

"OFF" 

1M* 

0««M 

M9* 

0OCC7 

200* 

0  "Of  * 

2M* 

?o-»* 

onr  «« 

207, 

0"«F« 

?0A* 

oocct 

?ns* 

OOFTfl 

20S* 

0"F  *1 

207, 

oort* 

20M 

"OF? 3 

200* 

0  OF*  T 

210, 

nncTT 

2p* 

nnnn 

212. 

ORTA  "* 

21T. 

0  OTA  * 

21*  • 

007*  £ 

7tS. 

0O7A7 

?1  F  • 

0"7«p 

217* 

OfJTC? 

71K. 

0"7<t 

?J«» 

O07^f 

2?  o* 

po7«;a 

77|. 

007FC 

7?T* 

"OTSF 

22  3* 

rn*«c 

22** 

27«i* 

0"7SR 

22F* 

nfl7*F 

77*, 

POTSc 

22*. 

0"7SF 

?2«* 

0O*<F 

730* 

nn7«fc 

731* 

007S7 

2*7* 

"07ST 

73** 

00**.  0 

71*  • 

O07CO 

?!*• 

0n*F? 

23F* 

OOTfit 

7*7, 

0O7G? 

73A* 

0  0715  7 

23*. 

0"7F3 

?*o. 

Ooti;7 

®AT* 

00-*F* 

2*o, 

007c, 

7*3* 

O07C5 

2**. 

OGTF« 

7*s* 

*»0*CF 

7*  R  » 

OO’FT 

?*7* 

"OTT** 

7AA» 

00771 

7*0* 

qo-»77 

750* 

00*7? 

25»» 

007*7 

7R-», 

00777 

2C7. 
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DEFINITIONS  Of  COMPUTER  TERMS 


SYSTEM  PHYSICAL  CHARACTERISTICS: 


MACHN  - 

MACH  NUMBER  OF  VELOCITY  COMPONENT  NORMAL  TO  BLADE  SPIN  AXIS  ff 

wss  - 

ROTOP  BLADE  TIP  SPEED  RATIO  C0RRESP0N0IN6  TO  MACK  NUMBER  OF 

AIRFLOW  NORMAL  TO  BLADE  ...  Q 

S5o 

W6T  - 

WEIGHT  OF  8UDE  (LBS)  ..>  w9t 

W6TS  - 

WEIGHT  OF  STORE  (LBS)  .  .  .  «s 

ALTO  - 

INITIAL  LAUNCH  ALTITUOE  (FT) 

VINF  - 

INITIAL  FREE  STREAM  VELOCITY  (FT/SEC)  ...  V_ 

m  - 

OEfSITY  OF  AIR  (CONSTANT  OR  VARIABLE  WITH  ALTITUOE)  (SLUG/FT3)  p 

s  - 

BLADE  PLAN  FORK  AREA  (FT2)  . . .  S 

c  - 

BLADE  CHORD  (FT)  ...  C 

L 

BLADE  LENGTH  (FT)  ....  t 

6  - 

LOCATION  OF  BUDE  NIOrOINT  FROH  HUB  SPIN  AXIS,  MEASURED  ALONG 

BUDE  (FT)  B 

GH  - 

INITIAL  FLIGHT  PATH  ANGLE,  OFF  HORIZONTAL  (DEG.)  .... 

IH  - 

MONE.'T  OF  INERTIA  OF  HUB  ABOUT  HUB  SPIN  AXIS  (SLUG-FT2)  ....  IH 

11  - 

MOMENT  OF  INERTIA  OF  BUDE  ABOUT  BUDE  SPIN  AXIS  (SLUG-FT2)  I, 

10  - 

MOMENT  OF  INERTIA  OF  BUDE  ABOUT  BUDE/HUB  PIVOT  POINT  (SLUG-FT2) ...  I„ 

SY  - 

REFERENCE  /J !£A  OF  STORE  -  BUTT  END  CROSS  AREA  (FT2) 

Sy 

DY  - 

REFERENCE  DIAMETER  OF  STORE-BUTT  END  DIAMETER  (S  INCHES)  (FT)  .:..dy 

CLSS 

maximum  blade  lift  coefficient  c. 

«0 

COSS 

MAXIMUM  bUOE  DRAG  COEFFICIENT  Cn 

»0 

CYSS 

MAXIMUM  BUOE  SIDE  FORCE  COEFFICIENT  ...  C. 

“o 

OHG  - 

INITIAL  HUB  SPIN  RATE  (REV/MIN)  .;..no 

OMGI  - 

INITIAL  HUB  SPIN  ACCELERATION  (RAD/SEC2)  ...  ii0 

DEB  - 

RADIUS  DISTANCE  OF  BUDE  PIVOT  POINT  FROM  HUB  SPIN  AXIS  (FT)... 3 

CMY1  - 

TURBINE  VANE  SPIN  MOMENT  COEFFICIENT  C 

evaluated  In 

teiw  of  one 
vane 

CW2  - 

TURBINE  VANE  SPIN  DAMPING  MOMENT  COEFFICIENT. .  .C. 

v2 

OCT  - 

BUOE  ROTATION  DAMPING  COEFFICIENT  IN  B  -  PLANE  .... 

% 

CMO  - 

MAXIMUM  BUOE  SPIN  MOMENT  COEFFICIENT  C, 

0 

CM’  - 

BUOE  SPIN  TRANSIENT  COEFFICIENT  C^ 

CM2  - 

BUDE  SPIN  TRANSIENT  COEFFICIENT  C- 

c 

Cll 

BUOE  Lin  TRANSIENT  COEFFICIENT  .» 

CL2 

BUOE  LIFT  TRANSIENT  COEFFICIENT  C, 

TH)  - 

INITIAL  BUDE  DEPLOYMENT  ANGLE  (DEG.)  90 

THI  - 

INITIAL  BUOE  DEPLOTMENT  RATE  (RAD/SEC)  .:..eg 

U  - 

INITIAL  B’.ADE  SPIN  RATE  (REV/MIN)  ...u0 

WI  - 

INITIAL  BUDE  SPIN  ACCELERATION  (RAD/SEC2)  ag 

CDS  - 

STORE  BOOT  DRAG  COEFFICIENT  ....  Cn 

J 

THI  I  - 

INITIAL  BUTE  OEPlOTMENT  ANGJUR  ACCELERATION  (RAO/SEC2)  60 

Nil M  - 

BLADE  SPIN  RATE  LIMIT  (REY/MIN) 

HINT  - 

0IS7ANCE  FRO*  HUB  SPIN  AXIS  TO  INITIAL  BLADE  STATION  FOR  EVALUATION 

OF  AERODYNAMIC  LOADS  -  MEASURED  ALONG  BUOE  (FT) 

•  r0 

•’  **'»»«*•  M 


*-■*  'r^VO  <SM8B«WBB“S!E,r  '*«' 


■  C*'*T*-~  .*v  '  V  rtVI  S'  '. 


OR  -  DISTANCE  8ETWIEN  BLADE  STATIONS  -  MEASURED  ALONG  BLADE  (FT)  ,.;.  Hr 

AN  -  NUMBER  OF  BLADES  IN  SYSTEM  . . .  N 

AN  -  INTERFERENCE  FACTOR  (OEG)  . ...  K 

IL  -  BLADE  STATION  AT  WHICH  LINEAR  TIP  LOSSES  INITIATE 

AEROOYNAMIC  CONDITIONS  OVER  BLADE  SPAN: 

NOTE:  THE  SUBSCRIPT  (I)  REFERS  TO  THE  PARTICULAR  BLADE  STATION 
NUMBER  FOR  WHICH  THE  TERN  IS  COWUTED. 

TINE  -  TINE  AFTER  STORE  Rt.ZASE  (SEC)  t 

I  -  BLADE  STATION  (INTEGER)  1 

R  -  BLADE  STATION  AT  WHICH  AERODYNAMICS  ARE  EVALUATED  -  MEASURED 

ALONG  BLADE  (FT) r, 

VB  -  TOTAL  VELOCITY  OVER  8LA0E  AT  STATION  I  (FT/SEC)  ....  7, 

WSS  -  BLADE  TIP  SPEED  RATIO  AT  STATION  I  u, 


OB  - 

LB  - 


BLADE  DRAG  FORCE  AT  STATION  I  (ACTING  OVER  A  Ar  SPAN 
AREA  AROUND  STATION  I)  (LBS)  0^ 

BLADE  LIFT  FORCE  AT  STATION  I  (LBS)  L„ 

■t 


OS  - 

TORQ  - 
MB  - 

SFB  - 

a  - 

FN  - 
ALPHA  - 
NOTE: 


BLADE  DRAG  FORCE  (IN  OlAECTION  OF  FREE  STREAM  VELOCITY)  AT 
STATION  I  -  (LBS)  ...0, 

SI 

BLADE  TORQUE  FORCE  (IN  HUB  DISK  PLANE)  AT  STATION  I(LBS)  .... 


BLADE  SPIN  MOMENT  AT  STATION  I  (FT-LBS)  ...  Nj 

BLADE  SIDE  FORCE  AT  STATION  I  (LBS)  ...  SF. 

"t 

BLADE  LIFT  COEFFICIENT  AT  STATION  I  ....  C. 

L1 


BLADE  NORMAL  FORCE  (DOE  TO  AERODYNAMICS)  AT  STATION  I  (LBS)  . ...  F„ 

"l 


ANGLE  OF  ATTACK  AT  STATION  I  (OEG.) 


THE  SITUATION  OF  THE  TERMS  AT  THE  VARIOUS  STATIONS  REPRESENTS 
THE  TOTAL  LOAD  ACTING  ON  THE  BLADE.  THESE  TOTAL  OR  NET 
LOADS  APPEAR  BENEATH  THE  CORRESPOWING  COLUmS  FOR  DB.  LB. 

DS.  TORQ.  MB.  SFB.  FA.  ANO  FN 


SYSTEM  PERFORMANCE: 


X  - 
2  - 
TH  - 
OMSS  - 
W  • 

XI  - 
21  • 
THI  - 
OMG  - 
WI  - 
XIG  - 
2IG  - 
THII  - 
OMGI  ■ 
Q  - 


HORIZONTAL  OISTANCE  TRAVELED  B’  STORE  (FT) ...  X 

VERTICAL  OISTANCE  TRAVELED  BY  STORE  (FT).:...  2 

DEPLOYMENT  ANGLE  OF  BLADE  (OEG)  0 

TIP  SPEED  RATIO  OF  BLADE  TIP  ABOUT  HUB  0 

BLADE  SPIN  RATE  (REV/MIN)  u 

HORIZONTAL  VELOCITY  COMPONENT  OF  STORE  (FT/SEC)  ....  \ 

VERTICAL  VELOCITY  COPTONENT  OF  STORE  (FT/SEC)  ...  2 
OFtOYMENT  RATE  OF  BLADE  (RAD/SEC)  .»  9 
SPIN  RATE  OF  HUB  (REV/MIN)  ...  B 
SPIN  ACCELERATION  OF  BLADE  (RAO/SEC2)...  i 
HORIZONTAL  ACCELERATION  OF  STORE  (j'j)  ..;.'x/9 
VERTICAL  ACCELERATION  OF  STORE  (9's)  ..,.'i/9 
DEPLOYMENT  ANGIAAR  ACCELERATION  OF  BLADE  (RAO/SEC2) .  ...  g 
ANG1AAR  SPIN  ACCELERATION  OF  HU8  (RAO/SEC2)  fi 
FREE  STREAM  OYNAHIC  PRESSURE  (IBS/FT2)  ....  q,  q_ 

TOTAL  AXIAL  FORCE  ALONG  BLADE  SPIN  AXIS  ACTING  AT  BLADE  BEARING  (LBS) 

TOTAL  NORMAL  FORCE  TO  BLADE  SPIN  AXIS  ACTING  AT  BLADE  BEATING  (LBS). 

MMua  T0  STWE  water  line  acting  at  blade 
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Fj  -  TOTAL  FORCE  COMPONENT  PER'EICICULAR  TO  STORE  CENTER  LINE  ACTING 

AT  SLADE  BEARING  (LBS)  ...  Fj 

F  -  TOTAL  FORCE  ACTING  AT  BLADE  BEARING  (LBS)  ...  F 

MACH  -  FREE  STREAM  MACK  NUMBER  M„ 

NK  -  TOTAL  (NET)  SPIN  MOMENT  ACTING  ON  BLAOE  SPIN  BEARING 

(FT-LBS)  ...  K 

U 

KT  -  TOTAL  (NET)  MOMENT  ACTING  OH  BUDE  SCARING  IN  DEPLOYMENT 

ANGIE  PLANE  (FT-LBS)  M0 


HO  - 

DRAG  - 
VINF  - 
ALT  - 
CD  - 
GY  - 
Df  - 
ATT  - 

AMO  - 
COT  - 

RBT  - 
R80  * 


TOTAL  (NET)  MOMENT  ACTING  ON  BLADE  BEARING  in  HUB/DISK  PLANE 
(FT-LBS)  Mg 

STORE  DRAG  OUE  TO  ALL  BLADES  (LBS)  ...  D 
FREE  STREAM  VELOCITY  (FT/SEC)  .>.  Y^ 

ALTITUDE  (FT)  ALT 

TOTAL  0RA6  COEFFICIENT  BASED  ON  STORE  OIAICTER  CQ 
FLIGHT  PATH  ANGLE  OFF  VERTICAL  (DES)  ...ry 
STORE  DRAG  OUE  TO  BOOT  (LBS) 

HOMENT  ABOUT  BLADE  BEARING  IN  BLAOE  DEPLOYMENT  PLANE  DUE  TO 
AEROOYNAMIC  LOADS  OH  BLADE  (FT-LBS)  ...  M. 

flA 

*A*DUtT0  •aamnc 

“a 

TOTAC  DRAG  COEFFICIENT  BASED  ON  TOTAL  SWEPT  AREA  ...  C_ 

CENTER  Of  PRESSURE  FOR  BLADE  DEPLOYMENT  MOMENT  (FT)  ....  r, 
CENTER  OF  PRESSURE  FOR  HUB  SPIN  MOMENT  (FT)  ..> 
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APPENDIX  B 


INPUT  CONDITIONS  AND  SOLUTION  PRINTOUTS  FOR  THE  ANALYTICAL 
SOLUTIONS  OF  THE  STEADY  STATE  CONDITIONS  FOR  THE 
FULL-SCALE  AUTOROTOR  DECELERATOR  IN  THE  UNIVERSITY 
OF  MARYLAND  SUBSONIC  WIND  TUNNEL 


Canputer  Pro yrm  Input  for  the  Analytical  Computation  of  the  Perforate* 
of  the  Full  Sell*  Autorotor  Decalerator  it  Tastad  In  the  University  of 
Maryland  Vfnd  Tunnel  : 

Ntr»M  V«* 

a""0"i)"O"  .10010000*01* 

•5"""0OAi«no  ,inoiiooo*oi 
.3i«ioon*H  .5«oiP001*iO 
.  crvi/w)no',»'v  .e’onooorucc, 

.t on ooo nr* #00  .53000000  30 
.•nn^oao't*"/*  .*iawooi*00 
.*i"000n*n  .?eonopi*00 
.  liyingnna.Ol  .?oonno01*rC 
•f*O"00lO*"l  .na^hOOn^CO 
.7HnoOiO*ni  .13000000*00 

wr.tr  .3*«omio*oo  v«73z  .*oOiocao*"*  »tro:  .30000001*03  vinf=_  .ioooooi0*O3  *"0z _ .?3T7qo'ti-p2 

t:  .2*7|*m0*03  C:  .25000000*""  l  r  .515*n0i0*f*0  «:  .7*30i""0*00  «Hz  .**0000000 

|MI  .*7500000-0?  ||Z  .305*ni"0-i*  TOr  .5t??f>Ci0-07  7V=  .13020300*00  "  5vz  .eiSloOO"* Otf“ 

CISC  .13511100*01  r«S$  =  .17510000*11  CTS<=  .*7500000*00  0?GZ  - ,7SOiOD«0*O3  0**1 .00000001 _ 

0C«=  ."1111000  C«V|:  -.3*3*5100-01  C*»2:  -.17**1 7*0*00  CW»:  -.$OOOeO"0*O?  C*C»Z  /* 00 00  000 

O' nr  .717000000  *00  f|:  .1 50°  0000*  ^  C»*z  -.  3*100000*00  Ptl  =  ,171**i00*0*  0.7=  ‘-TTlOOOff^OO  ■  ’ 

TMnr  . *«ionno0*0?  tmI  =  .loOiOOiO  v:  .••P0000"*0*  viz  .000010*1  COS”  .5O00O0W*«O 

TMTtz  .01501111  VLtNZ  .!**330000**^  ®»VTZ  .7lCI50"0*no  "•=  .51*00"00-01  OTKOz  •  •  100000 00* J2 

*v=  .*nooiOO«oi  a«z  ."0010000  ttz  ln 


Steady  State  Conditions.  K  •  0  : 
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.0*5* 320**00 

•»r= 

.0?I57*M*0" 

•noz 

•35315275-07 

Steady  State  Conditions.  K  •  10  : 


t 

• 

W 

V55 

05 

t.1 

05 

T"*e  "• 

3*0 

rn 

at>na 

*.155 -"1 

1.035*0* 

1.153*10 

3.«C4-0| 

*.750-m 

5. li"-0» 

1.553-OJ  -5.*5"-03 

3.000-07 

1  .*7?*00 

5.157-01 

25.71  - 

*.0*5 -m 

1 .05«*1* 

1  .1  3**"o 

3.570  -"1 

*.*•5-11 

5.*7*-01 

1. <33-01  -5.55«-07 

7.355-02 

1. *75*00 

5.5*0-01  *7.00 

*.107-11 

l.C«C.O* 

1  .105*1" 

3.071-11 

*.750-0* 

5.0*1-01 

1.7*7-01  -*.507-"3 

*.031-02 

1.575*00 

5.5**-01 

37.71 

*.*1*-"1 

1.110*0* 

1  -0  75^1f> 

*.0**  -"1 

5.057-"! 

5.173-fll 

1.015 -Cl  -3.775-03 

7. *00-02 

I. *07*00 

5.  *11  -01 

51.51  ' 

C. *35-11 

1.155*0* 

1  .0*7*10 

*. 355-"I 

5.*15-il 

5.070-0! 

5. <*5-0*  -I.7?a-03 

7.505-02 

1  ,*0Y»O6 

»V.M 

t.7S3-"I 

1.13«*0* 

l  .on? *on 

*.5«5-"l 

5. *05-01 

*.3*7-11 

7.535-0*  -*.55*-"5 

*.*55-02 

1. *51*10 

7.5*3-01 

55.15 

*.cTo-ni 

1. 7*0*0* 

5 .715-1 1 

5.051  -«1 

«.*$»-"! 

*.57n-oi 

-1.03*-0*  1. *37-03 

2.3*7-02 

1. *55*00 

0. 033-01 

51.15* 

e.*07-i| 

!. 7*7*0* 

5.357-11 

5.*50 

c.**e-"f 

0.771-01 

-*.*•7-0?  *.07*-03 

2.170-"? 

1.5*1*00 

O.**2-0l 

55.50 

1 .050*10 

1.377*0* 

5.011-11 

5.507-"l 

?.*15-"f 

5.2t*-0I 

*5. 0*0-0?  c. *30-03 

7.053-0? 

1.577*10 

5.515-01 

55.07 

H 

1.1*7*10 

1 .S5«*1* 

0  .«  50  -O  1 

f .1*5 --I 

7.0C*-Ol 

5.**«-0! 

- 1.57*-0!  0.5*  *-03 

2.015-0? 

1.35**0" 

1.001*00 

50.53 

*.«05*"0 

S.n*e.rvi 

*.775*00 

3.133-01  3. *5 *-07 

2.5?4-Ot 

e;»5T-oi 

"TV*  0*60 

fz  .701*7*75*07  7z  .DUJCIT.O*  »m:  .XC«I17*0?  -.l"0*5531*0l  v=  .5t«57?53*05 

ITT  .*01 3133**0*  71=  .<1  30*15*  *••»  =  .l*3|«63|-0*  0 "<  =  -.*7*07717*03  VI 7  .755a535*-0? 

in;:  -.*i«ient«*"0  »tt:  .3*013570*"**  tmfi:  .0"OOOCOO  0*7=  -.M3*?773-oi  •=  .11005000*07 

rtr  .3«<3«*l**n?  **'-  -.55  *5*** 75  *01  'I:  .  7??C«tr>*nt  r*  ?  -S**3*5"3*07  fz  .3a3031l3*07 

HCH7  .*0*05051-01  •»:  .***175*1-0*  "7=  .OidooOOn  -"Z  -.Cl 3*5015-0$  *•**:  .75105173**7 

*tvz  .inOOOOOO*nt  Ut;  .3517**33*0*  CO:  .t«5t51***17  <v=  .37*5**««*07  nr  7  .5710*037*00 

ITT:  .5**1 3*33*01  a»o:  -.  *3*n«5*7^|-*  C"T:  .3775155**00  »0T7  .70555*07*00  *50Z  -.  151 1557*-01 

HOI  S*P«00UC'8U 
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APPENDIX  C 

COMPUTER  SOLUTION  PRINTOUTS  FOR  THE  STEADY  STATE  PERFORMANCE 
PARAMETER  STUDY  RUNS  LISTED  IN  TABLE  III 


MJH  1 


I 

» 

v* 

9*3 

Of 

13 

0< 

TO*C  *■. 

5FS 

CL 

rn 

•  LPM0 

1  3.1  32-0  r 

l7o*2*02 

“  V.TW*oo 

3.21*  -0 1 

3.90  0- Of" 

3.700-01 

I. 690-0)  -5.720-P3 

7.375-02 

1.391*00 

3.026-01 

31.20 

2 

K.m-m 

1.665*02 

l.t  31**0 

3.3*7-01 

3.173-01 

5.11 *-01 

1.373-01  -3.330-03 

2.777-02 

1.303*00 

5.161-01 

75.33 

— 1- 

A.oco-m_ 

-I.fl33ttl2 _ 1.103*00 

3.551-01 

3.303-fl 

5.370-01 

1.15*-01  -7.336-03 

2.260-02 

1.307*00 

S.5Z3-0J 

30.67 

• 

5.030-01 

1.125*0? 

1  .071*00 

3.7*3-01 

3.673-oj 

5.073-01 

0.7? *-0?  -2.000-03 

2.105-02 

1.331*00 

5.372-01 

32.15 

5 

6.660-01 

1.1*0*07 

1  .033*00 

3.021-01 

3.30t-Ol 

6.305-01 

5.600-02  -1.9*3-03 

2.030-02 

1.353*00 

6.370-01 

35.32 

-I. 

1 .133*11? 

1  .005*00 

3.306-01 

S.J0O-O1 

*.77?-0I 

2.370-0?  5.231-05 

1.536-02 

1.330*00 

6.035-01 

39.39 

T 

0.336-01 

1.231*0? 

3.707 -ni 

0.625 -*1 

5.735-01 

7.773-01 

-1.053-02  1.737-03 

1.303-0? 

1.501*00 

7.355-01 

51.36 

9 

3.170-01 

1. *04*02 

3.3*1-01 

3.377 -nj 

*.330-01 

7.332-01 

-3.313-02  3.602-03 

1.777-0? 

1.536*00 

0.037-01 

53.33 

_ 9_i.0M*nn 

.L339*Q?  5*fll9rfll 

-3.3*3-01 

*.705rfU 

3 .330*01 

-0.733-0?  5.636-03 

1.710-0? 

1.521*00 

0.,533-ni 

56.66 

in 

1.093*00 

1.303*0? 

«. 703-01 

5.7*2-01 

7.037-01 

3.37  7. 01 

-1.733-01  7.035-03 

1.639-02 

1.335*00 

3. 033-01 

50.69 

3.230*30 

5. 33  5*  no 

f.6M*00 

3.219-01  -3.023-07 

2.033-01 

3.237-01 

6.76  3*00 

.  D  Ilf  ML  CHECK  f*UU  *T_.  OIMI0  ..  . 


X  = 

.27536669*03 

27 

.26330725*07 

TM  = 

.77312397*0? 

0*657 

-.3323?30?*OO 

9  = 

.72O65O0**O* 

—Xlz. 

.  ..2*023057*02 _ 21= 

..30037030*07. 

THI  = 

.-.1913Q3  36-03 

0"G: 

-.06332016*03 

VIZ 

-.75567630-0? 

JKz 

-.10300000*00 

716= 

.35037331*00 

TH 1 1  - 

-.33030336-10 

0*G  I  = 

.33376593-03 

Or 

.1 1085000* 0? 

r»: 

.57*67130*02 

FN  = 

-.61231335*01 

nz 

.6631 0397*01 

F7  = 

•  576r9t  72*0? 

Ft 

.57331977*0? 

«a tz. 

_  -.30903031-01 _ n*Z 

-.33333125-0* 

Jtr;. 

-.56033310-12 

■  r>r_ 

.5?5«5o«7-06 

O»06i 

.267*9397.0? 

fins 

.10000000*03 

•  Lt  = 

.23003775*07 

co= 

.00  0000  00 

CV= 

.1637*?7?*0? 

OF: 

•OOOOOOOO 

*TT- 

•52020761*01 

•  PK>= 

-.01009620-05 

COT  = 

.53539231*00 

00T  = 

.7571 ??7?40o 

•  "0  = 

-•259?7?99-P9 

I 


I  RUN  2 


..JL- 

J...  - 

-  Jfk.  _ _ Jf*S _ 

-  DB 

L3  _ 

05 

70*5  V 

SF6 

CL  . 

_ F* 

1 

7.332-01 

1.092*02 

1.163*00 

3.060-01 

3.739-01 

3.992-01 

1.612-01  -5.333-03 

9.099.0? 

1.91 1*00 

9. 612-01 

30.39 

_ 2-AJtfgfU _ 1 .0651112 _ 1 .139*00 

3. 203 -01 

3.379-nj 

9.905-0! 

.1.359-01  -0.629-03 

•  •.71 5-02 

1*917*0.0 

0,936-01  39,23 

7 

6.000-01 

1.032*0? 

1.110*00 

3.S31-OI 

9.209-01 

5.167-01 

1.030-01  -3.799-03 

9.596-0? 

1.923*00 

5.293-01 

37.95 

9 

5.039-01 

1.123*02 

1.073*00 

3.606-01 

9.071-0! 

5.532-01 

9.796-0?  -2.67*-03 

0.376-02 

1.930*00 

5.605-01 

91.39 

1.159*021- 

.1.096*00 

3.059-01 

0.779-01 

5.907-01 

5.336-0?  -1.930-03 

0.205-02 

1.937*00 _ 6.119-01 

99.63 

6 

7.5IW-01 

1.137*0? 

1.013*00 

9.135 -01 

3.127-0? 

6.39  7-01 

2.225-02  -3.101-07 

9.033-02 

1.995*00 

6.570-01 

*7.77 

7 

0.336-01 

1.230*0? 

5.706-01 

9.093-01 

5.517-01 

6.097-01 

-1.126-02  1.609-03 

3.960-0? 

1.952*00 

7.079-01 

SO. 6. 

_  -0 

.5.170-01- 

1.212*02— 

3.936-01 

9.737^*1 

S. 990-01 

.7.002-01 

-9.702-0?  3.333-03 

3.600-02 

1.950*00 

_  7.6 16 -01  53.17 

3 

1.000*00 

1.330*02 

9.112-01 

5.177-nj 

C. 005-01 

0.251-01 

-7.075-02  5.330-03 

3.362-02 

1.553*00 

0.503-01 

56.77 

1« 

1.009*00 

1.373*02 

0.707-01 

5.531-01 

7.103-01 

0.737-01 

-1.230-01  7.973-03 

3.299-02 

1.529000 

9.009-01 

50.77 

-  - - 

- - 

- 

9.127*90 

.5.101*00 

.6.360*00 

3.0*0-01  -1.955-07 

.9.092-01 

9.630-01 

6. 592*00 

- . 

1  = 

.25937031*03 

2= 

.19957115*0* 

TM= 

.71909155*02 

CW*5= 

-.9060*0 7? *00 

9  = 

.92503033*09 

JtU_ 

-Hismuaz 

--  U= 

•02070290*02 

*•  19155939.-09 

.  0-6= 

-.000009 50*03 

Ml  = 

-.760267*9-02 

*16  = 

-.26353025*00 

216= 

.92119*03*00 

TMTfz 

»T9??9?7S-10 

omxz 

•  95  36  S'*  5* -03 

0  = 

.110*5000*02 

F0  = 

.29909116*02 

rnz 

-.29670191*01 

Fir 

.63597057*01 

F?= 

.2*961709*02 

F  = 

•29553519*02 

NCM  = 

.9O3O5031tQ1 

T.22313099H16 

*T  = 

.22737367-12 

■  nr 

• 3069007?-05  0006= 

.  .25039193*02 

710  = 

.10000000*03 

•  LT  = 

.310*2005*0* 

C0= 

.OOOOOOOC 

677 

.29971990*0? 

PF  = 

.00000000 

ITTr 

.51097770*01 

01*0= 

-.09569907-0* 

COT  = 

.59676090*00 

wr: 

.75597507*00 

•90  = 

-.16309135-09 

KM  3 


.  I 

* 

M 

0*5 

OS 

IS 

05 

TP50  *« 

SF0 

a 

FN 

0LPH0 

i  T 

3*332-01 

1*092002 

f.TsiToo 

37279 -ni 

-9. 066- Of 

0.906-01 

1.729-01  -5.050-03 

1.005-02 

11576*00 

0.921-01 

37.50' 

2 

9.166-01 

1.066*02 

1.129*00 

3.930-01 

9*250-01 

5.239-01 

1.960-01  -5.029-03 

1.305-02 

1.512*00 

5.257-01 

35.32 

.  T 

_S  JMML-OL  _J.C9J*0Z _ 

1 .059*00 

3.S15.-0J 

r*.953-oi 

5.605.-01 

1.107-01  -9.070-03 

t. 326-02 

. 1.S10*"C 

.5  **20-01 

,3".96 

*  9 

5*039  “01 

1*125*02 

1.06**00 

S.033-OI 

9.753-01 

6.006-01 

9.901-07  -2.936-03 

1.269-0? 

1.516*00 

6.033-01 

0**.02 

5 

6.660-01 

1.161*02 

1 .036*00 

9.005-01 

5.066-01 

6.092-01 

5.910-02  -1.070-03 

1.211-0? 

1.510*00 

6.973-01 

95.67 

_1_ 

7.507-01  1. 700*07  1  -0*7*00 

_9 .371-03 

5.920-01 

_E.3IJ-01 

2.520-02  7.296-05 

1 .155-0? 

|.5?0*Q0 

6^999-Cl 

90*73 

7 

9.336-01 

1.292*02 

9.679-01 

0.690 -"1 

5.936-01 

”  7.51’-0I 

-6.959-03  1.707-03 

1.005-0? 

1.553*00 

7.556-01 

5!  .90* 

0 

9.170-0! 

t*207*O2 

9.339-01 

5.093-01 

6.2S7-01 

7.599-01 

-9.770-02  3.675-03 

1.092-0? 

1.532*00 

0.037-oj 

50.33 

— 3_l«flnft*flg. 

t.*  JIA*02._ 

9* COS-01 

5.929*01 

6*659-01 

0.997-0! 

-9.037-02  5.726-03 

1.001-02 

1.509*00 

9.535-01 

56.55 

10 

1*009*00 

1.399*02 

9.602-01 

5.099-01 

7.009-01 

t.991-01 

-1.369-01  7.992-03 

5.693-03 

1.905*00 

9.099-01 

50.59 

9.363*00 

5.397*00 

6.009*00 

3.296-01  -7.351-06 

1.199-01 

3.639-OI 

6.099*09 

7=  .23796276*03  2= 

•161033?S«O9 

TH= 

.00393329*02 

0S*5= 

-.99359006  *  00 

wr 

.91910779.0* 

11  = 

.32000797*02  71= 

. 77?76500*O? 

THJr 

-.1957*169-0* 

0*r.r 

-.0550O«9O*O3 

WT  = 

-«62*7607*-O| 

ilZZ- 

-.26091076*00  _.2I6= 

•^7112966*00 

THltr 

.519603*7-09 

006 1  = 

.0056  300*  -0? 

.«  = 

.11005000*0? 

r*= 

.11*13*62*09  rnz 

-.12933063*0? 

FI= 

.6006*9  09  *  01 

F9r 

.11*60775*03 

F  z 

.110009*0  *  03 

*CM  = 

.90*09091-01  *9= 

-.73209*35-05 

0T= 

.59639390-06 

.10057025-03 

0*66= 

.27233990*07 

-Hi: — 

.  10000000*03. 

__.3l9166T2*a= 

CO: 

.00000000 

67= 

.29575163*09 

OF = 

-OOOOOOOC* 

»TT= 

.52227953*01  0*0= 

-. 10*00355-09 

Cf»t= 

.59315739*00 

*6T= 

.75553?I0*OO 

9*0  = 

-.331950.“/-O3 

I 

* 

79 

955 

09 

I* 

06 

7090  *• 

5*8 

CL 

=  5 

H*NI 

.  1 

ITfV^T 

T.i  0«*"oo 

7.3*09-01 

2.912-0! 

3.591-01 

1.051-0!  -J.909-03 

1.699-02 

l.9tS**'0 

3.59n-**l 

39.00 

2 

3.792-01 

1.000*0? 

1  .1  ?1  *00 

?  .062-01 

3.053-01 

3.021-01 

0. *33-02  “3. =10-03 

1.599-0? 

1.000*00 

3.763-01 

39.10 

.  3 

*.375-01 

1.116*02 

.  1.093*00 

-2*590*01 

3, 21?-oi 

9.071-01 

7.003-0?  -2.609-03 

1.529-02 

1.991*06 

0,006-01 

•  *.25 

9 

9.9M-01 

1.107*0? 

1  *069*00 

2.795-01 

*3.009-01 

0.391-01 

5.730-02  -1.603-03 

1.970-07 

1.999*00 

9.269-01 

99.23 

5 

5.597-01 

1.106*0? 

1  .033  *00 

2.915-01 

3.615-01 

0.632-01 

3.293-02  -0.975-00 

1.011-02 

1.997*00 

9.561-01 

97.0s 

_ 6. 

—  6.1  26rfll —  1.717*0? 

.1.002  *oo 

_  3. 109.-01 

3.00S-D1 

9*90 7 -01 

1.100-02  1.077-00 

1.350-02 

1.300*00 

9.053-01 

99.73 _ 

7 

6.709-01 

1.256*07 

9.715-01 

3.312-01 

9.106-01 

5.?7*-Ol 

-1.030-0?  1.211-0  = 

1.296-0? 

1.503*oo 

«.179-01 

5?.?3 

S 

7.791-01 

1.296*0? 

9.909-01 

3.530-01 

0.509-01 

5.729-01 

-3.126-02  ?. 009-03 

1.215-0? 

1.503-00 

5.911-01 

35.00 

...5 

_  J.*75l-01 

. 1*335*02 

..9.106-01 

.  3 . 7B0 -0 1 _ 

.9.795-01 . 

6.000-01 

-5.052-02  =.705-03 

1.17J-02 

i.52S*00 

5.920-01 

56.97 

10 

9.050-01 

1.309*0? 

0.011-01 

0.090-01 

9.992-01 

6.367-01 

-0.790-02  5.000-03 

1.135-02 

1.503*00 

6.239-01 

50.70 

3.0*5*oo 

3.990*00 

0.0*0*00 

1.790-01  -2.552-07 

1.392-01 

5. Olt-Ol 

9.793*00 

-  oitiof-cnccK-Fmii.it- 

01  can.  . 

— 

xr 

.23291595*0= 

2  = 

.12232520*00 

THr 

.77303507*0? 

0*55  = 

-.9»*6=O<l*00 

9  = 

.05; 76510*00 

_ xl=_ 

-59121695*0? 

21= 

.75236296*0“* 

t*I  = 

-.10173350-00 

0*G  = 

-.11073100*00 

91  = 

-.7607*535-02 

016  = 

-.10162970*00 

216= 

•61633307.00 

THtI  = 

-•*6270910-06 

0*6  T  = 

.90500535-03 

0  = 

.H#053OO*oj 

F»: 

.50303655*0? 

TMZ 

-.71535907*01 

Ft = 

.00003776.01 

F  »: 

.505S096|.o? 

r= 

.50.72010.02 

-  Mftt- 

.JU=_ 

-*.11  266  799 -0£ 

sir 

-.19917010-00 

*0= 

-2900*027-05 

Q*6C= 

.L9S21095*?? 

tin: 

.10000000*03 

ILT 

.37767003*0* 

Cor 

•00000300 

,  r 

.301616*0*07 

OF  r 

.oooooroo 

m= 

*30759776*01 

mO= 

-.03231990-05 

roi= 

.»1«?»JST«00 

9*T  = 

.61000919*00 

o»0  = 

-.?0009»l«-o* 

NOT  REPRODUCIBLE 
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RUN  5 


I 

8 

V8 

V5S-. 

OB 

-  19 

OS 

7090 

SPB 

CL 

FN  . 

ALPHA 

1 

3.458  -01 

1.033*02 

1.173*00 

4.043-01 

5.014-01 

5.973-01 

2.810-01 

-7.919-03 

4.208-02 

1.859*00 

5.889-01 

29.15 

-2 

—4.542-01 

1.057*02 

iU47*Q0 

4 . 240  *0  1 

5.257-01 

6.436-01- 

2.049-01 

-6.9’6-03 

4.033-02 

1.063*00 

6.332-01 

33.46 

X 

5.625-01 

1.085*02 

1  .1  1.  »00 

4.489-01 

S. 567-01 

6.953-01 

1.675-01 

-5.642-03 

3.859-0? 

1.468*00 

6.829-01 

37.57 

4 

6.708-0! 

1.120*02 

1  .082*00 

4.792-01 

5.942-nj 

7.526-01 

1.281-0! 

-4.067-03 

3.687-02 

1.473*00 

7.379-01 

41.45 

7.792-01 

1.158*02 

1  .046*00 

5.148-01 

P.384-01 

8.156-01 

8.595-02 

-2.202-03 

3.516-02 

1.478*00 

7.985-01 

45.10 

6 

8.8  75  -01 

1.202*0? 

1  .009*00 

5.558-01 

6.891-0! 

8.844-01 

4.041-02 

-5.26  7-OS 

3.946-02 

1.483*00 

8.647-01 

48.50 

7 

9.458-01 

1 .248*02 

9.707-01 

6.020-01 

7.465-01 

9.590-01 

-9,045-03 

2.378-03 

3.177-02 

1.488*00 

9.365-01 

51.66 

a 

.1.104*00 

1.239*02. _a.131.--Ql.  _6.£36 -"1 

_8  422-QL 

1 .065*00 

-5.534-02 

5.075-03 

2.926-02 

_l.S51*0a 

1.034*00  AS.  16. 

9 

i.?!2*on 

1.352*02 

8.962-01 

7.105-01 

8.960-01 

1.137*00 

-1.164-D1 

8.056-03 

7.808-02 

1.523*00 

1.108*00 

57.53 

10 

1.321*00 

t. 408*02 

8.606-01 

7.7?7-Pl 

9.526-01 

1 .212*00 

-1.860-01 

1.130-02 

7.699-02 

1.493*00 

1.181*00 

59.68 

5.566*00 

6. *83*00 

8./62*00 

4.990-01 

-5.246-07 

3.426-01 

5.858-Ql 

.  1.570*00 

xr 

.20104146*03 

2= 

.  15700431  *0-* 

TH2 

.75529745*02 

GH«6r 

-.10317025*01 

W  = 

.97578031*04 

Xtr 

■  L!4?.75CLfi*02_ 

— 2U; _ 

^55819966  ♦O’* 

THL; 

-.12630048-44 

o**g; 

-.740254 86*03 

HL= 

t. 75761836-02  _ 

XK.r 

- .26784?  32*00 

?1G  = 

. 16259466*00 

THIIl 

.11091399-09 

OHGTr 

.95205728-03 

02 

•  1  !8050‘T*o2 

F»: 

•64  ?162n0*02 

ruz 

-,75227«62*01 

rxr 

.87620223*0! 

F  3r 

•64058876*0? 

F= 

.64655336*02 

HCRZ 

-.90909091 -01— 

HI,; _ 

-.53882043-06 

.  "T  = 

.10231815-11 

HOr 

.88358769-05 

09A6  = 

.35040090*02 

VTn: 

.inooaono*03 

«L1  = 

.34299069*03 

COr 

.00000000 

GVr 

.17105265*0? 

OF: 

•ooooouoo 

ATTr 

.81959686*01 

*  HO  “ 

-.12407079-0* 

COVr 

•53020114*00 

RBTr 

.00341644*00 

9B0: 

-.24863491-04 

_t _ 2 _ 


.  ¥2 _ '£.S _ _  O3 


I  *. ’32-01  1.P4?*0’  1.150*oo 


3  K.noo-'lt  t.P°3*C’  1.0®’*0o 

%  6.434-m  ).t’4«Q’  1 ,06«*o« 

_ 6  _ _$ *5£|<_-X»1  .1 , ico*Q’  l ,n  j j  ,o « 

6  7. *02-01  1  •  l  o<3-*0'>  B.qWf 

7  t  •  ?4  )  *07  9.C64-"! 

-J  _  *  .  *  ?n^l_.l.?«6*Q2_AU32U:n 

®  1  "0*00  ». 3*3*0’  8.4*9-°l 

1«  1.684*"’  1.343*0’  P.'Cr.oj 

nrvTnr"  o<r~K  t“*t“  oif074’‘ 


?.?17-’|  ’.’qq-’l 

■».44n-’l  5.0^6^’ 

3,7»'5-^i 

**  .7*7  }  ’.nn:."!. 

?.44’-n]  3.K4  8-01 

3.1<4  jii  . 

,3 ’8  -2J~_  4.162.-21  _ 
4.406-’! 
4.666-m 


7-01  1.168-01 

n-01_  9. *70.-0 2_ 
5-"l  7. ”6-0’ 

o-oi  *.0*6-02 
j.^pu-o"* 

o-O)  1.7*6-02 
"-"I  -5,747-03 
5=Q1  -3.24°-" 2. 
•-0*  -6.10’-02 
’-01  -9.1<6-0> 
0*00  2,272-01 


•7.679-03 
£.256-0?  . 

■1.266-** 

•5.*?*-P4 

3,q’4_P8 

0.051-04 
1.5»  *-01 
2.5*’-0’ 
3.547-03 
-2.734-06 


.409-03  1.532*00 

.221.-03 _ L.  5 13  ♦ID. 

.03P-rt3  1.534*00 

.855-03  1  .53»l*oo 

♦  677-03  .  U53»*00 
.604-03  1.536*00 

.  *26-03  1.547, HO 

.192-03  l_.522*n0_ 

.0*7-03  l.S00*00 

.95,-03  1 .477*0** 

.624-02  1.827-01 


3.329-01 
.1*563-01  : 

3,404-01 
4.070-oj  i 

4,3.63-0!  i 
4.680-nj  4 
6.C70-01  ! 
6.,  350-01  ! 

S. 6*2-71  « 
6.074-01  « 
4,588*on 


_ 1Z _ *IZ69.1527.*II2. _ Z2 _ .rjL?'’r}Q*77 _ TH=  k£.3£’33?1*02. 

»T  =  .4*650757*07  71  ,  1«  47*456*0-  Tmt-  - .  1 5 1<»1 64 *-04 

Xl"  =  -,l**?<»»73*Pn  7  rgr  . « 7«,cn«qc.o-  r».  ifir  -.*’1362414-04 

_  rjU. _ .37.99452VO? _ JFR7-  .^2n122r. 9^*(J.Z _ rv-  _  * 95 70 *$c?*"l 

"CH  =  .*0*04091-01  MVz  -^’SS-’O’-O*  **Tr  -.4’|”7S«-P6 

¥TK  =  .inono000*03  »tT=  .  1»*6o4! -.o-  C0=  .n"onOOUf) 

_JLlLr_. .  ^3AJJ^a?a*oi _ LS0-=_.-..7Xl4  4''l’-.0*_  fnj-.  .39280236*00 


RUN  7 


_ L__ 

_ O _ _  _  y  0 _ 

•»-  * 

04 

1 

9J- 

1 

X.XJZ-OI  l. 0*7*0’ 

.  .1  49*0" 

4.446-91 

5.S1 

6.673-01 

,  fc=ai  _i .  ifs.'iL 

.  u i’a*"o 

.4.6*5  -’1. 

5.76Q-C1 

_  T. l?o-91 

3 

♦;,onn-«1  1.093*0’ 

l 

4 .841  -4  ) 

6. 964-01 

7.6M-01 

M 

Cc«38.0f  1.174*0’ 

1  .965.90 

5.1*0-41 

6.4”-’! 

8.14f»-01 

_  5 

6.658--OJ.  1 . 1 r  0*02. 

-l.tTAJ**'? 

9  .5 14  »*  1_ 

-“.731:01 

6 

■*.5f*2wpl  l,loq«n’ 

9.993-oj 

4.4*3  — " 1 

.’oe-°i 

9.352-01 

7 

1.336-nj  1.241*97 

9.«  54-01 

6 . 5 1  7 -’  1 

7.431-01 

1 .094*00 

- 8 — 2_»J  7Q---TU —  1.236*112-  9V*  17-^1 

6. 7*5 -'•l 

_5  .I2L-R1 

_l.Qf 6*00 

4 

l,3’3*p? 

8.98"-") 

7  .  ?47 1 

4,7*7-01 

t.l’4*00 

10 

1.084*09  1. 3*3*0’ 

8.rf3-«i 

7.36fc_-,i 

9.  *1  7-of 

i.2"*’*on 

— 

..  . 

«  f  8*’ 

-7*215*on 

9. 17?*°0 

X  =  .10764746*07  zr  .  ?n 344*0* *0’  T«r  .•4^20434*02 

— JLlt-_^...36^inZ5a*.Q] _ Z.;___.*L’G5.6’9S*a_1  THI7_  ^-J3*Q7627-04. 

Ti«r  -.13570762*00  2tG  =  ,  o’oj  j  rj  7-n  1  *MlIr  .  1”7  7S9* -03 

6»r  .7*74?Voi,*n?  ruz  ,16u9’475*o»  ry-  ,Qi7l836**ni 

?:«:  ^-»qo-iQ?nai^.nL _ 1<*4 77046-0"  \  2LrJ  .  . ’TS46’Il-0S 

VIN  =  .10000oon*nx  *LT  =  . 7®61 5’0S*n’  "o=  .monOOOD 

»7Tt  .6°S01 34®*m  •  M0  =  -.10544481-0-  "0T  =  .74443446*00 


0Hf<: 

-.99161901*05  W: 

.  ,1  373*41’**5-  .. 

0,J  9  T 

-.84714’n4*0’  wir 

-.704621 c5-"l 

9MGT- 

.703?2’iq-0’  Or 

.1  184609O*"’2 

n-_ 

• *?7"8®95*C2  F- 

-^7*46033_*(>2  _ _ 

mnz 

'.7081 6* ’2-04  D*»G: 

.183f74"3‘*>;?  *” 

ovr 

.71691547*0’  OF: 

."OOOOOno 

R?  T  : 

.«44’49?*oa  -?A.a= 

--*72.9207*1-03 _ 

ho 

SfB_  _ 

CL 

FN 

*LPHX. 

2.*4’-01  -1.063-0’ 

6.817-03 

1.5*6*00 

6.676-01 

31.72 

1.074-0!  -4 .pj 5-0 3 

6.521T«3 

1.537*90 

7.124-01 

55*52- 

1.69’-01  -7.21(1-03 

6.737-03 

1.537*40 

7.6*18-01 

39.15 

1.7)1-01  -5. OK4-03 

5. 956-03 

1.538*00 

8.157-01 

4-.S8 

7,954-0?  -2.PJ15-03 

,5.680-4* 

i,5JA*no 

8,74 /-0| 

3.534-02  ).6’9-04 

5.411-03 

1.539*40 

9.375-oj 

48.86 

-l.l®’-0?  3.23’-03 

5.140-03 

1.540*00 

1.096*00 

51.69 

r6*534-02  6.S*«-03 

,_4*.939-P3. 

Z.J520*aCL. 

_A*QZn*00- 

?*,14_ 

- 1 • ’24-01  1.0’5-n? 

4.745-03 

1 ,49°*o« 

1.1’6*00 

56.38 

-l.H7*-OI  1.42O.0? 

4.566-0? 

!.47S*on 

1.204*00 

58.43 

4.4C6-21  -’,*55-06 

.5.,c.0Z-°2 

’94.-01 

9.184*"0, 

- - 

OP**r  - .991 1 34®4  *00  V  =  .6863681**04 

.  r\»r.z  -,84*9*60^*03 _ V\  -  1 976?98’-<M 

OHrr  =  ,c?ini’cj_o?  0:  .I188*0n"*’2 

Vz  . 7*744484.07  F=  .75801417*02 

“0=  ._  ,10270*19-03.  Jip_*Q r  ^36617347.0^ 

ov=  .*5O4n47*0i  of  z  .ooobonoo 

B4Tr  .75374775-no  o°0  =  -.23760853-03 


RUN  8 

I  *  78  w*$  08  18  5*  TOPQ  ~H*  "  SFB  CL  FN  TtPHA 


_  J  _ 3^. ^32_-0 1  _ _ 5 . U 2*_0J__ I. * 1 56 * P _* . 0 4JJ-0 2__9,9 70 -07  ^197-01 

7  4.  T66H>  1  9.327*01  1.131*00  8.817-02  1 .049-01  1.279-01 

3  5.000-01  5.965*01  1.103*00  1.878-02  1.101-C1  1.369-Q1 

— 9  S,A241fll_  1.0  71±fl  0__9.AZJ.-0Z-  _L.U3r  01  .  _  U  86  9-01 

5  6.648-01  5.802*01  1.039*00  1.005-01  1.287-01  1.576-01 

6  7.502-01  5.997*01  1  .005*00  1.077-0!  1.335-0!  1.693-lJl 

7  -8*3 J6_-01  6L.20ft*fl!  _1..7.a7-0Jl_  _1U56^01_. ZaSJArCl  t.81 8-OS 

•  “  770-01  6.8S2*01  9.388-01  1.798 -^tl  1.5(5-01  1.985-0/ 

9  1.000*00  6.670*01  94)39-01  1.391-01  1.677-  01  2.107-01 

_  1G — L*fl8tf-00—  f^AlAittl-_t_aHJL-dOZ  1.9.96  =OX_l*778rOL  _  2.233tqi 

1.079.00  1.336*00  1.673*00 


.  8. 284-02  1.9 30-03  6.J88-C3  1.8oi*00_  1 .207-01  31  .20 

3.582-02  -1.233-03  5.99  3-03  1.88'8«J0  1.290-01 '35.03 

2.898-02  -9.897-08  5.701-03  1.487*00  1.382-01  38.69 

2.181-0?  -7.001r08  ...5,862-pJ  J, 991900  1.983-0X  97.15 

t. 820-02  -3.658-P8  5.728-03  1.899*00  1.593-01  95.82“ 

6.078-03  1.301-05  8.990-03  1.898*00  1.711-01  88.49 

“2.649r03.  8,35.7-04 . 8.758-OJ  X-S0l*00  1.8V*M)1_5K36_ 

-1.103-02  9.008-04  8.442-03  1.586*00  2.009-01  58.98 

-2.188-02  1.809-03  9.275-03  1.521*00  2. 173-01  56.66 

-3.335-02  _  1.959-03,  9, 12QtoJ  .  1.89.5^0 _Z*Z8.trO! 
8.086-02  -7.361-07  5.110-02  1.059-0!  1.691*00 


,  .1: _ 

.19029092*03 

.2: 

.35641003*03 

TH= 

.-77312570*02 

QHS6= 

-.99294665*00 

W  = 

.46033008*04 

XT: 

.32614803*02 

71= 

« 13957520*0X 

THT  = 

-.14915292-04 

0H6  = 

-.43196709*03 

«I  = 

-.12790693-01 

xiGr  - 

-.J80M476-0I 

216= 

.83711050*00 

THItr 

-.66907714-04 

OHGT= 

.18577709-0? 

0= 

•297l?5U0*01 

_ FAr... 

-144172Jatfl2 _ TICZ. 

-.1531 089(1*01 

.16727944*01. 

n- 

.14401524*02 

.  F= 

.15490350*02  . 

HCHr 

.45454545-01 

HW= 

75055789-06 

HT  = 

—• 3840^^/8-06 

M0= 

.11J15271-04 

Dp  *G= 

.SS51|T9««ai 

VlHr 

.50000000*02 

*LT= 

C0= 

. 0  yunoouo 

6vr 

.13152416*02 

OFr 

.00000000 

..AIJ;... 

,13.90516  3*  QL_ 

AHO.=  . 

-.127156034-04, 

-COT=  . 

^•35538873*00 

38T= 

.  .75,712305*00 

_?  *0 = 

r*.l  57523-94jJJ3__ 

Appendix  C 


174 


II 


■T  ^  W  f*  .  *~- 


RUN  9 


T 

0 

V6 

V^S 

0? 

L6 

06 

1 

3.337-01 

?. 088*0? 

1  .161*00 

1.293*00 

J. 60^*00 

I .980*00 

? 

8.166-01 

?. 136*0? 

1.1  75*00 

1.355  *of) 

1.681*00 

?.07«*00 

3 

6.000-01 

?.t98*0? 

1.105*00 

1.832*00 

1.776*00 

2.277*00 

8 

5.538-01 

?.?6Q*P2 

1.0  73*00 

1.523*00 

1.859*00 

2.800*00 

6. *68 -01 

?. 338*0? 

l.0»5*"0 

1.628*00 

7.019*00 

2.588*00 

6 

7.502-M1 

7.818*0? 

1  .008*00 

1 .787*00 

?. 167*00 

?. 787*00 

7 

8.336-01 

?.50?*0? 

9.6  98-0] 

1.880  ♦'*0 

7.830*00 

3.077*00 

8 

9.170-U1 

2.595*0? 

9.386-01 

?.0?7*P0 

?.57S*00 

3.272*00 

9 

1 .000  *  00 

7.697*02 

9.005-01 

? .  1 08  »*n 

?. 729*00 

3.879*00 

to 

t. 088*00 

2.796*02 

0.5  73-01 

2.363*00 

1 .788  *01 

2.890*00 

2.176*01 

3.692*00 

2.763*01 

7080 

r» 

SFB 

CL 

FN 

ALPHA 

6.89«-Ql 

-2.392-0? 

9.626-Q2 

1,883*00 

1.9J5*nO 

31.55 

5.832-01 

-2.061-02 

9 .235-02 

1.886*00 

2.060*00 

35.88 

8.728-0! 

-1.653-02 

8.889-02 

1.888*00 

2.198*00 

35.16 

3.568-01 

-1  .169-02 

8.866-02 

1.803*C0 

2.362*00 

87.67 

2.725-01 

-fi.!0»-01 

8. US 8-02 

1.896*00 

2.58-1*00 

85.97 

9.939-02 

2.325-08 

7. 7J5-02 

1.500*00 

2.732*00 

89.07 

•1.787-02 

7.272-03 

7.183-02 

1.567*00 

3.015*00 

57.59 

•1.888-01 

1.506-02 

6.C59-02 

1.583*00 

3.209*00 

55.01 

•3.629-01 

2*356-02 

6.596-0? 

1.518*00 

3.811*00 

5T.23 

■5.535-01 

3.278-02 

6.3S2-02 

1.892*00 

3.618*00 

59.25 

1.316*00 

7.679-07 

7.993-01 

1.708*00 

2.705*01 

« 

..8*692713*01 

2= 

.15609717*0^ 

TH: 

.77078606*02 

0065  = 

-.10181813*01 

w: 

.18526598*05 

XT:  ' 

-.1*051329*07 

71  = 

.12610600*0’ 

TMT: 

.13709337-00 

0*G: 

-.17676839*08 

MI  = 

.15158137-01 

iKr 

-.72960702*01 

716= 

.25190670*01 

thtt; 

.79’262?6-10 

0  MG  I : 

-.18818609-02 

0  = 

.8  7539  9  99  *  02 

FAS.- 

,*20052992*0.3 

FN= 

-•  2523 ! 587*02 

Ftz 

.27537082*02 

F?= 

•  28027738 1()3 

F.= 

.281 88969*03  _ 

HCMI 

.18181818*00 

HW  = 

.00007916-06 

MT: 

.05078730-12 

mo: 

-.11228839-^8 

DRAG: 

.11013233*03 

.20000000*03 

«LT= 

.30090683*0* 

CO: 

.09000000 

6V: 

.173501 *5*07 

OF- 

.n00O0O0n 

* 

ATT: 

.71615013*0? 

•*o= 

.16900086-00 

COT: 

.61170677*00 

RBT: 

.7579791^*00. 

000  = 

.  .12151701^00  .  . 

RUN  10  . 

"  T 

* 

v*  6  0* 

t’ 

06 

TOO'' 

MQ 

9C0 

~cl  fn 

JL  3«.?J^1.. 

?  0 ,166-01 

3  c .nno-ni 

Jl.  «4«3»-oi  . 

*  fi.cGR-n^ 

6  7.*n?-P1 

_ X^fi-Jlfiirai  . 

8  *.T70-O» 

*  t.nrtn*nn 

_ln~  *  »QM  *ESL_ 


i  •  08  1 .1  53*<'o 

».'U6*o**  l.l‘***o* 


)  .04340* 

1  .125*0*! 

i  .?pn**? 


I  ,«  o«»pn 
l 

j  .n^.r** 
1  .(1324"i 


1.2iZiC2._9.67?-01 
i.?*7*o?  p.3'*'i-oi 

l.37*,*n-»  4.046-*! 

1  .!•%♦«»  3  .c',2-n1 


I --I 
I .7^7  -*1 

l.1*!*-** 

i.*»8«-*| 

7.!R8-«1 

?  .*5 ’4  •*•  1 
*\  7  «  8  -1 1 

■7  (  «•»«  .7  t 
7.1*8  4*0 


_'*#03F-P| 

*•  .  1  71-*  f 

?.?4*-'M 

**.  3CO_(\ f 
?,S»C.P| 

?.7l3-oi 

2*.V?-01 

3.’T8-*M 

*.6i4.gi 

2.702*04 


?,*<6-41  .  8,6S6-P2_-?,9?«* 03  __7. 7(9-03  1.510*90  2.863-01  XI. 50- 

’.R’l-ni  7.7P*-0?  ->.518-03  c. 968-03  1.512*00  ?. 632-PI  3S.32 

7.«np-oi  S.Q!*-0?  -?."14-oJ  6. *72-03  1.518*00  2. 817-01  38.96 

’.10F-P1  8.66*-*?  -1.870-03  6.379-03  „l,*I5.*nfl _ J*.020-;:.J_87.J£_ 

3#?7 r -O t  ?.*10-0?  -7,7*P-"8  s. 092-03  1.518*00  J.28«-«f  85.67 

*. 8*1-01  1.7*1-02  3.637-05  5.810-07  1.520*00  3.879-01  88.73 

. ’.7<l- Pi  -3.827-03  o.p83tH8  .5,o57-OJ,_l..S5.3*a?_ *JM_ 
3t9»w*01  -2.7*5-02  1.870-03  «. 739-03  l.53?-'*0  8. 028-PI  r0.33 

4.?»t-01  -4.c>4-02  ?. 867-0’  5.037-03  1.509*00  8.?73-oi  56.55 

0.*P1-*I  -6.8?7-n?  J,"76-P7 .  8,850-M _ 7,885*00  8.53o-ot  58. S* 

’.408*00  1.685-01  -8.959-06  5.977-02  1.825-0!  3. '26*00 


HZ 

^3l597R62*o* _ ZZ 

.2^26**09*0- 

T 

.-.0oTSq,9f.O? 

o»r  6= 

-•99?61080*OO. 

W: 

,-*5W1.77«*0J.  _ _ 

XT: 

.6*071  77]  ♦"•» 

’i: 

,  1P61P008.P’ 

TM»  ’ 

-. ’90500 ?n-04 

0“G: 

-.85*‘0*OO6*O7 

Ul: 

-.00*11107-01 

Xir.r 

-.17RJ7l73.no 

ZTG  = 

*MTT  = 

0HT. »: 

.17*06698-01 

9  = 

.$9099000*01 

FA: 

.56907137*22 

-c 

-.'•■»‘»7  0  7/\9.01 

rxz 

.  *00*31  J»  *01. 

F7= 

_.,S7?’40*e*O2_ 

F: 

_...5Z3258A8*D2 _ 

MTH  r 

.00009001-01 

h«: 

-,9«o**r0*5-(lr 

MT; 

. ’50*  9766-06 

**o: 

.769*0003-0* 

0*A6: 

. 1 3637207*0? 

VTN: 

.lOOOnoOO.OT 

«LT: 

.?77*»05!*f>- 

CO: 

.p^noonoo 

0v: 

„tt  !0Plt8*0*» 

nrr 

.oonooooo 

ATTZ 

^76107276*91. 

AMOZ.. 

3n??6  300  -9** 

..cor: 

.50719866*00 

09’  = 

.7£558798*90. 

.P  foz.. 

09202826 -oj.  _ _ 

•  RUN  11 


..i 

0 

_  U£6 _ 

— .0’-  .. 

L9  - .  . . 

0<  TP?*  ***  . 

5F«  . 

..CL,  _  .. 

FN 

ALPHA 

1 

’  •  ’  7?  — °  1 

1  ,00‘»*O'» 

1  .1R740" 

6.177  -’1 

T. 597-01 

*.0°r-oi  3.?’7-m  -1.060-07 

9.791-02 

1.011*00 

9.235-01 

30.39 

--  2- 

_8-t6«-"t 

l.pr  «*0? 

1 .L7**O0 

6.8-»5  -rj 

7.96  7-** 

9. 621.-01  2.77?-U1_-?,268-pt_ 

_9.*5I.t0.2.. 

_uaiit£s>_ 

••.’j 

3 

6.000-0] 

1 .097*C7 

1  .1  lO*on 

6.7*0-'') 

8 .*  t  8-P] 

1.031*00  2.1*9-01  -7.8*7-03 

9.118-0? 

1.023*00 

1.060*00 

77.90 

* 

6.838-01 

1.1?3*O0 

1  .079*nn 

». 719-01 

8.952-01 

!, 100400  1-66] -0]  -6. 367-03 

8.770-02 

1.030*00 

1.138*00 

81.38 

-5 

6.668  -21  ..  JL-.15.5tQ2 _ L.9 

_7.716-«1 

_9.558r/?!_ 

1.191*J39_  1^069-01  -2>fK3-o3 

■  l.?24»0P  «».H 

6 

1 .197*97 

1  .ni»4"-» 

0. 7’9 -” 1 

1 .077499 

1. 7*1*00  0.060-0?  -8.619-07 

8.386-02 

1.800*00 

1.317*00 

07.77 

7 

•.7.76-Pl 

] .^*0407 

9.70G-O] 

0.400-9] 

1 

1.37*400  -7.706-02  3.?2o-«3 

7.700-02 

1.051*00 

1.918*00 

50.67 

..  *_ 

«  ,J  70-2!. 

.  1x23.3*02  _ 

9  ,0 ,8/v-o  ] 

7.607  -9J_ 

ta2J*9o, 

-S.T-fcnz  5. -"<•-"5 

1.068400 

1 .6?$.Pt? 

57.37 

9 

1  .non* 00 

!. 370*07 

9.1  l3-nl 

! .0’?  *0 9 

1 .763*09 

1  .6C 2*09  -1.56«-01  1.06®-02 

6.70O-C2 

1.557*00 

1.703*00 

36.77 

10 

* .es**po 

1.770*0? 

P .7R7-01 

1.11940C 

] >*784*9 

1.709400  -2.087-01,  l.*06-n? 

6.611-02 

1.520*00 

1.803*00 

58.77 

_ _ _ 

_  _ 

*. 7C7.00 

,l.Qr:iQL-  1.  27  .••JlI_6..U'-Li?)_-\.?5il-U.7_5,  207-01- 

.»*Ul!i!U _ 

r: 

i7«o*]4n» 

7  = 

.  1)7*9^3402 

0M64: 

-.98600706400  Hr 

.97580 200*^8 

V  T  - 

.1  TR  *50  69*0*. 

_  u- . 

»*51«9]9.0’ 

THt  = 

-.10758818-00 

..8  ".L-l?®  ?P2PT. .  wi.=_ 

_i*.U50iRitr.v_  ... 

XTG: 

—.8*707**2—01 

’TG: 

?7n?8956409 

TH  3 1 Z 

-.7«’?0276-|O 

OHGT: 

#9061*001-03  0= 

.23800009*07 

fa: 

.SHH7in.ni 

fn: 

-.6®3Tt,:604?1 

CX: 

.1777*]]n*02 

F*: 

.57771773*0?  F: 

•591097’3*02 

."CH7. 

.00*04091 -?1 

-  “X=  . 

-.09  3.07000 

»T: 

r*  7*1*- 87^7-12 

wnr 

.  .509?Oi.9O-O6. po  86  =_ 

,.50936907*0?  __ 

V  IN  r 

.10*000*940  » 

UT: 

,*n?*7'>8n.0* 

*0: 

.rtPPPPOOO 

GVr 

.00100607*01  OF: 

.noooooo* 

ATT: 

.107700304*? 

*  **0: 

*7*C76T*-fj6 

r*T: 

.6»67R5?9400 

prt: 

•759«86?0*OO  P*0: 

-.15355739-00 

. 

RUN  12 

I 

P  V8 

ws 

Z  08 

t  n 

05 

TORO 

srs 

Ct  fN 

ALPHA 

1  7.732-01  l.08?*C?  1.130*00  3.221-01  '3.995-01  8.825-01  1.787-01  -2.201-03  7.781-02  1.888*00  8.763-01  31.22 

2  8.166-01  1.065*0?  1.105*00  3.372-01  8.107-01  5.t*2*01  1.887-01  -1.915-03  2.287-02  1.887*00  5.090-01  35.08 

3. S.DOOzOLl 1.053.102  _i«D.mQO  _  3.556  r01__4. 81 5x533,-01.  _  1*2 13 -01  tl4VJ.tlL3-._2*l9J=fl2 — 1*8911*00  _  S. 950^01  _36. 70 

9  5.839-01  1.125*0?  1.097*00  3.779-01  9.680-01  5.939-01  9.756-02  -1.139-03  2.101-02  1.993*00  5.895-01  92.16 

5  6.668-01  1.160*0?  1.015*00  8.026-01  8.992-01  6.381-01  6.358-02  -6.288-08  2.009-02  1.997*00  6.275-01  85.83 

_  ,6 I-x502-fll —  1 .199t02-  9.819-01  9^3lO,-jU...5^3>S-Jl — £*MA-IU— J I*18tta2._=l* 960=115 — 1.919-02  —1.500*08  .  i«l«gfll..9aJ»Jfl_ 

7  8.336-01  1.281*02  9.88S-01  0.629-01  5.613-01  7.275-01  -6.889-03  6.321-08  1.856-02  1.870*00  7.188-01  51.03 

8  8.170-01  1.286*0?  9.158-01  8.981 -'ll  5.816-01  7 .601-01  -8.982-02  1.392-03  1.814-02  1.819*00  7.501-01  53.16 

_ 9  1.0Q0*Q0_.  1.334*02  ..6‘328-OJ  5.366 tOI  ^.1131  ^01^4.1116^01-73.516-02  28 3r0 LJ48 *211  _  1^461^ illxll-. 

10  1.088*00  1.383*02  8.510-01  $.785-0]  6.256-01  8.800-01  -1.830-01  3.205-03  1.796-02  1.319*00  8.291-0]  56.90 

8.302*00  5.132*00  6.603*00  3.389-0!  -1.991-06  2.008-01  8.261<-0l  6.992*00’ 


X: 

.26570680803 

2= 

•  1*552839*0* 

THr 

.77565051*02 

0H6Sr 

-.99235535*00 

HZ 

•89985036*08 

XT  I 

.3051 37?0<  02 

71  = 

.79501170*02 

THI: 

•10598123-00 

0M6: 

-.86257108*03 

Ul  = 

-#32337699-01 

TIG? 

.2 62935  5Z*0  P. 

_.Z16=. 

_  .3512*152*0", 

_TH1X= — 

=.2«14tairfli. 

.  fl«5I=.. 

.  ,.llQ7J125z9r. 

..  s- 

-j,uuzsmi9z 

FA: 

.57657088*02 

fn= 

-.69293659401 

Ft: 

.67032850^01 

F?z 

.57*98076*02 

F: 

.57862053*02 

MCH: 

.90989 091-Ot 

mh: 

-.18975730-06 

«T: 

-.13901523-05 

HO: 

•11636781-09 

DRAG: 

•26912980*0? 

a ua=. 

j  m  ua  . 

. ..  JJA47>6iiM _ £0= _ »tt9aimap.« _ 

_ 

Of-- _ .gtUBiM"  _ 

ATT: 

.50958062*01 

AMO: 

-.96007766-06 

cor: 

.59600970500 

AST: 

.75097795*00 

«2&- 

“#  ?**6?277-08 

1 


NOT  reproducible 
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Appendix  C 


m 


I 


“  *v»  v>  .%•  .3  *  ‘•■*  *.^t’  -VV*-*  'fH'tfi*  ^’-V'^*'  *>V^4  ’V-SS*  *T*  J^*\\*^,^*vj?:f'^‘'£‘’'^' 


I 


I 


RUN  13/ 


T 

9 

V8 

v<5 

08 

19 

1 

7.732-01 

1 *08! *02 

8.9  71-0! 

3.290  -o  1 

9.018-01 

» 

9. !*6-ni 

f.Ofi  7*0? 

8.?8*-ni 

3.986-^71 

8.196-01 

7 

8.000-01 

1 .0°n*0’ 

8.091-01 

3.5*9-01 

9.919-«1 

8 

,5.839-01 

i.l’ixi’ 

7  .869-° 1 

3 . 7  78  -o  i 

8.650-01 

5 

".668-OJ 

1.155*0’ 

2.6  38-Oj 

8.0i6-O| 

9.980-01 

6 

T.*02-oj 

!.!®'*0’ 

?.39?-Oi 

8.291  -oi 

5.321-01 

7 

8-376-0! 

1.239*0? 

7.197-0! 

9. 599 -CM 

5.571-01 

* 

o. 1 70-o? 

*  .907-0  I 

8.938-^1 

5.763-01 

9 

1.000*00 

1.378*0? 

6.c«2-0l 

c  « 3 1 1  1 

5.97/  n] 

in 

1 .089*00 

1.372*0? 

6.927-01 

5 . 7 1 S  1 

6.I96-’! 

9.283*00 

5.113*00 

ruvror  CHTCK 

r*utr 

0160  79 

D< 

*.8?9-qi 

s.i'o-oi 

5.505-01 

5.89’-01 

*.3i'-oi 

6.768-0! 

7.1S?-ni 

7.507-01 

’.8K6-01 

B.738-01 

6.572*00 


T  09  O 

1.735-01 

1.873-01 
t • 20'-01 
9. 209-0? 
6.219-0? 
3.039-0? 
-7.958-03 
-9.98* -0? 
-9.T68-n? 
-1.80'-01 
3.  '88-'  . 


-2.106-03 

1.931-02 

1.897*00 

8.853-01 

31.07 

-1,878-03 

1.858-02 

1.899*00 

6.180-01 

39.89 

*1, *896-03 

1.778-02 

1,501*00 

5.539-01 

38.99 

-1.08  7-03, 

1.709-02 

1.509*00 

5.931-01 

91.86 

-6.037-08 

1.630-02 

l.S07*00 

6.368-01 

85.09 

-9.023-05 

1.557-02 

1.509*00 

6.820-0! 

88.13 

6.033-09 

I.S07-O? 

1 .877*00 

7.213-01 

50.62 

1.337-03 

1. <173-02 

1.827*00, 

7.567-ni 

S'*.  79 

7.189-03 

1.983-02 

1.777*00 

7.929-01 

69.70 

3.072-03 

1.917-02 

1.328*00 

8.299-OJ 

56.89 

-5.282-0* 

1.629-01 

3.987-01 

6.569*00 

T  r 

.20*99?nn*O3 

Z= 

.  U  ?M6’02*n’ 

TH= 

.7*850309*0’ 

0Hr  6  = 

tXz 

.85"58379T0? 

zr: 

.*68  *0791  .09 

THf  = 

-.18*0828’-08 

0*Gs 

TIG: 

-.36*20606*00 

?IG= 

.4*oj-iS02*on 

TMirr 

. 18  76351’-03 

OHGI  = 

r»r 

.$68*85*9*0? 

FHZ 

-.9>79rT60*0t 

rxz 

,  .*6’?089l*0l 

r  '= 

HTM: 

.O0oq9p91-Ol 

“U: 

-.52717*69-0*“ 

"YZ 

.929865S8-06 

mot 

VIN: 

•iocoobno*?i. 

»LT  = 

.38703?Oft«07 

ro: 

•ooonoooo 

GV= 

ITT: 

.90986«;o3*oi 

»«o: 

-.99873'! 

r.njz 

•  57  77689’*00 

PRT: 

- ,975051 *8*00  V  = 

-.88359583*0'  VI: 
.75X57696-02  Oz 

.55765987*0?  FZ 
.9*29  3719-08  0«U6  = 
.391867?0*07  OPS 
.789576*9*00  080  = 


.6736! 890*or 
-.89839285-01 
.11885000*02 
.55699861*07 
.26088  J«7*  02 
•ooooonon 

-.18778069-^3 


RUN  U 


t 

o 

VB 

w^S 

09 

L9 

O' 

T090 

1*9 

5*8 

CL 

FN 

ALPHA 

* 

3.7  32*01 

1.081*0? 

1  .39**00 

3  *  1 7  3  -n  | 

3.938-01 

8,702-01 

1.70?-01 

-? .1 5n-03 

2.999-02 

1.866*00 

8 . 765-01 

30.72 

? 

8. 16* -oi 

1  .06  3*0? 

1  .367*00 

3.317-01 

8.119-01 

*.019-01 

1.895-01 

-1. 867-1*3 

2.816-0? 

1.870*00 

5.077-ni 

79.89 

3 

6.nrxt-oi 

1. 0*9*0’ 

i.338*oo 

3.998  -01 

8.373-01 

5.36**01' 

1.197 

-01 

-1.518-03 

2,778-02 

1.878*00 

5.872-01 

3P.10 

8 

6  •  *  T9  -0 1 

1. 170*0’ 

1  .?97«nn 

3.703r^l 

8.592-01 

5.787-01 

9.082 

-02 

-1.099-0' 

2.6*2-02 

1.878*00 

5.821-01 

81.53 

5 

*. 666-01 

1.168*0? 

l  .269*0* 

3. *88  -Ol 

8.891-01 

C. 1*1-01 

6.119- 

-02 

-6.060-08 

2.552-02 

1.982*00 

6.289-01 

99.76 

6 

T. 502-01 

1.197*07 

1 .219*00 

9.217-01 

5. 270-01 

6.610-01 

2.99? 

-02 

-3.126-05 

2.882-02 

1.986*00 

6.701-01 

97.81 

7 

*.'76-01 

1.237*0? 

l  .1  ?9*nn 

9  .5?3-nl 

5.507-01 

7.012-01 

-6.68' 

-03 

6.212-06 

2.361-02 

1.962*00 

7.110*01 

50.38 

* 

». 170-01 

t.??6*0? 

1  .*  39*00 

9,860  -*  1 

S.696-oi 

7.356-01 

-8.830 

-02 

1.3S6-P3 

2*'10-02 

1  .,8 12*00 

7.960-01 

S’,  51 

9 

1 . onp.no 

1,322*0’ 

l  .099*09 

5.230-''! 

5.901-01 

7.709-01 

-9.218 

-02 

2.183-0' 

2*265-02 

1.367*00 

7.S17-01 

68.97 

P 

1  .0*8*00 

0 .3’0*0’ 

l  .061*00 

5.672-01 

6. 1 16-oj 

8.070-01 

-1.383 

-01 

3.108-03 

2.226-02 

1.715*00 

8.189-01 

66.27 

9.2'|9*oq 

6.031*no 

6.375*00 

3.289 

-01 

-6.106-0* 

2*588-01 

9.817-01 

6.960*09 

X 

=  .25188932*0! 

2=  .1*867.110*0' 

TH=  .75999818*02 

0  865s 

_ 

971921*6*00 

V: 

1 1099306*05 

XT 

=  .366*8755*.02 

?lr  .  ,*1638395*0’ 

THir  .19997276-09 

0MG  = 

- 

85025708*03 

WI=  - 

10187098*00 

*10 

=  -  .?*  13*8 ’2*0(1  ZtG  =  .91A59778*0''  THIT=  1  3R'3T37-03 

OHGir 

3671 7'22-0? 

Qz 

1 1885000*02 

r  a 

=  .6833O«*l*0? 

*M=  -.'2265263*01 

ftz  .*3750537*01 

F'= 

5683*537*02 

FZ 

55800897*0?  ' 

*CH 

=  .90*09091-01 

*V=  -.69777*66-05 

"T=  -.79635219-06 

*n: 

221879*6-08 

O’AGr 

25500218*02 

VIN 

=  .innononn*n!  «t  T=  .31 13?99«)*0' 

ro:  .ooopoooo 

GV= 

291979*3*02 

n*: 

none:*  non 

ATT 

=  .89S5 

?5’6*Oi  AHQz  -.21019951-0* 

C?T=  .57  786808 *00 

RPT  = 

76013691*00 

oo0:  - 

63921686-1 

08 

tJ 


$ 

I 

► 

? 


R 


l 


\ 

I 


£ 


RUN  15 


r 

0 

V9 

V" 

OR  . 

LR 

06 

TO’O  . 

MR 

6*6 

CL  ,FN 

ALPHA 

1  ! 

r.'Tj7-oi 

037*0’ 

I  ’.086*09 

7. 1 32  -0  1 

3.132-01 

8.172-01 

1.386-01 

-1.699-93 

2.960-0 2 

1.187*00  U. 202-01 

26.18 

i  2  « 

1.I66-OI 

t. 0*0*0’ 

1  .067*00 

'.786-91 

3.296-Ot 

8.397-01 

1.172-01 

-1.872-03 

2.883-02 

1,189-99  8.832-01 

29.58 

_s_  s.noo-oi 

_  1.071*0? 

1  .0  96*001 

3.385-ni 

3.385-01 

9.5*5-01 

9.525r02 

-1. 198-nj 

’.807-02 

1.191*00  8. 688-01 

3-.S6 

8  « 

i.878  -01 

1 .095*0? 

l‘.0?3*no 

3.599-"l 

I. 589-01 

9. 9J5-01 

’.255-0? 

-8.609-"* 

’.731-02 

1. 198*110  8.9-8-01 

36.01 

5  « 

;.*6f-oi 

1.1?S*0? 

9.975-m 

7.738-91 

3.778-91 

5.207-01 

8.882-0? 

-8  .'700-09 

2.656-02 

1.197*90  $.2>5-ni 

39.02 

6_  1 

'.602-01 

1.183*0? 

9.712-nj 

3*953-91 

3.968-01 

5.522-01 

2.387-0? 

-1.808-05 

’.581-02 

1.199*00  5.575-9J 

91.87 

7  l 

». '36-01 

1  ■  1*6*0? 

9.981-91 

8.193^1 

8.116-01 

6.897_01 

-5.789-03 

8.982-09 

7.528-02 

1.180*00  5.868-ni 

88.35 

8  « 

*.170-01 

1.227*02 

9.168-91 

8.958-01 

8.779-01 

6. 06R-O1 

-3.91 '-02 

1.070-03 

7.889-02 

1.196*90  6 . i ?7-01 

96-58 

9  1 

1.000*00 

1.260*0’ 

8.691  -01 

*.789-91 

9. 371-01 

6.3'5-Ot 

-7.920-0? 

1 .716-03 

7.897-02 

1.111*00  6.397-01 

98.57 

19  1 

1.088*00 

l.inn*o? 

8 •* 19-n  1 

5.069  -0 1 

8.509-9] 

6.510-01 

-1.111-01 

2.933-03 

2.9J9-02 

1.077*00  6.679-01 

99.87 

3.987*99 

3.8?8*0-j 

5.368*00 

2.669-01 

-5.610-07 

>.699-01 

8.273-01  5.817*00 

DIVIDE  CHFCK 

FAULT  A' 

016079  . 

, 

x: 

,  '12156 S6*0' 

2=  .29826139*07 

TH=  .765792864  ’ 

0H69S 

86152701*00 

N  = 

89567068*09 

xt  = 

.38055773*02 

71=  .10605711*0' 

TW'r  -.12251665-18 

0»6=  -, 

75 1*3'?8*03 

VIS  - 

10778978-9J 

*'<*- 

-  .18133' 1 1 *99  ?tG=  .99978?ll*np  ThII=  .17’69789-n8 

0B6ir 

95882105-07 

OS 

1  16850on*  o? 

F  A  r 

.83896978*0’ 

FH=  -.88529681*01 

FX:  .^3677579*91 

r  7-  , 

833*5777*92 

FZ 

93716570*0?  * 

hch: 

.90909091-01  I 

8U=  -.R^liACSB-O* 

MT=  .1 1957966-96 

**n=  , 

57625700-05 

0°  AG: 

21971033*02 

VIM: 

.lOOOOOOO*"'  »L7:  .20173867*0' 

CO:  .00000000 

6V=  , 

197'  89*1*0? 

OF: 

93000000 

A'T  = 

•  9 1 7*6*9  6*  0 1  AHO:  -.77786809-09 

fOT=  .88n99288*00 

RBT: 

79811083*00 

pro:  - 

29236208-08 

RUN  16 


T 

0 

V3 

v<5 

08 

l  * 

O' 

TORO 

MO 

668 

CL 

FN 

ALPHA 

1 

3. '3? -01 

1.095*02 

,  l  .182*90 

3.282-91 

8.639-01 

5.18’-01 

1.927-01  -2.358-0' 

2.256-02 

1.676*00 

<.2?'-Pl 

3'. 70 

? 

9. 166-01 

l.O'«VO? 

1.1  16*P*J 

3.8  08  -oi 

9.766-01 

5.56’-01 

1.639-01  -7.050-03 

2.198-02 

1.680«0Q 

5.615-01 

37.6* 

■? 

6  .nnfi-oi 

1. 0*9*0? 

1  .0*6*00 

3.60?  -O  ] 

S. 087-01 

5.5»’-0l 

1.3*0-01  -|.67'-03 

2.092-0? 

1  ,68'*no 

6.082-01 

81.39 

9 

8. *38 -01 

1.133*0? 

1.053*™ 

3.6(6-01 

5. 370-"1 

6.992-01 

1.026-01  -1.217-03 

1.937-02 

1.687*00 

6.510-81 

88.9* 

1 

* 

6.668 -01 

1.171*0’ 

l  .OJ9*00 

8 .105 -9 1 

5.78  7-01 

6.9* '-01 

6.915-02  -6.760-0* 

1.838-021 

l.69i*on 

7.070-01 

88i?7 

6 

',502-01 

t. 217*07 

9.8RC-01 

8.811-91 

6.176-01 

7.896-1)1 

3.331-02  -8.5*5-05 

1.733-02 

1.695*00 

7. 572-01 

SI  .39 

7 

*.'36-01 

,3  .767*07 

<).«R*-ni 

9.763-01 

6.508-01 

7.958-01 

-8.87T-03 

6.78 7-08 

1.662-02 

1.662*00 

8.0*8-05 

S'. 92 

A 

R. 170-0| 

1.305*3’ 

9.199-01 

S. 131 -91 

6.755-01 

8.3'8-01 

-5.S06-C2 

1 .890-23 

1.619-02 

1.601  400 

8.*59-rtI 

55.98 

9 

1 .900*00 

1.3*5*07 

8.810-01 

5.599 -ill 

7.015-91 

8.766-01 

- 1 .03 I -OI 

2. 8IO-03 

1.581-02  ‘ 

l.S8?*00 

8  •869-,'l 

57.86 

1° 

t ,0*8*09 

1 .808*0’ 

8  .8  *1  -0 1 

5.9*8  -^1 

8  •9'1?  *93 

'.7R7-01 

S.920*0?i 

R.I88-QI 

7.185*00 

-1.567-01  3.838-03 

3.723-01  -1.357-06 

1.598-02 

1.836-01 

1.888*00 

9.182-01 

9. 299-81 

7 • 265*00 

59.57 

< 


'  = 

.233)0691*03 

Z*. 

.178?'6R9*Q' 

Tm: 

.776073^6*0? 

OH'*: 

-.in28?1?l*01 

VS 

.811 970'9*  09 

xrr 

•2*9396 1 2*02 

71  = 

.79*89667*0’ 

THI  = 

.13888810-08 

0**Gr 

-.88359910*0' 

vis 

-.21629372-01 

Tir-r 

-.26*07119*00 

2TG: 

*  3'33’778 

y»rr: 

-.*56*0381-08 

0H6T  = 

'.1?09R9*3-P2 

OS 

.11885000*02 

FA  : 

.61*58*63*9? 

fn: 

-.6?3o»R96*ni 

F=: 

.11858998*01 

F'  = 

.61 735' 1 2*o’ 

Fz 

•6?l5207?*02 

HCH  : 

.90*99051-01 

mu: 

-. 1'69?1 16-0' 

mt: 

-.?6’?0$39-O6 

MO; 

.73818725-05 

0°  AGs 

,2  87*19*9*02 

vin: 

.lP909"nn-n' 

»L'  = 

. 3? 1 '6'96*0* 

'*0: 

.onnnaooo 

Gvr 

.21896659*0’ 

of: 

."ooaooon 

ATT: 

.56’97?97*91 

•  mo: 

-.537R31*?-9«; 

COT: 

.63781186*00 

pRT; 

.7579  7’nj.00 

0  no: 

-.19889778-08 

F 

i 
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1  T*5*  Wwf'MiW-xA:  f'^- 


fjafSTH^vrwv  *-t 


RUN  17 

I 

ft 

%' 

W<6 

08 

18 

08 

TO*e  88 

578 

Ct 

FM 

ALPHA 

t 

3.777-01 

!. 051*07 

1  .200*"9 

7.715-0? 

3.208-01 

8.508-01 

1.857-01  -3.0*8-03 

1.535-02 

l. 507*00 

8.635-01 

38.66 

7 

8.t66-m 

1.085*0? 

1  .1*0*90 

7.882-01 

•8.887-01 

<.087-01 

1.521-01  -2.711-03 

1.867-02 

1 .510*00 

5.078-01 

33.08 

7 

6.ono-ni 

1.173*0? 

l  .i?8*np 

3.108-0? 

8. '17-OJ 

5.528-01 

1.38-*-0I  -2.233-03 

1.330-0? 

1.513*00 

5.563-01 

83.23 

8 

5.838-01 

1.170*02 

l.0si«rtn 

i.m-oi 

3.71 *-ni 

5.053-01 

1.035-01  -1.635-03 

1.713-02 

!.515*"0 

6.110-01 

37.08 

3 

6.661-01 

1.228*0? 

1.0  3  ft*00 

3.G58-OJ 

5.570-01 

5.583-01 

7.235-0*  -3.381-03 

1.280-0? 

1.516*00 

6.703-01 

50.68 

t 

▼.502-01 

*.281*9? 

4.90-01 

7.99?-9| 

«. 187-01 

7.287-01 

3.718-02  -1.021-03 

1.167-02 

1.521*00 

7,333-01 

57.31 

7 

*.73f-01 

*.37«*o? 

9.616-01 

8.3*S~pi 

5.552-01 

7.788-0! 

-5.288-03  i. 651-08 

1.125-02 

1. 375*00 

7.857-01 

56,83 

0 

*  •  1 70  -0 1 

I. *00*02 

S.t^-"2 

8. 77* -*3 

5.888-01 

8.255-01 

-5.853-02  1. 078-03 

1.038-02 

1.311*00 

8 • 327-01 

S8.51 

« 

'.000*00 

t. 8*3*0* 

j  .70*- -01 

*•  ?3fl  Ol 

7.161-01 

8.721-0! 

-l.OSp-Ol  3.232-03 

1,075-0? 

1.333*00 

8.738-01 

60.33 

13 

i  .OM«no 

1.530*0? 

8.370-01 

5.721  1 

7.885-01 

3.187-01 

-1.501 -01  8.687-03 

1.055-02 

1.231*00 

3.277-01 

52.08 

3.882*00 

5.*5?*00 

5.315*00 

3.728-01  -1.363-06 

1.287-01 

2.361-01 

6.371*00 

t- 

.26033120*03 

2= 

.18l6S«02*0T 

TMr 

.73867113*07 

e*K68r 

-.12013177*01 

v: 

.37327303*08 

*t  = 

.37366660*02 

?T  = 

• 77ni 3576*02 

TMT  = 

.130?76t"-08 

0**Gr 

-.10353570*08 

UIZ 

-•21578809-i'j 

xir,; 

-.26708513*00 

ZIC: 

.36757163*00 

THTir 

.€7505026-05 

ewi: 

.13883027-0? 

Qz 

•11895000*^2 

rA= 

.83650760*0? 

FM  = 

5761*01 

Fxz 

•  63 1 37067  *0! 

F*  r 

.83750552*1 12 

fz 

•880?65S7*02 

MCHr 

.30308031-01 

-.17653305-06 

MTZ 

.*8787888-07 

. 3130’881-06 

0*  A  S3 

.276588?«*02 

*!H  = 

.10000000*03 

ALY  = 

.3! *01 038*0' 

CC= 

.00000000 

Gvr 

.27738183*0? 

OF  r 

.oonoooon 

ATT  r 

.5-*66S3?l*01 

AtfQ- 

-.73002832-05 

corr 

,60370280*00 

ft*Tr 

.76086100*00 

ft  *03 

-.13871833-03 

RUN  18  '  > 


1 

0 

V" 

V*6 

08 

L* 

0< 

T890  •• 

666 

CL 

6* 

ALPHA 

1 

7.332-01 

l.O?3*0? 

-  .061*00 

3.657-8  1 

3.779-01 

8.360-01 

1.561-01  -1.5*7-03 

3.609-02 

1.887*00 

5.011 -nl 

?7.88 

7 

3.166-01 

1 .038*0? 

i  .G85*"0 

7.776-01 

3. 307-01 

6.200-01 

1.318-01  -1.338-03 

3.697-0? 

1.886*00 

5.256-01 

71.06 

3 

6. poo-oi 

1.062*02 

1  .027*0* 

3.322-01 

8 .052-01 

5.862-0! 

1.061-01  -1.075-03 

3.785-02 

1.850*00 

5.525-01 

78.19 

8 

5.033-01 

1.088*02 

1.005*00 

8.038-01 

8.270-01 

5.787-01 

1.005-0*  -7.631-08 

3.878-0? 

1.858*00 

5.817-01 

57.20 

5.668-01 

1.108*0* 

3.M3-H1 

8.232-01 

8.876-01 

6.066-Ot 

5.218-0?  -8.176-08 

3.763-0? 

1.868*00 

6.173-01 

90.09 

6 

7.507-01 

1.175*0? 

3.603-0! 

8.517^01 

8.668-01 

6.788-0! 

2.838-0?  -3.770-06 

3.252-02 

1.862*00 

6.878-OJ 

81*. 06 

7 

3.776-01 

1.168*0? 

3.36/-P1 

8.7G3-01 

8.873-0? 

6. 7M-01 

-7.306-03  8.577-08 

3. 157-02 

1.052*09 

6.803-nl 

85.79 

* 

6.170-01 

1.136*07 

3.1  18-01 

5.087-01 

5.011-0! 

6.5*2-01 

-8.378-02  9.687-08 

7.096-0? 

1.817*00 

7.079-01 

87.81 

3 

1  .000*00 

!. 270*0? 

8.868-01 

5.751-OJ 

5.167-0! 

7.262-01 

-6.186-02  1.580-03 

3.081-02 

1.375*00 

T. 362-01 

89.75 

10 

1.083*00 

1.266*0? 

8.618-01 

5.633*01 

5.302-01 

7.583-01 

-1.218-01  2.172-03 

2.991-02 

1.336*00 

7.653-01 

St. 16 

3.511*80 

8.581*00 

6.27?*00 

2.967-01  -8.261-07 

3. 787-01 

5.750-01 

6.311*00 

xr 

.?8698I8**07 

Zz 

.27933761  *0-* 

'THr 

.78916721*02 

0H76T 

-.80577870*00 

V? 

•  8  33256 16*  n8 

xir 

.30768162*0? 

71  = 

.96086067*0? 

THTr 

-.87058011-05 

0**G- 

-.70997X27*03 

HI: 

-.76108360-9? 

XTOr 

-.lP7»61ft2*00 

?JG  = 

.8n5669?8*00 

THfTS 

.33030386-10 

omrz 

.9507*1*9-03 

0: 

.11805000*9? 

Fir 

.38856866*0? 

?«: 

-.87250877,01 

FX  = 

.62*16393*01 

f-»z 

.38198733*0? 

Tz 

.30699095*0? 

xTH; 

,9O«0?Q91-Ot 

»;y  r 

-%88G60Y*0-O^ 

PTr 

.6*.M 83810-!* 

.56*68808-06 

0*A6- 

.28926550*92 

V  TMr 

,innooonn*07 

ALT- 

.2?:.36?79*0’ 

ro: 

.09000000 

6Vr 

.1789t835*0? 

OF- 

."0000009 

ATT  z 

,8«88?660*01 

A  HO” 

-.*2628638-05 

rotr 

.56935089*00 

ftBTr 

.78 1268 *5*00 

9P0S 

-.28513703-08 

not  reproducible 
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APPENDIX  D 

COMPUTER  SOLUTION  PRINTOUT  OF  THE  DEPLOYMENT  TRANSIENT 
FOR  THE  NOMINAL  CASE-TABLE  III,  RUN  1 


mom  ti« 

.ownooni)  .inonioon.o! 

.inonnoon.ni 


.32900ogo«oo 

8«IS= 

.*0000000*"’ 

£LTO= 

.50000090.07 

vi*rr 

.10000000*03 

•wor 

.27770009-02 

S: 

.708*oooo.og 

Cr 

.75000000.00 

1  = 

♦***0"0"o.no 

•  r 

. 70309 n"o.on 

GH? 

."oooono" 

*M= 

.26000000-02 

11  = 

.$8680000-0* 

to: 

.57*ooo«o-n? 

*vr 

.ooo"no"0 

"V  = 

."000000" 

CL5S= 

.15500000.01 

cns*  = 

.12500000*01 

TY5  5= 

*575000^0*00 

o«c= 

-.75000000*03 

0*61  = 

•ooooooon 

or*; 

.oooonooo 

rnvi  = 

.Oooonono 

C*V2= 

. *0000000 

cmrz 

-.T50"09"0.0" 

Conor 

.0000000" 

.2oonoooo«no 

CM=. 

.1  5000000."" 

C"?=» 

-.*5000000*00 

CLI  =  * 

•i7i**nno.ot 

CL  2=  • 

-.?i*jooo"*oo 

THOr 

.inoonooo.oi 

THir 

.oooooooo 

v: 

•oooooono 

81  = 

•0O0oo"P0 

rr^-z 

.nooooocr 

*MTT= 

.OOOQonoO 

vimz 

.1n09O0"0.'* 

•  I**; 

.*3320000*00 

08  = 

.8**00000-01 

in; 

.•0000000*01 

i«r 

.1O0O00«0*02 

It  r 

10 

.00* 

I  *  *•  u<<  ok  i*  0?  fi»«  “>  ?r«  ri  nt  «Ipm» 


1 

*.332-01 

l.oon*o? 

i .357-na 

1.5*5 -o* 

•.939-0* 

1.5*5 -0* 

-2.091-0* 

7.51 *-05 

7.323-0* 

3.505-9* 

3.577-"* 

.7? 

2 

3.166-0! 

1.000*0? 

1  .357-0* 

1  ,G*7-o* 

8.939-05 

1.6* 1-0* 

-2.732-05 

7. ■•*-05 

3.553-0* 

3.G05-O* 

3.689-"* 

.*1 

7 

6.000-01 

1.000*0? 

1 .75?-0* 

1.713-0* 

8.9*9-96 

1.70*-03 

-S.**5-o$ 

3.125-05 

3.705-0* 

3.505-03 

3.M1-0* 

.50 

• 

6  .ft**  -0 1 

1.000*0? 

1.357-0* 

1.755-0* 

8.939-05 

1 .7*7-0* 

-••303-05 

a  8*51-05 

3.aa*-o* 

7.505-0* 

3.339-0* 

.58 

6 

G.saa-oi 

1.000*0? 

1  .357-0* 

l.«*l-o* 

8.97*-06 

i.a7*-o* 

-1.133-0* 

a. a*i-05 

3.031-0* 

3.505-0* 

3.977-n* 

.55 

c 

t.*0?-91 

1.000*0? 

l  .757-0* 

MM-"* 

8.9*6-"* 

1  .5*0-0* 

-1.357-0* 

3.275-05 

3.325-0* 

3.505-0* 

*.220-0* 

.7* 

7 

■•*35-01 

1.000*0? 

l  .357-0* 

7.115-0* 

i .535-05 

2.007-0* 

-I. *00-0* 

3.753-05 

9.586-9* 

3.595-0* 

3. 355-0* 

.8? 

a 

■  •170-01 

i.noo*o? 

1  .357-0* 

7.252-0* 

a. 335-05 

2.21 

-2.153-0* 

l.02*-05 

*.•75-0* 

3.505-0* 

3.523-0* 

.91 

• 

*. 000*00 

1.000*0? 

1  .3*7-0* 

2.398 -p* 

a. 375-05 

2.755-0* 

-2.5**-0* 

1.0**-05 

5.131-0* 

3.505-0* 

3.582-0* 

.99 

I" 

i .n**.oo 

1.000*0? 

1 .357-0* 

2.557-o* 

a.375-05 

2.S02-*** 

-2.357-0* 

1.152-05 

5.533-9* 

3.505-0* 

*.8*3-0* 

1.97 

l.W-o* 

1  .375-0* 

1.976-ft* 

-1.*?*-03 

3.251-05 

* -317-03 

1. 833-07 

3.173-07 

nwmr  CMfCn 

.‘*«LT  at 

01  GO?* 

t: 

.******t*.oo 

2= 

.»G*a6*17-0? 

TMr 

•13167*1* .*1 

0"<S  = 

-•20?G* 197-01 

nz 

.!33*70?*."J 

i*  = 

.99*99*56.0? 

?T  = 

.3?1739*2»0" 

THT; 

• 1 l*3919?.0l 

0*G= 

-.7*«52*oi.07 

8T  = 

•31225587.02 

*I«= 

- .1*7*8171-03 

216= 

.9*9*9**340'' 

thii= 

.  1 3?6??7‘».03 

0«GT  = 

.77096097.0? 

as 

.lli*«0on.n2 

FA; 

.5509*769-0! 

ntz 

-.1017573Y.01 

rtz 

.11612231-01 

f?r 

.191810*9.01 

r; 

•19I817F1.01 
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5. *73-0! 

1.156-01 

-3.357-03 

2.273-02 

1.387*00 

S.523-"! 

38.63 

« 

5. *36 -01 

1.1*5*02 

1  .0  71  *0f» 

3.770-01 

*.675-01 

5. *75-01 

•.713-0? 

-2.600-03 

2. 113-02 

1.331*00 

5.332-01 

37.16 

5 

«.668  -01 

l.J60«0* 

1  .033*On 

*.021 -"l 

*.3*7-0! 

6.3"6-0l 

5.675-0? 

-1.363-03 

2.013-02 

1.333*00 

6.371-01 

35.33 

6 

7.502-01 

1.136*0* 

1  •005*0'! 

*.307-01 

5.3*0-01 

6.773-0! 

2.*?* -0? 

S. 736-05 

1.935-02 

1 .338*00 

6.8*6-01 

38.50 

7 

•  .336-01 

1 .2*2*0* 

3.707-01 

*.6*6-oj 

5.736-01 

7.275-01 

-1.067-0? 

1 .735-03 

1.302-02 

1.502*00 

7.357-01 

51.36 

8 

•.170-01 

1.2*7*02 

3.367-01 

7.378-oj 

6.7*1-01 

7.3*7-01 

-3.3*3-02 

3.603-03 

1*776-02 

1.536*00 

8.038-01 

53.35 

3 

t.  "00*00 

1 • 37**0* 

3.033-01 

5.36* -"I 

6.710-01 

8.3*1-01 

-t. 735-02 

5.633-03 

1.703-02 

1.521*00 

8.573-ni 

56.66 

I" 

1 .088*00 

1 .3***0? 

8.708-01 

5.78* -oj 

7 .038-01 

•*37*-01 

-1.635-01 

7.830-03 

1.6*7-02 

1.335*00 

3.0*6-01 

58.63 

-.268*00 

5.3*5*00 

6.6*2*00 

3.217-01 

3.620-06 

2.033-01 

3.237-01 

fc. 763*00 

xr 

.23168251*07 

2  = 

.11708665*0’ 

THr 

.7731506  8  *  0* 

0H<5r 

-.397173  66*00 

Mr 

.920685*6* "* 

xir 

•5*?*6*8**02 

71  = 

•  6850OR5* *0’ 

TMTr 

-.67081113-03 

0*»0r 

-•863"89?9*03 

Ml  r 

•12816067«"0 

Xinr 

“-•*05  8  6  3  06*  00 

7TG= 

.*67319*8*0" 

TMTTr 

. 17*81633-01 

OHGtr 

•73532612-01 

•  r 

.  i  18850""*"2 

Ftr 

.57630*66*0? 

Fur 

-.61  25  7  011*01 

rxr 

•*63?fl*J**0l 

r  t= 

•57627516*0? 

Fr 

.5*0|3773*02 

■CHr 

.30609031-01 

■Mr 

.7520*682-05 

•Tr 

.967703*6-0* 

•or 

.3379*163-03 

0*A$r 

.26768175*0? 

VIN  = 

.10000000*03 

ALTr 

.7"261055*07 

ror 

.00000000 

6vr 

.373??865*0* 

OFr 

.0000000" 

A’tr 

.52"27658*01 

AMOr 

-.5*337018-0’ 

rotr 

.56535732*00 

86Tr 

.7571*8*9*00 

••Or 

-.16913723-02 

not  reproducible 
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